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INTERACTIVE FORTRAN IV COMPUTER PROGRAMS FOR THE THERMODYNAMIC 
AND TRANSPORT PROPERTIES OF SELECTED CRYOGENS [FLUIDS PACK]* 

R. D. McCarty 

Thermophysical Properties Division 
National Engineering Laboratory 
National Bureau of StandarJs 
Boulder, Colorado 80303 

The thermodynamic and transport properties of selected cryogens have been 
programmed Into a series of computer routines. Input variables are any two of P, 
p or T In the single phase regions and either P or T for the saturated liquid or 
vapor state. The output Is pressure, density, temperature, entropy, enthalpy for 
all of the fluids and In most cases specific heat capacity and speed of sound. 
Viscosity and thermal conductivity are also given for most of the fluids. The 
programs are designed for access by remote terminal; however, they have been 
written In a modular form to allow the user to select either specific fluids or 
Specific properties for particular needs. 

The program Includes properties for hydrogen, helium, neon, nitrogen, 
oxygen, argon, and methane. The programs include properties for gaseous and 
liquid states usually from the triple point to some upper limit of pressure and 
temperature which varies from fluid to fluid. Computer listings of the 
FORTRAN IV codings are presented. Copies of the programs may be obtained from 
either the Thermophysical Properties Division of the National Bureau of Standards 
at Boulder, Colorado, or from Walter Scott at the NASA-Johnson Space Center In 
Houston, Texas. 

Key Words: Argon; computer programs; density; enthalpy; entropy; equation of 

state; heat capacity at constant pressure; heat capacity at constant 
volume; helium; hydrogen; methane; neon; nitrogen; oxygen; thermal 
conductivity; velocity of sound; viscosity. 

1. Introduction 

Recent technological developments which made possible personal Interaction 
with large computer systems by remote terminal has proved to be a real advantage 
In situations where Immediate response to requests for data are Important. This 
computer program was developed In response to a need for thermodynamic and 
transport property data In an Interactive mode. The program makes possible the 

* Work sponsored by NASA-Johnson Space Center on Contract No. P0 T-6570C. 


acquisition of thermodynamic and In some cases the transport properties of seven 
different cryogenic fluids. Those fluids are: hydrogen, helium, neon, nitrogen, 
oxygen, argon, and methane. All of the properties are calculated from empirical 
equations which are derived from the existing experimental data via weighted 
least squares fit of mathematical models to those data. The reader who Is 
Interested In the details of the correlations Is referred to the sources cited 
for each fluid In the sections to follow. 

Five of the fluids (hydrogen, nitrogen, oxygen, argon and methane) utilize 
the same mathematical model of the equation of state, and four of the fluids 
(nitrogen, oxygen, argon and methane) utilize the same correlating function for 
the transport properties. No transport properties for hydrogen, or neon, are 
given. Using the same functional form for the properties of several different 
fluids greatly simplifies the computer programs, since all that Is needed to 

switch from one fluid to another Is a different set of coefficients to the 
equations. 

The properties of helium and neon are calculated from empirical functions 
which are different than those of the other five as well as each other. The 
addition of these two fluids represents approximately half of the total length 
(core storage) of the entire program. 

The programs have been written with a dual purpose in mind. As has already 
been mentioned the primary purpose Is to make available to the user In an 
Interactive mode, the thermodynamic and transport properties of selected fluids. 
The second purpose Is to provide the necessary subroutines to the user who has a 
very specialized need. Therefore the program has been written In a modular form 
which allows the user to extract almost any combination of fluids and 
properties. 

The remainder of this report includes a section devoted to each fluid and an 
additional section describing the general part of the program. 

2. Uncertainty of the Calculated Properties 

When calculating thermodynamic properties from an empirical equation of 
state, one should be aware of certain problem areas where It Is difficult to 
estimate realistically what the uncertainties really are In a given property, in 
the critical region (p c ± 0.2 p c and T c ± 0.05 T c ), the calculated density may 
be In error by several percent, while the calculated pressure or temperature will 
not be as Inaccurate. Specific heat capacities and the thermal conductivity In 
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the critical region become very large values and no realistic estimate of the 
uncertainty may be made. Saturation boundaries, gas, liquid and solid are 
potential areas of large uncertainties for derived thermodynamic oropertles. 
especially heat capacities. In the compressed liquid, calculated pressures will 
have an uncertainty of several percent. This Is a consequence of the nature of 
the surface and Is In no way the fault of the equation of state. 

In each of the following sections on the Individual fluids, uncertainty 
estimates will be given for that particular fluid. These estimates have been 
obtained from the source documents and do not reflect the potential large 
uncertainties of the problem areas outlined above. 

3. Computer Routines 

The thermodynamic properties produced by the computer program listed In 
Appendix A are all calculated using a mathematical model of the equation of state 
of the fluid and classical thermodynamic relationships. The reader who Is 
interested in the thermodynamics and mathematics of the problem is referred to 
McCarty (1975). Table 1 lists all of the subprograms with pertinent Information 
for each. 

In general eleven or twelve significant figures need to be carried in the 
property calculations to insure no round off errors in the result, and sample 
values for checking performance are included In each of the sections on the 
properties of a particular fluid. Performance should be checked perlodicrlly 
when running on a time- share system. 

The program, once executed, leads the user through a series of Input 
requests and the use of the program should be self-explanatory after the 
particular system requirements for access and execution are satisfied. 

Upper and lower pressure and temperature limits have been Imposed on each 
fluid according to the range of validity as claimed by the original source 
document. The lower pressure limit is an Indication of where the program fails 
due to the limitations of the computer calculations. All of the functions reduce 
to the Ideal gas in the limit of zero pressure but because of some terms In the 
equation which have very large exponents, the Iterative solution for density 
falls at very low pressures. Because these models are empirical and cannot be 
trusted to give even reasonable results outside their ranges of pressure and 
temperature, the user Is cautioned not to change the pressure and temperature 
limits originally set In the program. 
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4. Conclusions 

The preparation of this program has demonstrated the utility of maintaining 
functional form when modeling properties of several different fluids. With a 
single functional form, the change from one fluid to another becomes a matter of 
switching constants and even more Important the addition of more fluids to the 
scope of the program becomes very simple. 

The author wishes to acknowledge the support of NASA, without which this 
work would not have been possible, and In particular the support and 
encouragement of Walter Scott at the Johnson Space Center. 
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FUNCTION CV Calculates the specific heat capacity 

of constant volume of a fluid D,T 

FUNCTION FINDD Solves the equation of state for a 

density given pressure and temperature P,T 
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of a fluid D,T D,T A nW/ro-K 

FUNCTION EXCESV Calculates the "excess" viscosity of 

« fluid D,g/cm 3 ,T D,g/cm 3 ,T An **g/cm-s 
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derivative of the second vlrlal for 
helium 

FUNCTION VPNHE Calculates the vapor pressure of helium 

FUNCTION CPHE Calculates specific heat capacity at 

constant pressure for helium 
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temperatures below 300 K D,T D,T n rolcropoise 
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Table 1. Continued 
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FUNCTION ENTHNE Calculates enthalpy for neon - gaseous 

phase P»D,T P»D,T 

FUNCTION ENTRNE Calculate entropy for neon in the 

gaseous phase D»T 0,T 



5. Properties of Hydrogen 

The thermodynamic properties of hydrogen are calculated from a 32 term 
empirical equation of state. The source of the equation of state Is NBSIR 75-814 
(Roder and McCarty, 1975). The program allows the option of ortho, para, and 
normal hydrogen modifications. The properties for normal and ortho hydrogen are 
obtained by changing the Ideal gas thermodynamic properties from parahydrogen to 
ortho or normal hydrogen. This changes the derived thermodynamic properties but 
does nothing to the PVT surface. Therefore the PVT for all of the hydrogen 
options are the same. 
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Table 2. Properties and Associated Estimated Uncertainties for Hydrogen 


Property 

All Uncertainties In Percent (Except Enthalpy) 

Liquid (below 1 Q ) Gas (below T ) Fluid (above T ) 

^ c 

pressure 

5.0 

0.25 

0.20 

density 

0.1 

0.25 

0.2 

temperature 

0.1 

0.25 

0.2 

enthal py 

1 J/mol 

3 J/mol 

5.10 J/mol 

entropy 

1.0 

1.0 

1 

C P 

3.0 

2.0 

3.0 

Cv 

3.0 

2.0 

3.0 

speed of sound 

2.0 

1.0 

1.0 


The above are based on the estimated average difference between calculated 
and true value. Including problem areas as described In section 2. 
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Table 3. Pressure and Temperature Range for Hydrogen 


Temperature range 
Pressure range 

Fixed Points 

Triple Point Temperature 
Triple Point Pressure 
Triple Point Density-Vapor 

Liquid 

Critical Point Temperature 
Critical Point Pressure 
Critical Point Density 


13.8 to 400 K 
-434.83 to 260.43°F 
0** to 1200 bar 
0 to 17404.5 psla 

13.8 K* 

-434.83°F 

0.07043101 bar** 
1.0215154 psla 
0.0632229 moles/liter 
0.007956667 Ib/cu ft 
.21429 moles/liter 
4.809273 lb/cu ft 
32.938 K* 

-459.67°F 
12.83768 bar 
186.19481 psla 
15.556 moles/liter* 
1.95774 lb/cu ft 


*TH1s particular value is used to define the point In the sense that it was 
used as input Into the program. 

**See Section 3 for more detailed explanation of the lower limits of 
pressure. 


Table 4. Properties for Computer Program Performance Test for Hydrogen 
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6. Properties of Helium 

The thermodynamic properties of helium are calculated from a 32 term 
modified Benedlet-Webb-Rubln equation of state. The source of the equation of 
state Is a paper by McCarty (1973). The transport, properties for helium are from 
NBS Technical Note 631 (McCarty, 1972). Tabular values of the thermodynamic and 
transport properties of helium are given In Technical Note 'i31 and are identical 
to those produced by this computer program. 
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Table 5. Properties and Associated Estimated Uncertainties for Helium 


* 


Property 

All Uncertainties In Percent (Except Enthalpy) 

Liquid (be’ow T £ ) Gas (below T c ) Fluid (above T £ ) 

pressure 

10.0 

0.2 

0.2 

density 

0.5 

0.2 

0.2 

temperature 

.5 

0.2 

0.2 

enthal py 

2.0 

1.0 

1.0 

entropy 

2.0 

1.0 

1.0 

C P 

2.0 

2.0 

.5 

C v 

2.0 

1.5 

.5 

speed of sound 

2.0 

1.0 

.5 

thermal conductivity 

10.0 

10.0 

10.0 

viscosity 

8.0 

8.0 

10.0 


The above are based on the estimated average difference between calculated 
and true value. Including problem areas as described In section 2. 
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Table 6. Pressure and Temperature Range for Helium 


Temperature range 
Pressure range 


2 to 1500 K 
-456.07 to 2240. 33 # F 
0** to 1000 bar 
0 to 14504. psia 


Fixed Points 

Lambda Temperature 

Lambda Pressure 

Lambda Density-Vapor 

Liquid 

Critical Point Temperature 
Critical Point Pressure 
Critical Point Density 


2.1720 K* 

-455.760°F 

.04963285 bar** 
.71986364 psia 
.2904484 moles/liter 
.72575547 ib/cu ft 
36.53426 moles/ 11 ter 
9.1289672 Ib/cu ft 
5.2014 K* 

-450.307°F 
2.274640 bar 
32.990864 psia 
17.3987 moles/liter* 
4.347485 lb/cu ft 


*Th1s particular value Is used to define the point In the sense that It was 
used as Input Into the program, the program does not reproduce this number 

exactly. 

**See Section 3 for more detailed explanation of the lower limits of 
pressure. 
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Table 7. Properties for Computer Program Performance Test for Helium 




7. Properties of Neon 

The thermodynamic properties of neon are calculated from an 18 term 
empirical equation of state by McCarty and Stewart (1965). No specific heat 
capacities are calculated for neon because the accuracy of the model of the 
equation of state does not warrant the calculation. Transport properties for 
neon are not available In the form of a mathematical model as Is tK case for 
some of the other fluids. 



1 
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Table 8. Properties and Associated Estimated Uncertainties for Neon 


Property 

All Uncertainties In Percent (Except Enthalpy) 

Liquid (below T £ ) Gas (below T c ) Fluid (above T c ) 

pressure 

10.0 

0.5 

0.5 

density 

2.0 

1.0 

1.0 

temperature 

2.0 

1.0 

1.0 

enfal py 

20.0 

1.0 

10.0 

entropy 

20.0 

1.0 

10.0 


The above are based on the estimated average difference between calculated 
and true value. Including problem areas as described In section 2. 
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Table 9. Pressure 


and Temperature Range for Neon 


Temperature range 
Pressure range 

Fixed Points 

Triple Point Temperature 
Triple Point Pressure 
Triple Point Density-Vapor 

Liquid 

Critical Point Temperature 
Critical Point Pressure 
Critical Point Density 


28 to 300 K 
•414*67 to 80*33 F 

. 1 ** to 200 bar 
1.450 to 2901. psia 

25 K* 

-414.67°F 

.5102339 bar** 
7.4003172 psia 
.2528246 moles/liter 
.31855482 lb/cu ft 
61.45071 moles/liter 
77.42688 lb/cu ft 
44.40 K* 

-379.75'F 
26.57086 bar 
385.37775 psia 
23.93 moles/liter* 
30.151 lb/cu ft 


Ms particular ,.1ua Is used to define the point In the sense that It -as 

ed as input into the program. „ 4 . . 

See Section 3 for more detailed explanation of the lower limits 


pressure. 
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Table 10. Properties for Computer Program Performance Test for Neon 



8. Properties of Nitrogen 

The thermodynamic properties are calculated from a 32 term modified equation 
of state taken from Jacobsen, et al. (1973). The same functional form of the 
equation of state used for nitrogen ha also been used here for hydrogen, oxygen, 
argon and methane The transport properties are from Hanley, et al. (1974). 
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Table 11. Properties and Associated Estimated Uncertainties for Nitrogen 


Property 

All Uncertainties In Percent (Except Enthalpy) 

Liquid (below T c ) Gas (below T c ) Fluid (above T c ) 

pressure 

5.0 

0.3 

0.3 

density 

0.5 

0.3 

0.2 

temperature 

0.5 

0.3 

0.2 

enthal py 

3.0 

1.0 

1.0 

entropy 

2.0 

1.0 

1.0 

c p 

5.0 

5.0 

5.0 

C v 

5.0 

5.0 

5.0 

speed of sound 

2.0 

0.25 

1.0 

thermal conductivity 

4.0 

4.0 

6.0 

viscosity 

2.0 

2.0 

2.0 


The above are based on the estimated average difference between calculated 
and true value. Including problem areas as described In section 2. 
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Table 12. Pressure and Temperature Range for Nitrogen 


Temperature range 
Pressure range 


63.15 to 1900 K 
346 to 2960. 33°F 
0** to 1000 bar 
0 to 145040. psla 


Fixed Points 

Triple Point Temperature 
Triple Point Pressure 
Triple Point Density- Vapor 

Liquid 

Critical Point Temperature 
Critical Point Pressure 
Critical Point Density 


63.15 K* 

346 # F 

.1246399 bar** 
1.807749 psla 
.2396164 moles/ 11 ter 
.04190847 Ib/cu ft 
30.S7717 moles/liter 
54.178504 Ib/cu ft 
126.26 K* 

-232.402 # F 
34.10034 bar 
484.5836 psla 
11.21 moles/liter* 
19.606 Ib/cu ft 


*Th1$ particular value Is used to define the point In the sense that It was 
used as Input Into the program. 

See Section 3 for more detailed explanation of the lower limits of 


pressure. 


Table 13. Properties for Computer Program Performance Test for Nitrogen 



9. Properties of Oxygen 

The thermodynamic properties of oxygen are calculated from a 32 term 
mplrlcal equation of state given In NBSIR 78-882 by Weber (1978). The 
functional fonn of the oxygen equation of state Is the some as was used for 
hydrogen, nitrogen, oxygen and methane. The transport properties are from Hanley 

et al. (1974). 


Table 14. Propert'es and Associated Estimated Uncertainties for Oxygen 


Property 

All Uncertainties 
Liquid (below T„) 

w 

In Percent (Except Enthalpy) 

Gas (below T c ) Fluid (above T c ) 

pressure 

5.0 

0.25 

0.15 

density 

0.1 

0.25 

0.15 

temperature 

0.1 

0.2 

0.1 

enthal py 

0.5 

0.25 

0.5 

entropy 

0.5 

0.25 

0.5 

c p 

3.0 

5.0 

3.0 

Cv 

3.0 

5.0 

3.0 

speed of sound 

2.0 

0.5 

0.5 

thermal conductivity 

4.0 

4.0 

6.0 

viscosity 

2.0 

2.0 

2.0 


The above are based on the estimated average difference between calculated 
and true value. Including problem areas as described in section 2. 


Table 15. Pressure and Temperature Range for Oxygen 


Temperature range 
Pressure range 

Fixed Point s. 

Triple roint Temperature 
Triple Point Pressure 
Triple Point Density-Vapor 

Liquid 

Critical Point Temperature 
Critical Point Pressure 
Critical Point Density 


54.359 to 400 K 
-361.8238 to 26u.33°F 
0** to 1200 bar 
0 to 17404. psla 

54.359 K* 

-361 .8238°F . 

.001490085 bar** 
.021611856 psla 
.0003275488 moles/ liter 
6.5431796 x 10" 4 lb/cu ft 
40.81997 moles/ 11 ter 
81.542780 lb/cu ft 
154.481 K* 

-181 ,4242°F 
50.42218 bar 
731.31190 psla 
13.63 moles/liter* 

27.228 lb/cu ft 


*This particular value is used to define the point In the sense that it was 
used as Input Into the program. 

**See Section 3 for more detailed explanation of the lower limits of 
pressure. 
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Table 16. Properties for Computer Program Performance Test for Oxygen 



10. Properties of Argon 

The thermodynamic properties of argon have been calculated from a 32 term 
empirical equation of state from Hanley* et al. (1974). The functional form of 
the argon equation of state Is the same as was used for hydrogen* nitrogen* 
oxygen and methane. The transport properties for argon were also taken from 
Hanley, et al. (1974). 
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Table 18. Pressure and Temperature Range for Argon 


Temperature range 
Pressure range 


83.80 to 400 K 
-308.83 to 260.33 # F 
0 ** to 1000 bar 
0 to 14504. psla 


Fixed Points 

Triple Point Temperature 
Triple Point Pressure 
Triple Point Density-Vapor 

Liquid 

Critical Point Temperature 
Critical Point Pressure 
Critical Point Density 


83.80 K* 

-308 .83° F 

.6890708 bar** 
9.9941271 psla 
.1015395 moles/1 Iter 
.25322653 lb/cu ft 
35.40018 moles/1 Iter 
88.283522 lb/cu ft 
150.86 K* 

-188.122°F 
48.98050 bar 
710.40210 psla 
13.41 moles/1 Iter* 
33.443 lb/cu ft 


*Th1s particular value Is used to define the point In the sense that it was 
used as Input Into the program. 

**See Section 3 for more detailed explanation of the lower limits of 
pressure. 
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Table 19. Properties for Computer Program Performance Test for Argon 



* 


11. Properties of Methane 

The thermodynamic properties of methane are calculated from a 32 term 
empirical equation of state by McCarty (1974). The functional form of the 
methane equation of state Is the same as was used for hydrogen, nitrogen, oxygen, 
and argon. The transport properties are from Hanley, et al. (1977). 


Table 20. Properties and Associated Estimated Uncertainties for Methane 


Property 

All Uncertainties In Percent (Except Enthalpy) 

Liquid (below T c ) Gas (below T c ) Fluid (above T c ) 

pressure 

5.0 

0.25 

0.25 

density 

0.1 

0.25 

0.25 

temperature 

0.1 

0.25 

0.25 

enthal py 

1.0 

0.5 

0.5 

entropy 

0.5 

0.5 

0.5 

C P 

2.0 

5.0 

2.0 

Cv 

2.0 

5.0 

2.0 

speed of sound 

1.0 

0.3 

0.3 

thermal conductivity 

5.0 

5.0 

4.0 

viscosity 

3.0 

3.0 

2.0 


The above are based on the estimated average difference between calculated 
and true value. Including problem areas as described In section 2. 
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Table 21. Pressure and Temperature Range for Methane 


Temperature range 
Pressure range 

Fixed Points 

Triple Point Temperature 
Triple Point Pressure 
Triple Point Density* Vapor 

Liquid 


Critical Point Temperature 
Critical Point Pressure 
Critical Point Density 


90.68 to 500 K 
-296.446 to 440.33 # F 
0** to 350 bar 
0 to 5076. psla 

90.68 K* 

-296. 446 °F 

.1174350 bar** 
1.7032507 psla 
.01569094 moles/llter 
.015713995 Ib/cu ft 
28.15114 moles/llter 
28.192503 Ib/cu ft 
190.555 K* 

— 116.671°F 
45.98838 bar 
667.00507 psla 
10.23 moles/llter* 
10.245 Ib/cu ft 


*Th1s particular value Is used to define the point In the sense that It was 
used as Input Into the program. 

See Section 3 for more detailed explanation of the lower limits of 
pressure. 
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Table 22. Properties for Computer Program Performance Test for Methane 
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Appendix 1 


Example of the Use of the Program 


UiLl » IPKuPb 

HI t UM LHMI III I HI WUKLli HI 
PRlIMilO I I I rtlit I Mil h i. <» 


I I HIH iKilPlKill *• II Hill tiRL NM I Ktill I ' '•* wl *•' 


WHIM I Ht PROGRAM ASKS HIH A II U1U GHEl.tlUN* EMI EH iHI ill PKuPP I h 1 1 NHMHFk 
AM LNAPPKUPKI.'i 1 1; NUMBER WILL IE KM I MAIL I HE PROGRAM 
UIILM I HE PKUULAM AblsS I UK A PKfcbbUKL » Dt MU 1 1 V » AMU 
f LMF'LKA I IIK’fc * LMIIK ANY 2 UE IHL IIIKEt ANU A O MlK I HE IHIPli. 

I Hfc UK Ut k MUbl ill P» D» I » ANU UNL UE IHL I HKH MUM HE. U. 

IF ALL (HKLL AKL 0 IHL KKUUKAM ASK Si I UK A NLW FLUID 
IE YUU AKL LNILKtblkD IN A DkliMlIlUN UE IHL UAKXUUb 
MULUF ICAI LUNS OF HYUKOULN LMIER A l» IE NUI LNItK A 0 
Y 0 

IE TUU AKL LNILKtSTED IN IHL bOUKCLS OF IHkSE F'KUUKAfVi 
LNIE.K A 1 IE NUI LNItK A 0. 

<- 0 

btLLCI A ELUIli FROM IHL FOLLOWING L IS I 
FAKA HYDROGEN- t 
NORMAL HI DRUGEN "10 
UK THU HYDROGEN- 11 
HELIUM-I! 

NE0N-.4 
NIIR00LN«4 
OXYGEN- 5 
ARUUN*6 
ML THANE"/ 

r£KMINATIUN=8 


Y 4 

THE RANGE OF TEMPERATURE FOR NITROGEN IS 6.5.15 TO 1900K 
WITH PRESSURES TO 10000 UAR 

HO YOU PREFER ENGINEERING UNITS OR METRIC UNITS Y 
ENTER A 0 FOR ENGINEERING UNITS OR A 1 FOR Hfc I RIG 
? 1 

DO YOU WANT SATURATION PROPERTIES ? 

ENTER 0 FOR YES UR 1 FOR NO 
? 1 

DO YOU WANT A SINGLE POINT OR A TABLE OF PROPERTIES "* 
ENIER A 0 FOR A SINGLE POINT OR A 1 FOR A TABLE 

Y 1 

ENIER A PRESSURE * A STARTING TEMPERATURE* A FINAL 
TEMPERATURE AND A TEMPERATURE INCREMENT. IN BARS. DEGREES 
K. AND IN THAT ORDER 

Y 10.100.200.10 


PRESSURE 

DENSITY 

TEHP 

H 

S 

CP 

cv 

S SOUND 

CONli 

vise 

BAR 

MOL/L 

K 

J/HOL 

J/MOL-K 

J/MOL-K 

M/S 

MW/M-K 

hi \ 

10*00 

24.6342 

100*00 

- 2032.40 

94.64 

64.36 

26*23 

616 * 75 

98*53 

731*3 

10*00 

1.3169 

110*00 

2731.63 

140.55 

40.17 

23*18 

194*20 

12* / / 

79* / 

10*00 

1.1452 

120*00 

3111.38 

143 . as 

36.25 

22*33 

208*92 

13*35 

U5 * 5 

10*00 

1.0226 

130*00 

3462.67 

146.67 

34.19 

21*85 

221*53 

14*07 

91*5 

10*00 

.9282 

140*00 

3797.88 

149.15 

32.94 

21*56 

232*83 

14*73 

97*5 

10*00 

.8523 

150*00 

4122.94 

151.39 

32.12 

21*36 

243*20 

15 * 4 / 

103*4 

10*00 

. 7894 

160*00 

4441.10 

153.45 

31.55 

21*23 

252*87 

16*23 

109*2 

10*00 

. 7360 

170*00 

4754 . 3 S 

155.35 

31.13 

21*14 

261*98 

l /*00 

114*9 

10*00 

.6901 

180*00 

5063.95 

157.12 

30.81 

21*07 

270*64 

17*76 

120*6 

10*00 

.6499 

1 V 0*00 

5370.75 

158.78 

30.56 

21*02 

278*93 

18*53 

126*1 

10*00 

.6145 

200*00 

S 67 S .36 

160.34 

33.37 

20*98 

286*88 

19*29 

tj>i*5 

ENTER A 

PRESSURE. 

A STARTING TEMPERATURE. 

A FINAL 






TENPERATURE AND A TEMPERATURE INCREMENT. IN BARS. DEGREES 
K. AND IN THAT ORDER 
? 0. 0.0.0 


45 





HkIM.r A KUUi FKUM IHt FOU-CIUJMG US? 

F'AKA HVHKUUtN-1 
NOKMAI- HYUKllUf M- 10 
UK FHO HH'KCIGMMII 
Hfel. IUM~iJ 
NfcOM«< 

Ml 1 KI1UI-M 4 
DXYUEMi ti 
AKGQN-i 
t l. JHhMF - / 

If KM IMA I 10N*H 

• II 
I XI I . 
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Appendix 2 


Program Listings 


PROGRAM FLUIDS! INPUTS so » output) 

DIMENSION 0(32) * UP < 9 ) 

COMMUN/DA TA/O « K' * GAMMA * VP * 0 TP * PLC * P fP * ICC* T TP * TUI- * TU. »P UL * BCG 
COMMON /CON TV IF 

CQMMON/CR I T/EM * b OK *RM* TC* PC , X , PC , 8 I u 

TOO FORMA r < * HI* WELCOME TO THE WOK IP Of FLUID PROPER T IPO * IP YOU ARE NOT 
INOT FAMILIAR WITH IMP PPOORAM PI EAST ENTER A 0 HO 
PRINT TOO 
READ 101*1 
ip (I .EO.OCALL INTO 

101 FORMAT < 1 1 > 

99 PRINT 102 

PRINT 104 

102 FORMAT (* SELECT A FLUID FROM THE FOLLOWING LIST*) 

104 FORMAT < * PARA HYDROGEN ^1*/* NORMAL HYDROGEN- 10*/* ORIHO H 

lYDROGEN^ll*/* I IEI IUM"'. 1 '-/# NEUN«3*/* NITROGEN-4*/* OXYG 

2EN«S*/* ARG0N'-6*/* .TET HANE«/*/* 

4 TERM. t NAT I ON'-. 0*) 

READ *» IF 

GO TO < 1 * 2 * 3 * 4 * a * 6 * 7 * 999 * 999 * 10 * 1 1 ) * I F 

1 CALL DATA PH 2 
GO TO 13 

2 CALL DATA HE 
GO TO 99 

J CALL DATA NE 
GO TO 99 

4 CALL DATA N2 
GO TO 13 

5 CALL DATA 02 
GO TO 13 

6 CALL DATA AR 
BO TO 13 

7 CALL DATA CH4 
GO TO 13 

10 CALL NH2 
GO TO .13 

11 CALL 0H2 
GO TO 13 

13 PRINT 105 

105 FORMA! <* DO YOU PREFER ENGINEERING UNITS OR METRIC UNITS ?* 

1/* ENTER A 0 FOR ENGINEERING UNITS OR A 1 FOR METRIC*) 

READ 101 *IU 
PRINT 106 

106 FORMA T ( * DO YOU WANT SATURATION PROPERTIES ?* 

1/* ENTER 0 FOR YES OR 1 FOR NO*) 

READ 101 *IC 
PRINT 10/ 

107 FORMAT (* DO YOU WANT A SINGLE POINT OR A TABLE OF PROPERTIES 
1/* ENTER A 0 FOR A SINGLE POINT OR A 1 FOR A TABLE*) 

READ 101*10 

21 IF'< IC.EG»0)U0 TO 30 
IF < IV « EG* 1 >GQ TO 40 

22 IF < IU*EO. 0 >B0 TO 19 
PRINT 103 

READ * * P * D * T 
P=P/ 1.01323 
GO TO 20 
19 PRINT 100 

108 FORMAT < * ENTER PRESSURE IN LB/BO IN* DENSITY IN LB/CU FT* AND* 

1/* TEMPERATURE IN DEGREES F*> 

READ * * P * D * T 
P*P/ 14. 695949 
D=0*16. 01846371/EM 
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IF! I.EQ. 0.0)60 TO 20 
T»!T-32. >/l.fl+273. IS 
20 IF!F’,LE. 0.0)60 TO 14 

103 FORMAT!* ENTER PRESSURE IN BAR» DENSITY IN MOLf.S/LITfcR» AND TEMPER 
1ATURE IN KELVINS*) 

If !D.LE. 0.0)60 10 17 
IF <T .LE. 0.0)60 TO 12 
60 TO 22 

12 IF<P. LE.O. OR. D.LE. 0,0)60 10 99 
T»FIND T ! P » 0 ) 

CALL LIMITS!P»T » IL) 

IF! IL«LE«0)60 TO 22 
60 TO 16 

17 IF<T,LE.O.OR.P,LE.O)60 TO 99 
CALL LIMITS!P»T»IL> 

IF< IL.LE.0)60 TO 22 
D=FIND D!P,T> 

GO TO 16 

14 IF(D.LE.O.O.OR,T.LE.O)60 TC 99 
P=FIND P ! D » T ) 

CALL LIMITS! F‘» T • IL> 

IF! IL .LE .0 )60 TO 22 

16 CALL REPRO!P>D»T>IU>lV»ICr IP»TF»DELT> 

60 TO 22 

30 PRINT 109 

109 FORMAT!* DO YOU WANT SATURATED LIQUID OR SATURATED VAPOR* 

1/* ENTER A 0 FOR LIQUID OR 1 FOR VAPOR*) 

READ 101 t IP 

IF ! IV.EQ. 1 )G0 TO 40 

PRINT 111 

111 FORMAT!* DO YOU WANT TO ENTER WITH TEMPERATURE OR PRESSURE * 

1/* ENTER 0 FOR TEMPERATURE OR 1 FOR PRESSURE*) 

READ 101 , II 
IF! II.EQ, 1 )GQ TO 35 
211 IF! IU.EQ. 1 >60 TO 31 

112 FORMAT!* ENTER A TEMPERATURE IN DEGREES F*> 

PRINT 112 

READ *t TI 

T=!TI-32,)/l, 8+273. 15 
IF'! T.L.E ♦ 0 , 0 )G0 TO 99 
GO TO 33 

31 PRINT 113 

113 FORMA!!* ENTER A TEMPERATURE IN DEGREES K*> 

READ *»T 

33 IF!T.LE. 0.0)60 TO 99 
IF!T.GT.TCC.OR,T.LT.TTP)GO TO 34 
P»VPN!T) 

IF ! IP ♦EQ.O)P=‘P+. 00001 
D=FIND D ! P » T ) 

CALL RE PRO!P»D»T»IU»IV»IC» TPfTF fDELT) 

GO TO 211 

34 PRINT 114»TCC» TTP 

114 FORMAT!* YOUR INPUT TEMPERATURE IS OUTSIDE THE RANGE OF THE* 

1/* SATURATION TEMPERATURES FOR THIS FLUID.* 

2/* TC«*,F6.2>* TTP=*,F6.2»* TRY AGAIN*) 

GO TO 211 

35 IF! IU.EQ. t )G0 TO 36 
PRINT 115 

115 FORMAT!* ENTER A PRESSURE IN LB/SQ IN*) 

READ *»PI 

IF! PI ,LE.O.O)bO TO 99 
P»PI/14. 695949 
GO TO 37 
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3 6 PRINT 116 

116 FORMAT!* fc.NI KR A PRESSURE IN WAR#) 

READ *,Pl 

ifipi.le.o.ohm) ro 99 

P*PI/1 01 325 

67 ir<p.ur.PCc.0R.p,i-r.Ptp»uu ro 38 
! =F I NU fV!P> 

P~VPN! 1 > 

tr< ip.lg.o>p*pf,o<>oi 
d^find n<p»r- 

call re pro<p»ii» r* iu,iv* rc»ip. i> »delt) 

GO TO 35 

38 PRINT U7,PCC,PTP 

117 FORMAT!* YOUR INPUT PRESSURE 18 OUTSIDE THE RANGE OF SATURATION# 
1/* PRESSURES FOR THIS FLUID*# 

2/# PC=#,F6.3>* P I F‘~# » F 6 . 5 » # TRY AGAIN#) 

GO TO 35 

40 IF < IC.EG. 1 )G0 TO 50 
IF ! 1U.E0* l 'GO TO 41 
PRINT 118 

118 FQRMAK# ENTER A STARTING TEMPERATURE, A FINAL TEMPERATURE# 

1/* AND A TEMPERATURE INCREMENT, IN DEGREES F AND IN THAT ORDER#) 

READ *»TS,TF ,OELT 
IF!DEL.T.LE. 0.0)60 TO 99 
TS=! 1 3-32. )/l. 8+273. 15 
TF= ( TF- 32 . ) / 1 .8+273 . 15 
DELT=DEL 1/1.8 

1F!TS.LF.TTP.0R.TS.GI ♦ TCOGO TO 45 
IFITF.LT. TTP.OR. TF .GT.TCOGO TO 45 
GO TO 42 

41 PRINT 119 

119 FORMAT < # ENTER A STARTING TEMPERATURE, A FINAL TEMPERATURE# 

1/* AND A TEMPERATURE INCREMENT IN KELVINS AND IN THAT ORDER*) 

READ #,TS,TF ,DELT 
IF(DELT.LE.O.O)GO TO 99 
IF!TS.LT. TTP.OR.TS. GT.TCOGO TO 45 
IF!TF.LT« TTP . OR . TF «GT ♦ TCC )G0 TO 45 

42 T=TS 
F»VPN!T> 

IF!IP.EG.0.0)P=P+.0001 
D=FIND D(P,T) 

CALL RE FRO!P,D,T,IU,IV,IC,IP,TF,DELT> 

GO TO 40 

45 PRINT 120,TCC,TTI 

120 FORMAT!# EITHER YOUR STARTING OR FINAL TEMPERATURE IS OUTSIDE# 

L/# THE RANGE OF SATURATION TEMPERATURES.# 

2/* TC=#»F6.2,* rTP®*,F6.2,* TRY AGAIN*) 

GO TO 40 

50 IF<IU,EQ,1)G0 TO 51 
PRINT 121 

121 FORMAT!# ENTER PRESSURE, STARTING TEMPERATURE, FINAL TEMPERATURE* 
1/* AND A TEMPERATURE INCREMENT, IN LB/CU FT, DEGREES F, AND IN* 
2/# THE ABOVE ORDER*) 

READ #,PI,TS»TF ,OELT 
IF!DELT.LE.O.O)GO TO 99 
P»PI/1 4 .695949 
T«!TS-32.)/l, 8+273. IS 
TF=!TF-32. )/l. 8+273.15 
DELT=DELT/1 .8 
CALL LIMITS!P,T ,IL) 

IF! IL.LE.OGO TO 50 
CALL LIMITS!P»TF,IL> 

IF! IL.LE.O)GO TO 50 


50 


go ru s<? 

51 PHI INI 1.2;-' 

122 FORMAT <* FNTFk A PRESSURF.:* A START ING TEMPERA ) URL * A PINAL* 
1/* TEMPERATURE AND A ( 1 MPF.RA TURl. INCREMENT * IN BARG* DEGREES* 
2/* K * AND IN FHAT ORDER*/ 

READ **P.t»TB* FF*0F.LT 
IF (DELI .ll,,0.0>0U 10 VV 
r®rs 

P' a F'l/ 1 *01325 
CALL L1M1T8CP* I • XL > 

IF'< 1L.LE .0)00 TO 50 
CALL L1M1 TS<P* TF"» IL > 

IF < IL*LC*0)GO TO 50 

52 D®F I Nil D(P*r> 

CALL RE PR0<P*D*T*1U*1V*1C* JP*TF»DELT) 

00 TO 50 
999 continuf: 

STOP 

END 

SUBROUTINE DATA ET 
PRINT 100 

100 FORMAT <* ETHANE IS SCHEDULED TO BE ADDED TO THE FLUIDS* 

1/* PACK IN THE NEXT FISCAL YEAR*) 

RETURN 

ENTRY DATA PR 
PRINT 101 

101 FORMAT (* PROPANE IS SCHEDULED TO BE ADDED TO THE FLUIDS* 

1/* PACK AFTER FY78*> 

RETURN 

ENTRY DATA NH3 
PRINT 102 

102 FORMAT (* AMMONIA IS SCHEDULED TO BE ADDED TO THE FLUIDS* 

1/* PACK DURING FY78#> 

RETURN 

ENTRY DATA SH 
ENTRY DATA SO 
PRINT 103 

103 FORMAT (* THE SLUSH PROGRAMS WILL BE ADDED TO THE FLUIDS* 

1/* PACK IN FY78* ) 

RETURN 

END 

SUBROUTINE RE PRO( P*D*T» IU» IV* IQ* IP* TF * DELT) 

DIMENSION G<32)*VP<9> 

COMMON/DATA/G » R * GAMMA * VP * DTP * PCC * F TP * TCC * T TP * TUL * TLL » PUL 
COMMON/CONT/IC 

COMMON/CRIT/EM * EOK *RM » TC * DC * X * PC * SIG 
N=500 

IF< IV*EQ*0)TF=T-1 ♦ 

PRINT 100 
IFUU.EQ,0>G0 TO 1 
PRINT 102 
GO TO 2 

1 PRINT 103 
PRINT 104 

2 CONTINUE 

DO 10 1*1 *N 

IF( I «EQ. 1 )G0 TO 3 

D=FIND D < P » I) 

3 H“ENTHAL(F*V 

8®ENIR0P(D*T) 

W°80UND<D*T ) 

CPP=CP<D*T) 

CVV«CV(D*T> 
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Q-Ptl. 01325 

IF< IU»EQ,0)60 TO 21 

IF<IC.EQ.1.0R.IC«Ea.lO.OR,IC,EQ,ll>00 fO 20 

IF<IC.EQ.9)G0 TO 20 

THoTHERM(D'T) 

V-VISC(D.T) 

PRINT 101»Q»D»T»H»S»CPP»CVV»U TH»V 
00 TO 9 


20 PRINT lOltQ'DiT'H'SfCPPtCVV'U 
00 TO 9 

101 F0RMAT<F8«2*F8.4»F9.2f F9.2>F8.2»2F7.2*F8.2»F 7.2»F 7. 1 ) 


TEMP 


H 


S 


DEG F BTU/LB 


100 FORMAT < $ PRE88URE DENSITY 
1 SOUND COND OISCO 

102 FORMAT < * BAR MOL/L 
1M/8 MW/M-K Ml0/CM->*) 

103 FORMAT <* LB/8Q IN LB/CU FT 
1 F/S BTU LB/FT-8*) 

104 FORMAT (* 

1 FT-HR-F E+7*) 

21 M«H/ < 2 ♦ 324445BEM ) 

8«8/(4» 184001BEM) 

CPP«CPP/ < 4 . 184001 «EM ) 

CW-CW/ < 4 ♦ 18400 ITEM ) 

W«W*3. 280840 
P0«Ptl4. 695949 
D0«D*EM/16. 01846371 
T0-TB1. 8-459 *67 

IF< IC .EQ. 1 .OR, IC.GT .8)G0 TO 22 
TH-THERM (D*T )*. 000o78176 
V®VI8C(D»T) *.0067196897 
PRINT 105»P0»D0»T0>H»S>CPF’fC0VfU>THfV 

105 FORMAT < F8 ♦ 1 » F8 » 3 » 2F9 » 2 » F8 ♦ 2 » 2F7 ♦ 3 » F8 . 0 » F7 . 4 » F7 ♦ 3 ) 
60 TO 9 


CP 


CV 


J/MOL J/MOL-K j/MOL-K 


BTU/LB-F 


22 PRINT 101rP0»D0»T0»H»S»CPP»CW»g 
9 T»T+DELT 


IF<T,6T.TF+, 01)60 TO 11 
IFCIQ.EQ.O)P=VPN(T) 
IF<IP.EQ.O) P=P+ ♦ 0001 

10 CONTINUE 

11 CONTINUE 
PETURN 
END 


8 
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SUBROUTINE DATA CH4 
DIMENSION G(32)?VP<9)rUI<ll> 

DIMENSION GV<9>fGT<9).FV<4>?FT<4>?EV<B>?ET<8> 

DIMENSION A (20) 

C0MM0N/SEN/BETA?XQ?DELTA?E1 » 12 1 AGAM 
COMMON/CP I D/0 I 

COMMON/CRIT/ EM? EOK» RM? TC» DC? X ? PC? SIG 

CQMMQN/DATA/G?R?GAHMA?VP?DTP? PCC?PTP? ICC? T TP? I UL?TIL?PUI ?DCC 

CQNM0N/DATA1/GV?GT?FV?FT?EV?ET 

COMMON/ ISP/N ? NU ? NWU 

COMMON/ SATC/A?DTPV 

NWU-0 

PRINT 100 

100 FORMAT <* THE RANGE Of TEMPERATURE FOR METHANE IS 90.60 TO SOO K* 
1/* WITH PRESSURES TO 350 BAR*) 

N=0 * NW=0 

EM=16.042 * E0K=168. ♦ RM=4.101E-08 ♦ TC*1V0.555 

DC*0» 162/ * X*1.7l24 ♦ PC*45»95 * SIG*3.68 

X0-0.164 ♦ BETA=0.355 * DELTA=4.352 * El*2.03 ♦ E2*0.287 

AGAM= 1.190 

PARAME1ERS FOR TRANSPORT 
PROPERTIES - HANLEY ET. AL. JOURNAL 
OF PHY. CHEM. REF. DATA VOL 6 NO 2 
19/7 

DATA SUBROUTINE FOR METHANE 
GV ( 1 ) *-2 . 0909 74 /942E+2 
GV( 2 ) =2.6472692181E+2 
GV<3)=-1 .4728175613E+2 
0V(4)»4./16/401921E-F1 
GV<5> *-9. 4918721 789E+0 
GV<6>*1 .2199792872E+0 
GV(7)«=-9.62799355746E-2 
GV(8)»4.2741516570E-3 
GV( 9 ) *- 8 . 141530 7247E-05 
GT ( 1 >=-2. 1476213125E+5 
GT ( 2 ) *2 . 1904610575E F5 
GT(3)*-8.61809/3719E+4 
GT < 4 ) = 1 . 4960986936E+4 
GT < 5 ) *-4 ♦ 73066031 77E+2 
GT(6)*-2.33U779643E+2 
GT < 7 ) *3 . 7784390759E+1 
BH8)«-2, 3204806092E+0 
GT ( 9 ) **5 . 31 1 7637687E-2 
EV < 1 ) *-l . 0350606586E-FI 
EV < 2 ) =4 . 2903609488E-2 
EV<3>*1.7571599671E+i 
EV(4>“6. 1276818706E+3 
EV<5>— 3.019 5918656E+3 
EV(6)®1 .887301 1594E-F 2 
EV(7)*1.4529023444E+2 
EV<8)“0» 162 
FV< 1 >=1 .6969859271 
FV<2>*~. 13337234608 
FV<3>=1 .4 
FV<4>«168. 

FV(4)-*168. 

FT<1 >*-.252762921 
FT <2)“* 3343285931 
FT(3)*1. 12 
F I (4)=168. 

ET ( 1 ) *-7.0403639907 
ET< 2)-”. 74421462902 
ET(3)*12. 319512908 


S3 


ETT4)»2.2209758501EE3 
ET < 5 ) “-8 , 8S2597V9 33E f 2 
El T6>«72. 83589 7919 
ET< 7 >«~2. 970691 4540 
ETT8)=0.1628 

C PARAMETERS I- UR THERMODYNAMIC 

C PROPERTIES - MCCARTY' CRYOGEN.lt S 

C VOL 14 NO 5 MAY 1974 AND GOODWIN 

C MBS TECH NOTE 653 APRIL 74 
R“. 082056 16 
GAMMA*- ♦ 0094 
GIT 1 )*-l «8044/50507E+A 
A( 1) = .1 83603246 136E +02 

A< 2) = - . 182553840603E- 0 l 

A< 3>* .5866238071/8EE00 

A< 4)* - < 106005894683E+02 

A( 5)* » 16/157622432EE03 

A( 6> = -. 15556917321 7E+03 

A< /)* . 106609415022E+03 

AT 8)* -.341087933790E+02 

AT 9> = .973203073452E+01 

AT 10)* 388410018388E+02 

AT11)* . 5295734547 7 IE f 07 

AT 12)* - . 291075304738E +02 

AT 13)* » 58530764 7478L 1-0 1 

AT 14)* 6372765321 86E TOO 

AT 15)* ♦ 994444 1 09622E-03 

A T 1 6 ) * - . 39231 482 L657E -01 

AT 17)* ./ 28820880 748EE0v 

AT 18)* - .62582 181 531 5E f 01 

AT191* .295561390641L TO 1 

AT 20) = - . 763972A49504F t-0'' 

GIT2)*7. 7426666393E+4 
GI T 3) *-l *3241658 754L' 1 T 
61 T4)*l. 54381 49595F+1 
B1 T5)*-5. 14 79->05757E-2 
GI(6)*1 .0E091721V6E--4 
61 T 7 > —6 . 550 l /4343/L-8 
81 1 8 * * 6 . 74900561 71E+0 
’ lT9.**.3Ef4 
GIT 10)*G1 Til )*0 
10*298.15 
H0«19196. 1 
GIT10)=H0-HI ( TO) 

T0*300. 

80*186.4854941 
Gl \ 11 )*SQ-SI T TO) 

GT 1)* -.187027997685E-01 
GT 2)* .103387108009E+01 

GT 3)* -♦ 155387625619E+02 
GT4)* . 77231 1478564E+05 

GT 5)* - ♦ 377103300895E+ 05 
GT 6)* . 84681 8843475E-03 

GT /> = - . 496415884S29E +00 
GT 8)* . 86990935241 4E+02 

GT 9)* -.322821S92493E+05 
GT 10>* -. 3958430263 18E -04 
GT11)= .266 7723 18035E-01 

GT 12)= -♦ 30401 005 7839EE01 
GT 13)* . 191584507536E-0 1 

GT 14)* * .195587933458E-03 
GT 15)* . 60 74 799678 79E+01 

GT 16)* - . 529609525VU4E -03 
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0(17)“ ,1522642860046-04 

G(18>* -. 109952182842E-01 
G<19>* .19139S549929E-03 

G<20>* .386470003744E+05 

G<21 )* - . 157930582612E+07 
G<22>* . 1952701 44401E+03 

G<23>* .1659940814296+07 

0<24>* .6030511467116+00 

G<25)« .37648S162808E+02 

G<26>* .125593480622E-02 

G<27>» -.3435700325136+02 
G<28>» -.5409450941396-05 
G<29>» . 185622284663E-02 

0 < 30 ) » .7707869792456-08 

G<31 >» - . 2868683 18650E-05 
G<32>* .3723769616476-04 

DTP®. 28151143814236+02 
DTPV* . 15900415451606-01 
VP<1 >*4.77748580 
VP<2>*1. 76065363 
VP< 3 >*-.56 788894 
VP(4)*1 .32786231 
VP < 5 ) * 1 . 5 
VP( 6) *.1158993 
VP(7)*90.68 
VP(8 )*190.555 
VP<9>*0«0 
TCC*VP<8) 

PCC»VPN< rco# 1.0 1325 
PTP»VP<6)*1. 01325 
TTP»VP(7> 

TUL*500. 

TLL*TTP 

PUL*350. 

DCC* 10. 23 

RETURN 

END 

SUBROUTINE DATA N2 
DIMENSION G ( 32 > . VP <9).GI<11) 

DIMENSION GV(9> »GT<9) »FV<4> » + T<4> rEV(8> »ET(8> 

DIMENSION A(20> 

CONMON/CRIT/ EM. EOK. RM» TC» DC. X » PC. SIG 
COMMON/SATC/A » DTPV 
COMMON/ DATA1/GV » GT » FV » FT » EV » ET 
COMMON/SEN/BETA. XO. DELTA. El . E2. AGAM 

COMMON/DAT A/G » R » GAMMA . VP . DTP . PCC .PTP.TCC. TTP.TUL. TLL . PUL » DCC 

COMMON/CP I D/G I 

COMMON/ 1 SP/N » NU » NUU 

NWU=0 

PRINT 100 

100 FORMAT < * THE RANGE OF TEMPERATURE FOR NITROGEN IS 63.15 10 1900K* 
1/* WITH PRESSURES TO 10000 BAR*) 

N*0 * NU-1 

EM=28.016 * E0K*118. * RM*3.933E-08 • TC*126.24 

DC*. 3139 * X*l. 67108 $ PC*33.98 ♦ SIG*3.54 

X0*0. 164 * BETA*0.355 ♦ DELTA*4.352 * 61*2.17 * E2*0.287 

AGAM* 1.190 

C PARAMETERS FOR TRAN8P0RT 

C PROPERTIES FROM HANLEY ET. AL. 

C JOUR. PHY. CHEM. REF DATA V0L3. NO 4 

C 1974 

GV(1 >*-1.82242400006+2 * GV<2)*1 .99153273/4E+2 

0V<3>*-9. 1542324494E+1 * GV(4>*2.3255484059E+1 
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GV < 6 > =3 . 64575068 11 E- 1 

♦ GV<8)«7.00539O4895E~4 

* GT< 2 ) =4 • 9765746684E+3 
GT<4)“-5.5022716888E+3 
GT < 6 ) o-2 . 2974737237E+2 
GT <8)“-8*5677385768E-l 


♦ 

♦ 

* 


GUO) «-3 *6307214228 * 

GV < 7 ) »-2 . 22618808 1 7E-2 
GU<9)«-1 , 1894029104E-S 
GT < 1 >*-2.0029573972E+4 
GT<3)=8«018B959378E+3 
GT < 5 > = 1 * S363738965E+3 
GT < 7) »1 . 9360547346E+1 
GT<9)»1 .5564670935E-2 

EV< 1 ) o-9. 8255690362 • EV(2)“S« 7156092139E-001 

EU<«! , ’ 1 a 6 ftB2oiil}? 8E+001 * EU<4) t1 3*695408615BE+003 
^ * “-e*088980U80E+002 « EU<6)«6.8464435640E+001 
cn / ? » ®“2* 1241135912 ♦ EV<8>»0.313 

hUl)— 1.1217739623 % FV<2)=0. 32912317244 
FU<3)“1 .4 $ FV<4)=118* 

* FT <2)= *06102791 1104 
FT ( 3) =1 *2 * Ff (4)=118* 

! ET(2>*3. 72017, 3333E+1 
3*9013509079E*f 1 $ ET<4 )« m 3« 1S2610948SF+1 

THERMODYNAMIC PROPERTIES FROM 1S-6109485E+1 

NBS TECH NOTE 648 - DEC 77 
JACOBSEN ET* AL. 

R-8.20539E-2 


GAMMA»- * 0056 

GI<1)=> -0. 735210401 1S7252E 03 

GI(2)a 0. 3422399804 11978E 02 

GI<3)e -0.557648284567620E 00 
GI (4)= 0 , 350404228308756E 01 

GI<5)= -0. 173390185081005E-04 

GI(6)* 0.1 74650849766463F -0/ 

GI<7)» -0 . 3568920335443 48E- 1 1 

GI<8>« 0. 100538722808834E 01 

GI <9) = 0. 3353406 10000000E 04 

T0=298. 15 * S0=191 .502 * H0=8669.0 

GI(10)=GI(11)=0 
GI<1C)=H0-HI<TQ> 

GI < 1 1 ) "SO— *51 < f 0 ) 

G< 1)= 0.136224769272827E-02 

G< 2 ) = 0. 10703246990859tE 00 

G< At- -0. 24390072 1871 41 3E 01 

G< 4> = 0, 341007449376470E 02 

G( 5> = -0 « 4223743094661 67E 04 

G< 6> = 0. 105098600246494E-03 

G< 7) = -0 « 1 12594826522081 E-01 

G< 8> = 0.1 42600789270907E-03 

G< 9)» 0.1 84698501 60900 7E 05 

G< 10)“ 0*8111 40082588 776E-07 

G(ll)“ 0. 23301 1645038006E-02 

G< 12)= -0.507752586350986E 00 

G< 13)= 0.485027881931214E-04 

G<14>= -0,1 136567641 15364E-02 

G< 15)= -0 « 707430273540575E 00 

G< 16)= 0,751 706648852680E-04 

G< 17) = -0*1 116141 19537424E-05 

G < 18 ) = 0 . 368796562233495E-03 

G< 19)= -0.201 31 769 1347729E-05 

G ( 20 ) = -0*14971 7444755949E 05 

G<21)« -0* 1 19719240044192E 06 

G(22)= -0.975218272038281E 02 

G<23)= 0. 55463971 3151 823E 05 

G<24)= -0*1 79920450443470E 00 

0 •, 25 ) » -0. 2565829260771 84E 01 

G<26)= -0. 413707715090 709E-O3 

G<27)= -0. 2562454 15300293E 00 
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Q<28>* 
U<29>* 
G<30>« 
0(31 >■ 
A# 1 >a 


A< 8>« 
A< 9)» 


A<10) = 

A<11)» _ 

A(12)a ~. 21 6767601 780E+02 

A<13)» .496S5822647tE+01 

A<14)“ .21 830792847 7E+02 

A(15)= -.126493309807E+00 

A(16)= « 241544188633E+01 

A(17)a -.245256871 794E+02 

A<18)« .9359252071 24E+02 

A<19)« - . 360938251632E+02 

A<20)*= . 757453271989E+01 

0(32) » 0,5853671 7206952 IE-0 7 

VP < 1 >=5, 1113192094 ♦ VP(2)=6.48266/539E-1 

VP(3)a-l . 5108730916E-1 $ VP<4)*7*402849334?E-1 

. w(7,=«;i 

DTP** .309771 7741 4 77E+2 
DTPV« . 242822085710E-1 
TCI>VPt8> 

PCC*VPM rCC)*l « 01325 
PTP»VP<6>*1, 01325 
TTP«VP(7) 


-0, 124222373740063E-06 
0. 103556535840165E-04 
-0.538699146558303E-09 
-0, 7574 l 54 1 28395968-08 
-.158453465S07E+02 
. 41913691 1423E- 01 
-.101965371660E+01 
, 134763743799E+02 
-♦ 109930399087E+OJ 
.925518835497EF02 
-. 95623383 1320EH>2 
.100104366710E+03 
-. 701857937398E+02 
♦900076998647E+01 


286981 120347E+02 


TLl -FTP 
TUL®1900, 
F’UL =>10000. 
DCC=>11.21 
RETURN ♦ 


SUBROUTINE DATA AR 
DIMENSION G<32>,VP<9)»GI<11> 

DIMENSION GV<9),GT<9).FV<4)»FT(4)»EV<8).ET<8) 

DIMENSION A (20) 

C0MMQN/SEN/BET..,XQ*DELTA»E1 * E2» AGAM 
COMMON/SATC/A » DTPV 

COMMON/CRIT/ EM. EOK. RM» TC» DC. X * PC* SIG 

COMMON/DATA/G * R * GAMMA » VP » DTP . PCC * PTP . TCC » T TP » TUL . TLL . PUL » DCC 

C0MM0N/DATA1/GV.GT .FV.FT.EV.ET <lu.»-ul,DCC 

COMMON/CP I D/G I 

COMMON/ 1 SP/N , NU * NUVI 

NWWO 

N>»0 ♦ NW 1 

EM=>39.948 t E0N=152.8 • RM=>3 . 669E -08 ♦ TC®150.725 
DC*0.S33 ♦ X=1 ♦ 7124 $ PCM7.983 $ SIG®3.297 


X0«0. 183 * BETA»0.35S * 

AGAM- 1.1 90 

TRANSPORT PROPERTIES 
FROM HANLEY ET. At . 

J. PHY. CHEM, REF DATA VOL 
NO 4 1974 

GV< 1 )«6. 1 145472787E+1 
GV<2)«-1 ♦ 03943903 12E+2 


DELI A®4, 352 


El®2.27 »E2®0.287 


Usd 
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c 

c 

c 

c 


GV<3)=6. 759461 4619E 4 1 
C» V ( 4 > « -2 . 2536509380E 4 t 
GV ( & ) -4 . 2&939501 38 


OV < 6 ) * - 4 . 72526 7 1 093E - 1 
UV( 7 > -»3 , 1 7V527542SE-2 
0V<9)-i-l .1429083 'HOE -3 
GV(9>= 1.8043010592E-S 
0T(1)*6.2777703742E43 « 

OT (3) 'S.8887549191E4J 4 
0r(5>«3.4884371437E42 4 

OK 7) =-2*5207283167 4 GT(8>«-9. 1098744478E-2 
OT < 9 ) a 1 , 3990842942E *3 
EV( 1 )«-l . 001099399 3E4 I 
EV( 3)=1 *6029145122£41 
EV < 5 > « -S ♦ 6995898780E42 
EV( 7)«3*9870122403E41 
Fallal, 4653452433 4 FV(2)a~0* 77487424965 
FV(J)el,4 $ FV(4)=152*8 

* I! I! - * 2414210327 4 FT < 2 ) « . 075696234255 

F f ( 3) = 1 , 4 FT (4) >152*8 

ET ( 1 >a-2»4116686960EM 


0T< 2)o~9. 6096376637E43 
0 T < 4 > «- 1 . 8920926320E 4 3 
61 (4>*- 3.80167B6193E41 


EV< 2 ) n 2 *0694685 7 12E-1 
EV( 4>®1 . 171 7461 351 E43 
EV<6)a4 ,01 36071 933F41 
EV<8)a,537 


3 » 06948599"’ t £4 1 


ET(2)= 

ET(3>®2. 29S6551674E43 

tT<4) “ -3.5559415848E42 

TRANSPORT PROPERTIES FROM 
HANLEY ET* AL. J* PHY* CHEM* REF 
DATA VOL 3 NO 4 1974 AND 
60SMAN ET. AL* NSRD8 MON 27 1944 
Re. 08205416 
GAMMA®- ♦ 0055 
ei(l)sOK2)BQI(3)oO 
Gi(4)»GI(9) a 2.5 
GI<5)*GI<6>«BI<7)«0 
GI(8)«0 

GI(10)=GK11)*0 
T0®87*28 4 S0“129« 1786 

GI(10)*H0-HI(T0) 

GI ( 11 )»80-SI(T0) 

G < 1 ) *3 . 4342657242351E-3 
G < 2 ) =5 . 7857036681 387E-2 
G < 3 > » - 2 . 69824 708 1 2264E40 
G(4)*>1 ♦ 4481655285291E42 
G<5)«-1 .2849472420416E44 
G < 4 ) »-3 ♦ 2636490894686E-4 
G<7)“2*4629470t90841E-l 


4 H0»9504.8916 


G ( 8 ) »-6 . 9585445697842E4 1 
G(9)«l *9196156939788E44 
0 < 1 0 ) « 1 . 6403909805594E -5 
G< 1 1 ) »- 1 , 0860316345366E-2 
G< 1 2 > a 3 « 323 1 759004885E40 
G < 1 3 > »2 . 1 77636 1 947053E-5 
G< 14)*5* 1615085812771E-3 
G< 15)»-1 . 1 366705407293E40 
G< 16) ®-2« 901851761885VE-4 
G < 1 7 ) »3 . 7898289698060E-6 
G( 18)=1 ♦ 1 03048979098 7E-3 
0(19) =»-l* 4674092942955E-5 
G(20)»-l . 147961 0716 179E44 
G<21 )“-3*9393312963830E45 
0(22) »-9 . 9620084307336E4 1 
G ( 23 ) — I . 857534704601 1E44 
G(24) a -2* 9393483871 136E- 1 
G ( 25 ) » 1 ♦ 6408588086762E41 


j 
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I»l .** > - 4,044/1 / 47 'S’ + 4 

<-»<//; :» , 08/000 / 1 A V/ / 6 l HI 

r*CH) fi.<!/Av<;*oA<>VA»;»f t 
0(29) /«{*| 4 .*405/ / .'A 6 /+ - 5 
0<40> .*,H<V. /4.'Vi182<>/t |o 

Util) .1,400 4//2UAA834f -M 
U< CM ;'.5:*H/41344093AF A 
a< l) •./702429.'3 ///e+07 
A< ?> .) 4104024 1 HAAt +00 

A( .O .2A/4HA448I28I tOl 

A< 4*- #.1001 /AH0440AE+02 

ft*. '1 > .8/58991 494/61 to.* 

A( A' . 40H2A/44A456E+02 

A< /*•• . 104:*A80AA4ML tO.l 

A* ' * 1 « A / I .V»' i' > 4 <'VI. 1 y.‘ 

O' ♦' 1 < I 1 1 <>0. , 9.V#/».| l I <>.• 

O' 10*' . .Ml. "I V.4AA </* 10 * 

AM J > - 4 ,. I <*a:*i •■•♦mi H* 

A< I / > .42/149 /Ojsmvvi 1 11 • 

A* Hi- , 1005999 <km*.m>I (••*/ 

A* 14* . I .</A8/0847O'#l »o.> 

A < lb*- .AA46.<0.lA<|vj| 

A(lA*» . 133.448/117 ‘401 TOl 

All/)- 1 444/1 4432441 H*.' 

A\18>- , 6019484 //OOwM 0 * 

A< IV ) - -#240a8846488/E-l0/ 

A ^ 1*0 > - .46541845888/1 tOl 

VP<1 )-- 4.4151 115519 ( VP ( 7 ) • I .19108(2519 

VP < 4 * -3 . 40/632444+ - 1 % VP( 4>-U. 9555855251+ 1 

V+* ( 5 > -• 1 .5 $ VP (6) .68005 * 01(7 1-84.80 % VP(R) - 1 bo. HA 

VP(9>=0.0 

OTP- ♦ 454002761 yi88t +07 

D1PV-. 102922/02/429 
PRINT 100 

100 EORMAM* ) HE tlMPERAIUKE RANGE FOR ARGON IS 84.8 10 400 K# 
l/# W1IH PRESSURES 10 100O bAKt) 

ICC-0P<8) 

PCL-OPN< 100*1 .01325 
P7P*VP(6>*1 .0142b 
I ?P=VP(/> 

TUI-400. 

tll*ttp 

PUI »1000. 

DCC=l3.4t 
RE I URN I END 
SUBROUTINE DA I A 02 
DIMENSION G( 42)»VP<9)#U1(11> 

‘J I MENS I ON G9<9) »G7 (9) .1 V<4.' #1 T (4) »EV(8> »EI (8) 

DIMENSION A (20) 

COMMON/ SATC/A # D TPV 

COMMON/DAI A 1 /G V # G T » F V » E I #EV#Et 

COMMON/SEN/BE TA# XU# DELTA# El # E2» AOAM 

COMMON/CRIT/ EM# EOK# RM# It# DC# X # PC# SIS 

COMMON/DAT A/G » R # GAMMA » VP # D IP » PCC » P I P > ICC # I I P » TUL » 1 1 1 # PUL # DCC 

COMMON/ I SP/N»NW#NWW 

C0MM0N/CP1D/GI 

NUW=0 

N=0 4* NW=1 

X0=0. 184 * BETA-0. 355 ♦ DEL I A "4. 452 * El=2.21 * E2-0.28? 

EM-31. 9988 * EOK =114.0 ♦ RM*3.8896E-08* TC=154.5/S 
DC-0 . 4362 * X*2. 210636 ♦ PC= 49.7/0 ♦ HtG=4.44/ 

AGAM = 1.190 

C TRANSPORT PROPERTIES 
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FROM HANLEY ET. AL. 

VOL 3* NO 4r 19/4 

GV< 1 )»-V • 7076378593E+1 ♦ GV<2> “B.2B012S4201E+1 

GV<3)-~ 2.4668758803E41 * OV< 4 )«2. 1324360213 


GV < 5 ) »3« 7851049522E-1 
QV<7)*1 . 1134441 304E-2 
GV<9)“1.0279379641E-5 
GT< I )«-2.039S05219JE+S 
GT<3)«- 1.20141 7518 JE+5 
GT(5)“-S • 4244239598E+3 


• GV<6)«"1.Q487216090E-1 

• UV<8>"”5»36760937Q7E"4 

* GT<2)®2.4O081417O?E+S 

* GT<4)“3.295494919EE4 
« GT < 6 > “5 • 4734865540E+2 


GT ( /)“-3«28S4821539E+l • GT<B>™1 .0733572103 

GT<9)“- 1.4610986820E-2 

EV<1)“- 1.215238701 7E+1 • EM<2>“-3. 142172B994E-1 

EV<3)“1 .8201161468E+1 * EV<4)“2.7390429S25E+2 

EM<5>— 2.7498956948E+3 * EV<6)«2.434O609667E+2 

EV(7)“1.1911504104E+2 * EV<8>=.435 

FV<1 >*4.3526515153 • FV<2>« -2.0361263878 
FV<3>“1 .4 * FV<4>«100. 

FT<1>*.3060 ♦FT<2)“«2785 

FT<3>®1.12 4 FT<4)“100« 

ET<l)«-1.2310400765E+t « ET<2>*1 .6799S04261E41 

ET(3>“-2.9944878721E+3 ♦ ET<4)»4.7350508788E+2 


THERMODYNAMIC PROPERTIES 
FROM WEBER - TO BE PUBLISHED 
R“8.20539E-2 
GAMMA"- ♦ 0056 

OKI)" -0 . 49819985371 1943E 04 

GI<2>« 0 . 23024777999S218E 03 

01(3)* -0.345565323S10732E 01 

GI<4>» 0 ♦ 3521876773671 16E 01 

GI<5>“ -0. 43S420216024420E-04 

GI<6>“ 0. 13463S34S013162E-07 

GI<7>» 0.16205982S959105E-10 

GI (8)“ 0.1 031 4685 1S72S6SE 01 

GI<9>* 0.223918105000000E 04 

GI(10>“GI<11>*0 

TO* 298. 15 % H0*8682 . * S0*205.037 

GI<10)»H0-HKT0> 

G<1)» - . 4308768468E-03 

G<2>» . 1979591095E+00 

G<3>» -.4143014968E+01 

G<4>“ .18536S4396E+03 

GlS>“ -.1270637452E+05 


G<6>« » 1S36388737E-04 
G<7>* .1326068945E-02 
G<8)“ -♦ 21992751 23E+01 
G<9>» .4705445127E+04 
G<10>“ ♦ 4728198017E-06 
G<11>- . 24304081 98E-02 
G<12>» - » 1 89675961 5E+00 
G< 13)“ - . 6887067207E-05 
G(14>“ -.61 32885 1 80E-03 
Q< 15)* -» 1836518694E+00 
IM16>® .257566387 IE-04 
0(17)“ - . 2415604646E-06 
G< 18)“ » 1438660831E-03 
G<19)“ -» 17039159B6E-05 
G(20)“ - . 23537059 1 7E+04 
G<21)“ -.2271707669E+06 
0<22)° -. 27538154/ 1E+02 
G<23)“ . 927764872£t+05 
G<24)“ -.411 4926856E-0 1 
G<25)“ . 1982233262E+01 


i 

i 

i 

i 
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B<26>« -.1239651142E-03 

0(27>« - • 6322588664E+00 

Q(28) n " * 2443207666E"07 

Q<29>« .1328704370E-04 

Q<30>« ~.1146313ei2E~09 

0(31)® 102U69308E-07 

0(32)® . 2334998237E-06 

01(11 >»8C~8I (TO) 

A( 1>® .8813947S3076E+02 

A( 2>» -♦ 490241 1 96 133E”01 

A( 3>« ♦ 168328893252E+01 

A( 4)® - » 32S161223390E+O2 

A( 5)® . 550300989872E+03 

A( 6>® - . 5109685061 15E+03 

A( 7)® ♦ 315091859049E+03 

A( 8)® - . 2325666592S8E+02 

A( 9>« - . 48842S479359E+02 

A(10>® ~. 150624217523E+03 

A ( 1 1 > ® .280441 60385 1E+03 

A(12)“ -. 176693896861E+03 

A(13)» .403247747449E+02 

A(14)« .25219868B365E+01 

A(15>» -♦ 136098316472L-01 

A(16>» .282316159403E+00 

A(17>® ~ .286645905341E+01 

A(18>® .617024212284E+01 

A(19>» -.810220795462E+00 

A (20)® -♦ 279601 068969E+00 

VP < 1 > ®7 » 568956 ♦ VP(2>«5.004836 

VP( 4) ®3. 454481 ♦ VP(5>®1.514 » 

VP(7>=54.359 * VP(8)»154.581 

VP(9)»0.0 

DTP® ♦ 4081997364372E+02 
DTPV® . 3318894767078E-03 


» VP(3>®-2. 137460 

VP (6)® .001 4606 


100 


FORMAT!# THE TEMPERATURE RANGE FOR OXYGEN 18 
17* WITH PRESSURES TO 1200 BAR#) 


54.359 TO 400 N* 


TCC®VP(8> 

PCC®VPN( TCC) #1.01325 
PTP«VP(6)#1 .01325 


rTP®VP(7> 
TUL«400 . 
TLL®TTP 
PUL® 1200. 
DCC»13.63 
RETURN * 


END 


SUBROUTINE DATA P H2 
DIMENSION G(32)»VP<9)»GI<11) 

DIMENSION GV<* > )»GT<9)»FV<4)»FT(4)»EV(8)»ET<8> 
DIMENSION A(2V) 

COMMON/SATC/A » DTPV 
COMMON/CPID/GI 

SSSK^"~ 

COMMON/DATAl/GV»GT»FVfFT»EV>ET 

COMMON/PARA/PERCEN 
COMMON/ ISF/N 
N®1 

GO TO 1 
ENTRY N H2 

N®2 

GO TO 1 
ENTRY 0 H2 


fil 
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N*J 

uu m i 

ENTR t Hi' 

N*4 

ou io t 

fcNfhir 1 Hi! 

N r -!j 

UONUNUH 

NO I KANlil'fJK' I PROPER TIES FOR HYUROUEN 
I HtKHOU YNAM1C PROPER i'lfe'B FROM 
K'OOFR A Nil MUUARIY XK' 7U-814 ANU 
XR /4~3b7 
R<°.0820t>AlA 

UArtMA«“ . 004,1 

u < 1> - 4.61438/7Sb6Sj4E-4 

0< 2) a 4.2331B4f>b608AE-2 
G( 3) *“£5 * 096!3b622Aij03E“l 

0< 4) 2 * 9230S9738269E+0 

0< b> 2 . 98760914721 1E+1 
G< A) a 1 ,803148601 4.tOE-b 
G< 7) ““l . 3222569S4639E-3 

(3( 8) a 3 , 016504431 701 E-l 
G< 9) « 5.093705560851E+1 

G(10> a 1 *97382832491 ye- 7 
G<11) a 2.858492039828E-4 
G(12) *-2.2282 79239 123E-2 

0(13) “-2.2574811367A4E-6 

G(14) « 2 » 414272369746E-5 

G<1S) *~1 ♦ A9S713398S88E-3 

G<1A) *-5 . 393676391275E-7 

G<17) a 3 . 998955244328E-9 
0<18> « 1 . 1424575A1274E-A 

“-1 • 25256622589AE-8 
G< a 0) *-4 . 917861934882E+ 1 

G(21> “-1.S85AAA0173A8E+2 

G<22) *-l ♦ 901602946272E-1 

G<23) a 9.1 980208 A2S00E+0 
G<24) *-3. 180455518810E-4 

B(25) a 1 . 1910S779192AE-3 
B<2A) *-3. 7913S277322SE“7 

3(27) ““3.983377A99095E-S 

?. ( 28) *-l » 234510854688E-10 

j<29) = 1.9502AA293499E-9 

3(30) ““2. 38034391 71 09E~ 13 

j<31) ““4.07357AA08192E-13 

*< 1>* » 91 661 77201 87E+02 

H 2>* - . 1 7949252444AE+00 

K 3>» . 454671 158395E+01 

'< 4)* ” . A58499S89788E+02 

K 5)* . 7344AA80453SE+03 

*< A)* - . 682501 045 175E+03 

»< 7>* .631 783A74710E+03 

i< 8)* - . 5394088 732B2E+03 

i( 9>« . 43092381 1 783E+03 

i(10)a - . 30029573801 1E+03 

<11)« .156567165346EF03 

i ( 12 ) a -.504103608225E+02 
<13)* . 720706926514E+01 

< 14)* - . 12394444031 BE +03 

(15) = . 140334800142E+01 

(16) * - .21 102380431 3E+02 

< 17)* . 1 7325462281 7E+03 

(IP)* ~ • *■ 44294S80871E+03 


Cr- -■ 


A <1V ♦ 1386993653551 +03 

A ( 20 > * - ♦ 2357 74 1 6 1 0 1 5E + 02 

U< 32 > « B ♦ 80U54930//7L ~ 12 

VP a >^3*05300134164 

VP ( 2 ) *2 ♦ 808 1 092581 3 

VP < 3 ) «— 0 ♦ 65546 3 216567 

VP<4>»1 #59514439374 

VP <?*)«* I . SB 14454428 

VP < 7 ) « 13*8 

VP <6 >*0*0693 

VP< 8 ) *32*938 

VP<9>*0 

DTP* *382 1428945 438E +02 
U f PV~ ♦ 6322296353698E • 0 1 
EM*2*01S94 
PRINT 300 

100 FORMAT OK THE TEMPERATURE RANGE FOR HYDROGEN 18 13 *H 10 400 R* 
1/* WITH PRESSURES TO 1200 BAR* ) 

TCC*VP<8> 

PCC-VPN < TCC ) * 1 .01325 
PTP«VP ( 6 ) *1 ♦ 0x325 
T TP«VP < 7 ) 

TUL«400 ♦ 

TLL--TTP 
PUL* 1200* 

DCC-15 ♦ 556 

RETURN 

END 

SUBROUTINE PROPS (PPfDDtTT) 

C THE 32 TERM EQUATION OF STATE ? INPUT IS DENSI TY (MOLb S/I.. ) * 

C TEMPERATURE < K > j OUTPUT <PP> IS PRESSURE IN ATM* OR DP/DU IN 
C LITER-ATM/MOLE OR DP/DT ATM/K OR SfHiOR CV AT ONE I iMJ .1 )+ 

C INTEGRATION 

DIMENSION X ( 33 ) 

DIMENSION &<33>>G<32) 

EQUIVALENCE <B*X> 

C0MM0N/DATA/6 * R > GAMMA 

COMMON/l/B 

0ATA<ID*1> 

DATA ( IZ-1 ) 

1 CONTINUE 

IF ( IZ ♦ L.E ♦ 0 )GU TU 2 
IZ*0 

2 CONTINUE 
0*00 

P« PP 
T*TT 

GM-GAMMA 
02*0*0 
03*1.12*11 
04*03*0 
05*04*0 
06*05*0 
07*06*0 
08*07*0 
09*08*0 
010*09*0 
011*010*0 
012*01 1*0 
013*012*0 
T8»8QRT < T ) 

T2«r*r 

T3*T2*T 



non 


T4«T3#T 
T5«T4*T 
F n EXP <GM*D2> 

GO TO <100*200*300*400*300*600*700) *K 
ENTRY PRESS 

c ENTRY FOR PRESSURE* INPUT IS DENSITY 

C AND TEHP * IN MOL/L AND K* OUTPUT IS IN AIM. 

K«1 

GO TO 1 
100 CONTINUE 
B< 1)bD 2*T 
B< 2)“D2*TS 
B< 3 )»D2 
B( 4>=D2/T 
B< 5) ra 02/T2 
B( 6>=D3*T 
B( 7)=D3 
B< 8)=D3/7 
B< 9)=D3/T2 
B<10)«D4#T 
B<11 )«D4 
B<12)=D4/T 
B< 13)=D5 
B<14)»D4/T 
B(15)»D6/T2 
B<14)»D7/T 
B<17)*D8/T 
B< 18)*D8/T2 
B<19)*D9/T2 
B<20)=D3*F/T2 
B(21)*D3#F/T3 
B<22)«D5*F/T2 
B(23)=DS#F/T4 
B<24)«D7#F/T2 
8<25>«D7*F/T3 
B<26)=D9#F/T2 
B<27)»D9*F/T4 
B<28)=D11#F/T2 
B<29)=D11*F/T3 
B<30)=D13*F/T2 
B<31)=D13*F/T3 
B<32)«D13#F/T4 
IF(ID*GT.O)GO 10 102 
B<33)-*P-R*D*T 
RETURN 
102 P®0 
M=32 

DO 101 1-1 >M 
101 P=P+B<I)*B<I) 

PeP+R*D*T 
PP*P 
RETURN 
ENTRY DPDD 

PARTIAL OF PRESSURE WITH RESPECT TO 
DENSITY - SEE PRESSURE 
ENTRY FOR UNITS 
K-2 

GO 70 1 
200 CONTINUE 

F1-2.00#F#GM*D 
F21=3.000#F*D2 H 1*D3 
F22»5,000*F#D4 +Fl*Db 
F23=7.000*F*D6 +F1*D7 
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n n o 


F24«9*000*F#D8 tFl*D9 

F25»l 1 ♦ 00*F»D10f Fl*01 1 

F26* : 13*00$F#D12+F1#D13 

D< 1)«2*00»0*T 

B< 2>«2.00*D*TS 

B< 3)<>,00*D 

B< 4)«2«00*0/T 

B< 5 ) *»2»00BD/T2 

B< 6>»3»00*D2*r 

B< 7>»3»00«D2 

B< 8>»3#00*D2/T 

B< 9>«3*00*lJ2/T2 

B< 10>«*4*00#D3#T 

B( 11 ) «4 *00*03 

B< 12 )»4 *00#D3/T 

B<13)*5*00*D4 

B < 14 ) ®6 ♦ 00BD5/T 

B<1S>*6*00#D5/T2 

B<16)*7*00#D6/I 

B<17)«8*00#D7/T 

B<18)«8*00*D7/T2 

B< 19 >«9 *00*08/12 

B<20)«F21/T2 

B<21>*F21/T3 

B(22)*F22/T2 

B<23)»F22/T4 

B<24)»F23/T2 

B<25 ) -F23/T3 

B<26)«F24/T2 

B(27)«F24/T4 

B<28>*F25/T2 

B<29)~F25/T3 

B<30)^F26/12 

B< 31 )«F26/T3 

B< 32 ) “F26/T4 

M*32 

IF < ID ♦ GT ♦ 0 ) 60 TO 202 
B<33)«P-R*T 
RETURN 
202 P*0 

DO 201 I*1>M 
201 P»P+B< I > 3|tl3 < t ) 

P«PFR*T 

pp«p 

RETURN 
ENTRY DPDT 

PARTIAL OF PRESSURE WITH RESPECT 
TO TEMPERATURE - SEE PRESSURE 
ENTRY FOR UNITS 
K*3 

GO TO 1 
300 CONTINUE 
X< 1 )»D2 

X< 2 ) »D2/ < 2 ♦ 00*TS ) 

X( 3)«0 
X< 4 ) ~~D2/T2 
X< 5 ) ~~2 ♦ 00*02/1 3 
X< 6 ) «D3 
X( 7)« 0 
X( 8)»~D3/T2 
X< 9>»~2#00*D3/TJ 
X<10)«D4 


n n 



x ( i :» > =- - u4/w 

X<13>=0 

x(i4>— u6/rr: 
x< ib)“-;*.oo*WA/| 3 

X< 16) =-07 / 12 
X<17>=-DU/T2 
x < ib > - ■ 2* oo*uu/ 1 ,1 
X < 19 )=-2*00*DV/ I i 
X(20;--2.00*D3*F/T3 
X<21)=-3.00*D3*F/T4 
X(22)*-2.00*D5*F/I3 
X < 23 > =~4 « 00#D5*F/ T 5 
X < 24 > = - 2 . 00*D7*F / 1 3 
X < 25 ) =~3 * 00*D7*p-/ 1 4 
X <26 > ~— 2 . 00*09*F/T3 
X<27)*-4.00*U9*F/T3 
X(28)=~2.00*D11*F/T3 
X<29)=-3.00*D11*F/T4 
X<30)*-2,00*D13*F/T3 
X(31)=-3,00*D13*F/T4 
X < 32 ) = -4 . 00*D1 3*F/T5 
1F< ID*GT *0)60 TO 302 
X<33)=PP-R*D 
RETURN 
302 P=0 

DO 301 l-lii 32 
301 P=P-FG< I >*X< I ) 

PP*P+R*D 
RETURN 
ENTRY DSDN 

PARTIAL OF ENTROPY WITH 
RESPECT TO THE 0 COEFFICIENTS 
K-4 

60 TO 1 
400 CONTINUE 

S=S0-R*L06F ( D*R*T/PO ) + ( DSDN < D ) -DSDN < 0 ) > * 1 0 1 . 325 +CP0S ( T ) 
Ol-F/ <2« 00*GM ) 

32* ( F*D2-2 ♦ 00*61 ) / ( 2 , 00*GM ) 

63* < F*D4-4 ♦ 00*62 ) / < 2 . 00*GM ) 

64* ( F#D6-6 . 00*63 ) / < 2 . 00*GM ) 

65* < F*D8-8 . 00*64 ) / < 2 ♦ 00*GM ) 

64= < F*D 1 0-1 0 ♦ 00*G5 ) / ( 2 . 00*6(1 ) 

X( 1 )=-D 

X< 2>=~0/<2.00*TS) 

X( 3)=0.D0 
X< 4)*+D/T2 
X< 5)=2*00*D/T3 
X< 6>=-D2/2«00 
X( 7)=0*D0 
X< 8)=D2/<2»00*T2) 

X( 9>=D2/T3 
X( 10)*-D3/3*00 
X< 11 )=0*D0 
X<12)*D3/<3.00*T2> 

X(13)*0.D0 
X < 1 4 ) =D5/ ( 5 ♦ 00*T2 > 

X< 15)* 2«00*D5/<5*00*T3) 

X<16)=D6/<6*00*T2> 

X< 17)=D7/(7*00*T2) 

X < 1 8 ) *2 . 00 * 07 / ( 7 ♦ 00*T3 > 

X<19)=D8/<4.00*T3> 

X<20)=2»00*G1/T3 
X (21 >=3.00*G1/T4 


fi£ 


& 


o 





Xt 22 >*2.00*02/ 1 3 
X < 23 )«4. 00*02/ 

X< 24 >-2.0041)3/1 3 
X ( 25 1*3* 00*63/ ( 4 
Xl26 ) * 2 * 00*04/ 13 
X< 27 >^4 .00*64/15 
X < 28 ><=2. 00*65/ T3 
X< 29 >*3.00*65/1 4 
X < 30 ) “2 ♦ 00*66/ 1 3 
X < 31 >*3.00*G4/T4 
X<32)«4.00*6A/T5 
IFUD.GT.OIGO III 402 
RETURN 

402 P*0 

DO 401 I»1 .32 

401 P=P+0(I)*XU) 

PP=P 
RETURN 
ENTRY DUDN 

TERMS NEEDED FOR ENTHALPY CALCULATION 

K=5 

GO TO 1 

500 hSho+<!*D»DN<D)-IWON<0) >*101.3284 <DUDN<D-0UDH«0> >*101. 325+CP0H<T) 
+(P/D-R*T>*101.325 
6 1 *F / < 2 ♦ 00*6M ) 

62 , »<F*D2--2.00*G1 >/<2.00*GM> 

03® < F*D4-4 ♦ 00*62 > / < 2 . 00*6M ) 

64= ( F4D6-4 ♦ 00*63 ) / ( 2 * 00*GM ) 

G5» ( F4D8-8 ♦ 00*G4 > / < 2 ♦ 004GM ) 

66* < F'*D10-10 *00*65 ) / ( 2 * 00*GM ) 

X< 1 )=D*T 
X< 2>=D*TS 
X< 3>*D 
X ( 4>*D/T 
X( 5) =D/T2 
X< 6>»D2*T/2.00 
X< 7>«D2/2.00 
X< 8)*D2/<2*00*T> 

X< 9 ) «D2/ ( 2 ♦ 00*T2 ) 

X(10)«D3*T/3.00 

X(ll)«it3/3.00 

X < 1 2 > *D3/ ( 3 ♦ 00*T ) 

X<13>=D4/4,00 
X ( 14 ) «D5/ < 5 ♦ 00*T ) 

X<15>=D5/<5.00*T2> 

X<16)»D6/<6*00*T> 

X<17)«D7/<7*00*T> 

X( 18)®D7/<7 *00»T2> 

X < 19 ) =08/ < 8 ♦ 00*T2 ) 

X<20>*G1/T2 
X(21>=G1/T3 
X(22)“G2/T2 
X<23>=G2/T4 
X<24>*G3/T2 
X< 25 >-63/13 
X(26>*G4/T2 
X<27>»G4/T4 
X<28)=G5/T2 
X<29>-65/T3 
X<30>-G6/T2 
X(31)«G6/T3 


X<32>*G6/T4 



n n n on 


IF < ID«QT*0>00 TO 502 
RETURN 
S02 P-0 

DO 501 1-1*32 
501 -P+Q<I)*X<I> 

PP-P 
RETURN 
ENTRY Tl'SDT 

TEMP* TIMES I HE PARTIAL OF 
ENTROPY WITH RESPECT TO TEMP. 

K-6 

00 TO 1 

600 CONTINUE 

CV-CV0+ < TDSDN ( / ) -TDSDN < D ) ) * 1 0 1 . 325 
G1«F/<2.00*GM) 

02- < F*D2-2 . 00*G1 ) / < 2 . 00*GM ) 

03- < F*D4-4 . 00*02 ) / < 2 . 00*GM ) 

04- <F*D6-6 *00*03 )/(2*00*GM) 

05- < F*D8-8 ♦ 00*04 ) / < 2 ♦ 00*GM ) 

06- < F*D 1 0- 1 0 . 00*05 > / < 2 . 00*GM ) 
X<1)»0 

X< 2>— D/<4.00*TS) 

X<3>»0 

X< 4)>2.00*D/T2 

X< 5)«6*00»D/T3 

X<6)-0 

X<7)-0 

X< 8)»D2/T2 

X< 9 ) -3 ♦ 00*D2/T3 

X<10>»0 

X<11)»0 

X<12)-<2.00*D3)/(3.00*T2) 

X< 13)«0 

X<14)-<2* 00*D5 ) / ( 5 ♦ 00*T2 ) 

X(15)-(6* 00«D5 > / ( 5 ♦ 00*T3 ) 
X(16)-D6/<3*00*T2> 

X(17)-<2* 00*D7 ) / < 7 ♦ 00*T2 ) 

X(18)-<6* 00*D7 )/<7*00*T3) 
X(19)-C3*00*D8>/(4.00*T3> 
X(20)-6*000«01/T3 
X<21)»12.00*G1/T/ 

X C 22 ) =6 ♦ 000*G2/T3 
X < 23) -20 « 00*G2/T5 
X < 24 ) -6 * 000*G3/T3 
X< 25) =12.00*03/14 
X ( 26) -6 « 000»04/T3 
X ( 27) -20 « 00*G4/TS 
X ( 28 ) -6 « 000*05/ T3 
X < 29 ) -12 « 00*05/T 4 
X< 30) -6 .000*06/13 
X(31 )-12«00*06/T4 
X ( 32 ) -20 » 00«G6/T5 
IFdD.OT *0)00 TO 602 
RETURN 
602 P-0 

DO 601 1-1*32 

601 P»P+0<I)*X<I) 

PP-P 
RETURN 
ENTRY DP2D2 

SECOND PARTIAL OF PRESSURE WITH 
RESPECT TO DENSITY SQUARED 
K-7 


n o non n n 


oo ro 1 { 

700 CONTINUE 

Fl«2.*F*GM*D * 

F 12*2 . *F l*UM*0+2 < *F*GM 4 

F212»3.*F1*D2+3.*2.*D*FFF12*D3FF1*3.*D2 « 

F222«=5.*F1*D4 +S.*4.*D3*PF5.*U4*F 1+F12*D5 f 

F232-7 . *F l*D6f 7 . *6 , *DS*F+7 . *DA*F 1+F 1 2*D f \ 

F242»9 , *Fl*D8+9 . *8 . *D7*F+9 . *D8*F 1FF12*D9 * 

F252«»11.*F1*D10+10.*11.*D9*F+11.*D10*F1+F12*D11 
[*262*13. *F1*D12+13. *12. *D11*FF13.*D12*F1+F12*D13 
B< 1 )“2«*T ♦B(2)“2.*T8 * B<3>*2. 

B(4)»2./T ♦ B(5)»2*/T2 * B<6)»6.*D*T 

B<7)*6»*D ♦ B(8)-4.*D/T 4 B(9)®6.*D/T2 

B( 10)=12.*D2*T * B< 11 )*12»*D2 * B< 12>*12.*D2/T 

B(13)»20.*D3 * B< 14>«30.*D4/T ♦ B< 15)=30.*D4/T2 

B(16)»42.*DS/T * B<17)»56.*D6/T * B< 18>»56.*D4/T2 

B(19)®72.*D7/T2 * B<20)-F212/T2 * B<21 >®F212/T3 

B ( 22 ) “F222/T2 

B(23)=F222/T4 * B(24)*F232/T2 ♦ B<25)*F232/T3 

B(26)*F242/T2 * B(27)*F242/T4 * B(28)*F252/T2 

B(29)*F252/T3 ♦ B< 30 )=F262/T2 * B<31 >*F262/T3 

B(32)*F262/T4 
M*32 

IF < 1D.0T .0)00 TO 702 
B(33)*PP 
RETURN 
702 P*0 

DO 701 1*1 »M 

701 P*P+B< I )*6< I ) 

PP*P 
RETURN 
END 

FUNCTION VPN(TT) 

CALCULATES VAPOR PRESSURE IN ATMOSPHERES 
FOR AN INPUT TEMPERATURE IN KELVIN 
DIMENSION 0(32) f VP (9) 

COMMON/DAT A/G t R » GAMMA t VP 
T*TT 

X* < 1 ♦ -VP (7)/T)/(l. ~VP < 7 ) /VP ( 8 ) ) 

VPN*VP < 6 > *EXP < VP ( 1 ) *X+VP ( 2 ) *X*X+VP < 3 ) *X**3+VP < 9 ) *X**4+VP < 4 ) *X* 

1 < 1 * -X ) **VP ( 5 > ) 

RETURN 
END 

FUNCTION FINDTV(POBS) 

ITTERATES THE VAPOR PRESSURE EQUATION 
FOR A TEMPERATURE < IN KELVIN) 

GIVEN AN INPUT PRESSURE IN ATMOSPHERES 
CQMMON/DATA/G * R » GAMMA » VP t DTP 
DIMENSION G(32)»VP<9) 

T*VP<8) 

DO 7 1*1 t 10 
P«VPN(T) 

IF < ABS < P-POBS ) ~ • 00000 1 *POBS > 8 » 8 » 6 

6 CONTINUE 

CORR= < POBS-P ) /DPDTVP < T ) 

7 T*T+CORR 

8 CONTINUE 
FINDTV*T 
RETURN 
END 

FUNCTION CV < D » T ) 

CALCULATES SPECIFIC HEAT CAPACITY 
AT CONSTANT VOLUME FOR AN INPUT 
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c 


c 

c 

c 

c 

c 


OF DENSITY AND TEMPERATURE IN MOL/l- AND In 
DATA<R‘*8. 31434) 

DD-'D 

TT»T 

CALL TDSDf <CD*DD*TT> 

DD«0 

CALL I D8DT < CO • DD»T1 ) 

CV-CPK T T > ♦ < CO-CD > # 1 0 1 * 323 

cv*cv-r 

RETURN 


END 

FUNCTION FIND D<P»T> 

ITTERATE8 EQUATION OF 8 TAIL 
FOR DENSITY * GIVEN PRESSURE 
AND TEMP. IN ATM. AND KELVIN. 
I ITERATION FAILS TRY USING 
f UNCTION CALLED FIND M 
DIMENSION G<32)*VP<9> 
COMMON/DATA/G * R * GAMMA f VP * D TP * 
TT«T 

IF<TT.GI.VP<8))G0 TO 100 
IF < P.GT .VF'N(TT) )G0 TO 101 


IF 


PCC * PI P * TCC » T T P f TUL r PUL * DCC 


DD*SATV(TT) 

GO TO 102 

100 PC*PCC/ 1.01325 

X»<l.l/<9.*PC))*P+.7/y. 

DD»P/<R#1*X) 

IF(P/PC»GT .20. ,AND.T/VP<8).LT.2.5)DD»DTP 


GO TO 102 

101 DD“*3ATL<TT) 

102 CONTINUE 

DO 10 I-lrSO 

IF(DD.LE.0.0.0R.DD.GT.5O. )G0 TO 11 
CALL PRESS <PP*DD*TT> 
IF(F‘P.LE.0.0)G0 TO 11 
P2=PP 

IF ( ABS < P-P2 ) -1 . E-7*P >20*20*1 
1 CALL OPDD < PP * DD » TT > 


C 


DP-PP 

CORR“ <P2-P>/DF‘ 

IF < ADS <C0RR)-1.E-7*DD>20*20*10 
xO DD=DD-CORR 
11 CALL REGULA<P*DD»T> 

20 FIND D*DD 
RETURN 


END 

SUBROUTINE REGULA < PI » DD * TT > _ ,, _ 

ITTERATES EQUATION OF STATE FOR DENSITY WHEN FIND D FAILS 

DIMENSION 8<32>*VP<9) 

COMMON/DATA/G * R * GAMMA * VP » DTP * PCC * PTP » TCC * T TP » 1 UL. * 1 LL * PUL * DCC 


T“TT 

P«PI 


D2»0 

IFIT.LT .TCC)GO TO 10 
D0»DCC*TCC/T 
GO TO 20 

10 PP=VPN(T) 
IF(P.GT.PP)GO TO 13 
D0“3ATV< T ) 

DO 11 1=1*150 
CALL F’RESS ( PO * DO * T ) 
IF (PO.GE.P)GO TO 12 

11 DO^DOE .0001*00 


C?* - — * - — - v 
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UP 10 42 

12 !U*liO 

13 CALL PKESS<P1>01>T) 
iF<Pl.Lr.P)U(l I Cl 14 
XF<D1.LE..1#PTP>00 10 42 
00*01 

2«<P1~P)/P 

if<z.lt..i>z*.i 

IF<Z.GI . .9)2*.? 

01 * 01 - 2*01 

00 TO 13 

14 CAt L PRESS <P0> DO >1 ) 

00 140 1*1 >50 

0*01 

P3*P1 

IF < AB8<P“P1 ) .LT ♦ ♦ 00001 *F') GO TO 40 
P2-P-P1 

01*01 + <01-00 )*P2/< PI ~P0> 

IF < AOS <0-01 ) .LE. « 00001*0)00 TO 40 
I F < ABB < P-Pl ) 4 LT ♦ ♦ OOS*P ) D2*F I NO M<P>1 >01) 
IF <02. GT .0.0. AND. 02 .LT. 50 .>01*02 
02*0 

CALL PRESS < PI >01 fT) 
IF<PQ.GT.P.AND.P1.GT.P)U0 TO 120 
IF<PO.L r.P.AND.Pl.LT.P)GO TO 120 
GO TO 140 
120 P0*P3 
OC~D 

140 CONTINUE 
GO TO 41 

15 DO*SATL< T) 

DO 16 I*l>10 

CALL PRESS <P0>00>T) 

IF<PO.LE.P)DQ TO 17 

16 00*00-. 0001*00 
GO TO 42 

17 01*00 

18 CALL PRESS < PI >01 >T) 

IF <01 .GE .50. )GQ TO 42 
IF<P1.GT.P)G0 TO 14 
00*01 

Z«<P-P1)/P 

Z«Z*10 

IF ( T/TCC .LT ♦ .6)2*1 ♦ 

IF(Z.LT . 1 ♦ )Z*1 ♦ 

IF<Z.GT.9.)Z*9. 

01 * 01 + . 01 * 01*2 

GO TO 18 

20 CALL PRESS<P0>00>T) 

IF<P.LE.PO)GO TO 30 
01*00 

21 CALL PRESS<P1 >01 >T ) 

IF<P1.GE.P)G0 TO 14 
IF<01.GE.50.)G0 TO 42 
00*01 

Z*<P~P1)/P 

Z«Z*10 

IF < Z.LT ♦ 1 ) Z*1 
IF<Z.GT.9)Z*9 
D1*D1+.1*01*Z 
GO TO 21 

30 01*00 

31 CALL PRESS < PI >01 >T) 
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if<pi.le.p>oq to 14 
iftdi.le. .i#ptp>qq ru 4? 

D0»D1 

Z*<P1-P>/P 

z*z#io 

iF<Z.l.T.nZ«l 
IF<2.GT.?>Z*'9 
D1»DI-,1#D1*2 
00 TO 31 

40 DD»D1 
RETURN 

41 PRINT 101»P»T/D 

102 FORMAT <$ REOULA FAILED AT P»#>r /.2»* AND T»*»F7.2> 

101 FORMAT <# DENSITY 1TTERAT ION FAILED AT P«=*»F7.2»# AND T»#»F/,2» 
1/* DENSITY RETURNED I8#»E17.8> 

RETURN 

42 PRINT 102fP»T 
RETURN 

END 

FUNCTION CP < D f T > 

CALCULATES SPECIFIC HEAT CAPACITY 
AT CONSTANT PRESSURE FOR INPUT OF 
DENSITY AND TEMPERATURE IN MOL/L AND K 
CP IS IN JOULES/MOL -K 
CVEE-CV<D»T> 

CALL DPDT <DPr»D»T> 

CALL DPDD<DPD»D»T> 

CP=CVEE+ < T/ ( D*#2 > * < DPT**2 > /DPD t * 1 0 1 ♦ 325 

RETURN 

END 

FUNCTION DPDTVP(TT) 

CALCULATES THE DERIVATIVE OF PRESSURE 
WITH RESPECT TO TEMPERATURE AT 

SATURATION. INPUT IS TEMP. IN K» OUTPUT IS ATM/K. 

COMMON/DATA/G , R t GAMMA » VP 
DIMENSION <3<32>*VP<9> 

T=TT 

IF(TT.GT.VP<8) )G0 TO 1 

X*U.~VP<7)/T>/<1.-VP<7>/VP<8>> 

DXDT=<VP(7>/T**2>/<1.~VP<7>/VP<8>> 

DPDT»VP ( 1 ) $DXDT+2 ♦ *VP ( 2 ) *X*DXDT+VP ( 3 ) *3 . *X#*2*DXDT+VP < 4 ) # 

1 < ( 1 . -X ) **VP < 5 ) ) KtDXDT +VP < 4 > *X# < < 1 ♦ -X > ## < VP < 5 > - 1 . > ) 3RVP ( 5 > * < -DXDT > 
DPDT =DPDT#VPN < T ) 

DPDTVP»DPDT 
RETURN 
1 DPDTVP*0 
RETURN 
END 

FUNCTION FIND M<P»T»DD) 

ALTERNATIVE FOR FIND D» INPUT IS 

PRESSURE IN ATM. » T IN KELVIN AND 

DENSITY IN MOL/L. INPUT DENSITY 

IS A STARTING VALUE FOR ITTERATION 

OF EQUATION OF STATE FOR SOLUTION FOR P AND T 

TT“T 

DO 10 I=*l >50 

CALL PRESS (PPtDDfTT) 

P2“PP 

IF (ADS <P-P2)-1.E~7#P)20»20>1 
1 CALL DPDD<PP»DD»TT) 

DP»PP 

CORR* < P2-P ) /DP 
D=DD 
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If (ABB <CORR>M»E /*p>20>?0*10 
10 PP'~PD~C0RR 
MNP M O 
HE TURN 

:»•) MNP M* n PP 
NT TURN 
f NP 

FUNCTION tNTHALCPfDfT) 
n CALCULATES ENTHALPY FOR INPUT OF 
0 PRT SSURE r DENSITY ANP TEMP# IN 
C ATM* * MOL/l ANP K* UUIfU! IB IN 

0 JQULES/MUL 

R ***0820fi416 
PPM* 

TT T 

CALL DSDNCSDtDDt f f ) 

CALL DU»N(UUfPPr TT ) 

DD^O 

CALL UPON < SO > PH t TT) 

CALL PUDN< UO»PP *T T > 

LN rHAL «T* < 8D-80 ; * 101 ♦ 325+ < UD-UO ) *101 ♦ 325+HI < T ) + ( P/D-R*T ) *101 ♦ 325 

RETURN 

END 

FUNCTION ENTRQP< Dt X ) 

C CALCULATES ENTROPY 

C FOR AN INPUT OF DENSITY ANP 

TEMP* IN MOL/L ANP K* OUTPUT IS IN 
JQULES/MOL-K 
R ==*08205616 
DD«D 
TT-T 

CALL DSDN<SDrDD> TT) 

DP«0 

CALL DSDN<SO>DDtTT> 

ENTROP* <SD~SO> *101 * 325-R*AL0G < D*R*T ) Ht 1 0 1 * 325+8X1 D 
RETURN 
END 

FUNCTION SATL(TT) 

CALCULATES THE DENSITY OF THE 
SATURATED LIQUID AT TEMP* $ T IN KELVIN* 

OUTPUT IS IN MOL/L* 

DIMENSION A (20) 

DIMENSION C<32) *VP<9) 

COMMON /DAT A/0 > R t GAMMA > VP > DTP * PCC rPTP > TCC > TTP > TUL r TLL * PUL > DCC 
COMMON/SATC/A t OTPV 
ls>14 
KK~7 

GO TO 10 
ENTRY SATV 
K~1 
KKM3 

10 IF(T*GE*TCC)GO TO 20 
T«f T 
t i r*rrr 

IF* ( ITT+l-T *LT ♦ 1 ♦ >T«ITT 
X^( r~TCC)/< T TP-TCC) 

D*A<IO*ALOG<X> 

DO 11 I«2fKK 
K S; K+1 
MM - 1 

IF <MM*GE*5)MM^MM+1 
1 1 D-D+A < K ) * < 1 ♦ --'AM ( < MM-5 > /3 ♦ ) ) 

IF <K*LT*14)G0 TO 12 
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D ’ICCPKXP < U > *< P I P-DCC ) 

GO ro 1.4 

i ? n pcc ffxp < o ) « < dtpv-dcc: ) 

13 SATL"D 

l r < irr n - t r . lt , i . > s a r i - u- <w ucc >*ci r-r> 

PE TURN 
.TO U8ATL UCC 
RETURN 

ENll 

FUNCTION 80UND< D* l > 

CALCULA I E8 THE SPEED OF SOUND 
FOR AN INPUT (.IF DENSITY AND TEMP. 

IN MOt/L AND KELVIN, OUTPUT IS IN 
METERS/SECOND. 

COMMON/CRIT/W 
CALL DPTHl ( DP r D r T > 

SOUND" < ( UP < D , r ) /CV ( 0 » r > > #DP* 1 01320 . /W ) #* , 8 
RETURN 
END 

FUNCTION VISC<DD»T> 

RETURNS VISCOSITY IN <G/CM--S)*E+6» 

T IN K» D IN MOL/L 
COMMON/CRIT/EM 
D-DDfcEM/lOOO , 

VISC"D TLV(T) +FDCV < D t T ) +EXCESV < D » T > 

RETURN 
END 

FUNCTION THERM <DD»T ) 

C RETURNS TC IN MW/M-K.- T IN K, D IN MUL/L 
COMMON/HAN/CR » TCI 
COMMON/ I SP/N . NU 
COMMON/CRIT/EM 
D"DD*EM/1000, 

IF(NU.EQ.O > GO TO 3 
CR*CRITC(D»T) 

THER*DILT < T ) +FDCT < D » T > *1 00 . +EXCEST ( D , T ) +CR 
TCI=THER~CR 
THERM"THER 
RETURN 

3 CR~CRITC(D»T> 

THERM«DILT < T ) +FDCT ( D » T ) +EXCEST ( D t T ) +CR 

TCI=THERM-CR 

RETURN 

END 

FUNCTION EXCESV(D»T> 

C CALCULATES EXCESS VISCOSITY 

C0MM0N/DATA1/GV >BT»FV»FT »EV»ET 
COMMON/ I SP/N»NW 

DIMENSION GV<9)»GT(9)*FV<4>»FT<4>»EV<8>»Ei<8) 

R2"D##< ,5)#<<D~EV<B)>/EV<8>> 

R"D##< .1) 

X=EV<mEV<2>#R2+EV<3>*R+EV<4>#R2/<T#T>+EV<5)#R/T**<1.5>+EV<<S>/T 
1+EV(7)*R2/T 
X1"EV< 1 ) +EV(4)/T 
EXCESV=EXP (X)-EXF* (XI) 

RETURN 

ENTRY EXCE8T 

C CALCULATES EXCESS THERMAL CONDUCTIVITY 
IF<NW.E0.0 ) GO TO 3 
R»D##< ,1 ) 

X“ET< 1 ) +ET < 2 ) *R+ET < 3 ) *R/T*# (1.5) f ET ( 4 ) /T 

X1=ET(1)+ET(4)/T 

EXCESV«(EXP (X)-EXP <X1))/10. 

RFTURN 
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i PI = 1 1 + + t * b )*< ( U-t I ( 0 > ) /X T (II) ) 

K“U$1< , I ) 

X =1- I ( I ) H: r < ..')*R2+K | ( s ) *R Ft It'D *R2/ <( * I > It I ( u ) *R/ I ** < l . b ) ft r ( A ) / T 
i+t r< />*R2 / 1 

X I L I < 1 ) I t. I ( A ) / | 

EXI.ECV’EXI' (X)-fcXP (XI) 

RE I UK’N 
END 

I I INCH. UN IOCVUI.I) 

1 1 ivy i ntNS 1. 1 y cokrec r i on 

l-o R VISCOSITY ANU I HblVMAL l.'UNUUC I I Oil Y 
COMMUN/DA'I A t. /(IV . U I *1 V » I- T . F V*E I 
DIMENSION GV<9) »m<V) »l V<4> »F I (4) »EV<8> »EI (1.1) 

f dcv® < rv a ) it v < 2 > * < rv < n > auig < 1 /fv < 4 >>>**;> > *u 

Rfc I URN 
FN I IVY FOC I 

n DCV •» < FT < 1 > +FT < 2 ) * < F I < 3 > -ALGO < I /F t < 4 n ) **2 > *D 
RE I URN 
ENri 

FUNCTION UKirC(D»r) 

0 CALCULATES CRITICAL ENHANCEMEN I 

C FOR THERMAL CONDUCIIVITY 

COMMON/CR] r/ EM » EOK* RM» TC» DC» X » PC» S1G 
COMMON/ CHECK/DELD » BELT t DSTAR f TSTAR 
COMMON/H JM/EF'S 1 rCPCV » RRR , AK T 
COMMON/ ISP/N f NW r NWW 
C 0 IN G/CM3 f I IN K 
(. OUTPUT UNITS ARE MW/M.K 

AV=6 ♦ 0225E +23 $ BK«1 .38054E-16 
DELD-ABS <D -DO/DC * DELT“ABS (T-TO/TC 
C CALCULATE DISTANCE PARAMETER 

R= ( RM**2 . 5 ) * < D**0 . 5 ) * ( AV/EM ) **0 , 5 
R=R* < E0K**0 . 5 > *X/ < T**0 . 5 ) 

RRR»R 

C GENERAL EUUAT10N 
DX-D*1 OOO.O/EM 
CALL DPDT (DPT »DXf T) 

C DPDT IN ATS PER DEG. 

DPT”DPT*1 .01325E4A 
C DPDT NOW IN DYNES PRR DEG 
CALL DPDD ( DPD t DX t T ) 

C DPDD UN ATS» MOL/L 

DPD~DPD# < 1 . 01 325E+6 ) * ( 1 000 . 0/EM > 

CDPDD NOW IN DYNES t GM/CM3 

IF ( DPD. LT, 0.0) DPD=1 . 0 
94 VIS*VISC<DX»r>#< 1 . 0 E- 06 ) 

C VISCOSITY IN GM/CM.S 

IF < DEuD , EQ ♦ 0 ♦ 25 . OR, DELD.LT , 0.25) 8>10 

8 1F<DELT.E0. 0.025, OR. DELT.LT, 0.025 ) 9r 10 

9 COMPRES=SENG(Df D 
GO TO 12 

10 COMPRESS 1 , 0/ ( D#DPD ) **0 . 5 
12 E X«BK*T**2* ( DP ) #COMPRES 

EXB”R* < ( BK*T > #*0 . 5 ) * ( D**0 ♦ 5 ) * < < AV/EM ) **0 ♦ 5 > 

CRI T”EX/ ( EXB*A . 0*3 , 1 41 59* V I S ) 

C PUT IN DAMPING FACTOR 
BDD= < ( D-DC ) /DC ) **4 
DTT«< (T-TC)/TC)**2 
FACT«EXP ( ~ 18 ♦ 66*BTT - 4.25*BDD> 

DELC"CRIT*F ACT/ 100.0 
CRTTC-DELC 
Ah r =COMPRES*COMPRES 
EPSI~R*R*BK*T* ( AV*D/EM ) *AKT 
EPS t ”EPS I **0.5 
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CALC CP-CV 

CF'CV^T* ( DPT ##2 ) #AN T/D 
RETURN 
END 

FUNCTION 8EN0<D»T) 

COMNqE/c?i?/ 1 ?h °LS Trt e r £ 1=0,5 f-^riCAL REGION 
LunnUN/CRIT/ EMu EQK» RM* TC» DC# X • pc . urn 

CQMMON/fJEN/RETA > XO r DEL TA r El t F,2, ^fiAM * U BK 
COMMON/CHECK/DELD» DEL T » D8IAR » f 8 IAR 
DbTAR- D/DC » T8TAR™ T/TC 
DETOn./PETA 
XX®DELT/PELD#*BF TO 
AG®AQAH"1 .0 
BET2'' 2,0* BE TA 
AGB‘ AG/DE » 2 
DEL1 'jDEIIA~ 1 ,0 
AOBB- ( AG -BE T 2 ) /BET 2 
,VXOt.<XX+ XO)/XO 
XXB=«XX0**BET2 
»RAK«1.0 + E2*XXB 
BRAK1*BRAK**AGB 
H“E1*XX0*BRAN1 

RCOh« '( WLD* J nJi ^2u X ° J * E1 * E2 * < XX B > * < BRAK**AGBB > 
rsuun \ LitLUffDEUl )# ( DELTA#H - < XX/BETA i iUpprH \ 

RCOMP= 1 « 0/ < RCOM*DST AR**2 ) < **'«ETA>*HFKIM ) 

RCM=RC0MP/ ( PC*1 , 0 1 32SE+06 ) 

RCM IN DYNES 

RCM«RCM**0,S 

SENG«RCM 
RE TURN 
END 

FUNCTION DILV(T) 

GIVES DILUTE GAS VISCOSITY AND THERMAi 
CONDUCTIVITY FOR AN INPUT TEMP. ?S 
OUTPUT UNITS ARE SAME AS 
THOSE IN VISC AND THERM 

COMMON/ I SP/NfNWfNWW 
COMMON/DATAl/GV»GTi>FV,FT»EV»ET 
SUM=S SI0N 0V<9> » eT <’>»™<4>»FT<4>,EV<8>,ET<8> 
TF=T**<l,/3. > 

TFF«T**<~4./3.) 

DO 10 1*1,9 
TFF»TFF*TF 

10 SUM^SUM+GV < I ) *TFF 
IF(NWW*EG»7> GO TO 9 
DILV*SUM*1000. 

GO TO 11 
9 DILV»SUM 

11 RETURN 
ENIRY DIET 
TF«T#*<l*/3.) 

TFF=T**<~4,/3. ) 

SUM'=0 

DO 20 1*1,9 

TFF»TFF*TK 
20 SUM*8UMM3T < I )*TFF 
DILVaSUM 
RETURN 
END 

FUNCTION CPI \ T ) 

CALCULATES IDEAL GAS THERMO PROPERTIES 
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C FOR ALL FLUIDS EXCEPT H2. INPUT IS IN 

C KELVIN* OUTPUT IS IN JOULES* MOL/L AND K 

CQMM0N/CPID/G<11> 

COMMON/ I SP/N 
IF(N.NE.O)OQ TO 5 
K-l 

1 U=G(9)/T 
EU=EXP <U) 

TS=1./T*#4 
60 TO (2*3*4), K 

2 CPI=G(8)#U*U*EU/(EU-1. >##2 
DO 10 1=1*7 
TS=TS#T 

10 CPI=CPI+G<I)*TS 
CP I =CP I #8. 31434 
RETURN 

5 CPI=CPO(T *N) 

RETURN 
ENTRY SI 

IF (N.NE.O)GQ TO 6 
K=2 

60 TQ i 

3 CPI=G(«>*vU/(EU-l. >-ALOGU.-l./EU>> 

l-G(l>*TS*T/3.-G(2)*TS*T*T/2.-G(3>/T-FG(4>*AL0G<T)+G(5)*T+6(6)*T#T/2 

2.+G(7)#T**3/3. 

CPI=CPI*8. 31434+G (11) 

RETURN 

6 CPI=CPOS<T *N> 

RETURN 
ENTR HI 

IF < N . NE * 0 ) GO TO 7 
K=3 

GO TO 1 

4 CPI=G(8)*U#T/(EU-1 ♦ )-G( 1 >/(2.#T#I >~G<2)/T+G(3>#ALGG(T>+G(4)#T 
1 +G ( 5 ) *T*T/2 ♦ -FG ( 6 > *T**3/3 . +6 ( 7 ) #T##4/4 . 

CP ■ -=CPI#8 . 31434+G ( 10) 

RETURN 

7 CPI=CPOH(T *N> 

RETURN 
END 

FUNCTION lPO(TI*N> 

DIMENSION T (58 ) »CPP<58> *CPN(58) *CP0<58) *CPE(58) 

COMMON/PARA/PERCEN 

CALCULATES IDESL GAS SPECIFIC HEAT TOR H2 BY INTERP0LA1 IPG 
DATA TAKEN FROM RP 1932* UNITS OF THE TABLES ARE CAL/MOL DEG 5. 

UNITS OF OUTPUT ARE JOULES/MOL DEG K. THE INDEX N DETERMINES THE 

SPECIES*FOR N=1 *PARAHYDR0GEN*N=2 NORMAL* N=3 ORTHO* N=4 EGUILIB 
C N=5 * SOME ORTHO-PARA FIXTURE SPECIFIED BY COMMON /PARA/ * PERCENT 

C RANGE OF TEMP IS FROM 10 TO 5000K. 

DAI A(T= ^ ^ _ , 

' 1 10.0* 12.0* 14.0* 16.0* 18.0* 20.0* 25.0* 30.0* 35.0* 40.0* 45.0* 

* 7 50.0* 55.0* 60.0* 65.0* 70.0* 75.0* 80.0* 85.0* 90.0* 95.0*100.0* 

3105.0*110.0*115.0*120,0*125.0*130.0*135.0* 140.0*145.0*150.0*160.0* 
41 >0,0* 180.0* 190.0*200. 0*210.0*220. 0*230. 0*240. 0*250. 0*260. 0*270.0* 
5280 . 0 * 290 . 0 * 300 .0*350.0* 400.0*500 .0*600 ♦ *700 . * 1000 , * 1500 . * 2000 . » 
63000 . * 4000 . * 5000 . ) 

; DATA( <CPt< I ) * I - 1 *58 >=4.968 ,4.96884*4.97647*5.01 153*5.0/451 *5.208 

: 11 *5. 83508*6. 81282*7. 87989*8. 60613*9, 00231 *9. 08005*8. 93278*8. 65894* 

: 28 . 33603 *8,0120/*/. 71009 ,7. 4416*/. 2 11 09 /. 01 858 , 6 . 85857 , 6 . 72557 ,6.6 

32055 * 6 . 53555 * 6 . 46904 * 6 . 42003 , 6 . 38403 ,6.36151*6. 34602 *6.33 753 * 6 . 34C 
401 *6. 345// *6, 372/6*6. 413*6. 45925*6 .50975*6 .5605*6. 6195* 6. 65/24* 

V / 6.692*6.734*6. 7/1*6.804*6.832* 

■: 56. 856, 6.877*6. 895*6. 950*6. 9/4 *6. 993 *7.009* 7.036* 7. 219* 7. 720*8. 195* 
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68,859,9.312,9. /48> 

DA!A<CPO- 

14.965,4.968,4.968,4.968,4.968,4.968,4.968,4.968,4.968,4.968,4.968, 

24 , 968 » 4*968*4. 960 y 4 . 969, 4 , 9/2, 4 » 9/5 , 4 » 981 y 4 • 990 f Si , 002 » 5 .018,5.039, 
35 . 064 f 5 * 0y4 , 5 . 129 ,5 ♦ 169 , 5 * 2 1 3 , & ♦ 261 » 5 • 3 1 3 , S ♦ 369 , 5 . 42/ ,5 . 487 , 5 . 6 1 2 » 
45. 741 ,5. 868 y 5 . 992 1 6 ♦ 1 09 r 6 , 2 19 , 6 , 320 » 6 ♦ 4 1 1 * 6 . 493 » 6 . 566 , 6 . 629 * 6 . 684 , 
56 * ? S.l ,6.773,6. 808 , 6 « 91 7 , 6 . 962 »6 • 993 , 7 ♦ 009 » 7 . 036 , "7 , 219 » 7 » 720 ,8.195, 
68 « 059 , 9 ♦ 342,9 * 748 ) 

DATA'.CPP- 

1 4 . 9o8 y 4 , 968 , 4 . 968 , 4 . 968 » 4 . 968 , 4 . 968 , 4 . 968 1 4 . 968 1 4 . 969 , 4 . 972 , 4 . 983 » 
23 . 006 , 5 . 040 » 5 ♦ 1 14 ,5 • 207 » & ♦ 328 ,5.475,5. 646 , 5 ♦ 835 ,6.036,6. 245 » 6 ♦ 454 , 
.56*659 , 6 .854, 7 .037, 7. 203,7 ♦ 351 » 7.480,7 . 590 » 7 .681,7. 753 , 7 , 807, 7.870, 
47 • 883,7 « 858,7 « 800 ,7. 742, 7.667 ,7. 591 ,7*516, 7.445, 7. 380 r 7. 322, 7 . 270 , 
57.225,7.186,7. 152 » 7. 050 ,7. 010, 6.998 , 7.010 , 7.037 , 7.219 » 7. 720, 8.159 , 
68. 859 , 9 .342*9.748) 

DATA<CPN~ 

14 . 968 , 4 ♦ 968 , 4 . 968 » 4 . 968 » 4 ♦ 968 , 4 ♦ 968 , 4 . 968 , A . 968 , 4 . 968 1 A ♦ 969 » 4 ♦ 972 , 
24. 977 ,4. 988 >5. 005 ,5. 029 ,5. 061, 5. 100 ,5. 147, 5. 201, 5. 261, 5. 325, 5. 393, 
35, 463, 3.534, 5.606 ,5 .677 ,5.748,5.816,5.882 ,5.947,6.008,6.067 ,6. 177, 
46.276y6« 366y 6 ♦ 446 y 6, 517 y 6.581 y 6. 638y 6.687 y 6 » 731 y 6 ,769 ,6.802,6.831 , 
56. 855 ,6.876, 6.894 ,6.950,6. 974,6.993,7.009, 7.036, 7.219,7.720,8.195, 
68,859 »9 ♦ 342 r 9 ♦ 748 ) 

GO r0<l,2,3»4,5> ,N 

1 CPG=ATKINT < TI ,CPP» T, 58, 6»NES» .01 >*4 . 184 
RETURN 

? CPO*ATKINT ' TI ,CPN»T »58y 6»NES, .01 >*4. 184 
RETURN 

3 CPQ«ATKINT <TI ,CPO,T y58y 6,NES, .01 >#4. 184 
RETURN 

4 CPU- ATK1NT (TI ,CPE,T ,58, 6,NES» .01 >*4 . 184 
RETURN 

5 TUP^II+.S 
TDN=TI-.5 
HUP-CPOH < TUP / 5 ) 

HDN=CPOH < TDN , 5 ) 

CPG- 1 HUP-HDN ) 

RETURN 

END 

FUNCTION CPOH(TI»N> 

DIMENSION T ( 58 ) y HP < 58 ) y HN ( 38 ) y HO ( 58 ) y HE ( 58 ) 

COMMON/PARA/PERCEN 

C CALCULATES THE ENTHALPY OF THE IDEAL GAS FOR H2 BY INTERPOLATION 

C DATA TAKEN FROM RP l?32y UNITS OF TABLES ARE CAL/MOL 

C UNITS OF OUTPUT ARE JOULE/MOL. THE INDEX N DETERMINES THE SPECIES 

C SPECIESyFOR N=1 yPARAHYDR0GEN»N=2 NORMAL, N=3 OR THO, N=4 EQUILIB 

C. N=5,S0ME ORTHO-PARA MIXTURE SPECIFIED BY COMMON /PARA/ » PERCENT 
C RANGE OF TEMP IS FROM 10 TO 5000K. 

DATA(T~ 

1 10. Oy 12 ^ f 14. 0- 16. 0» 10. Oy 20. Oy 25. Oy 30. Oy 35.0, 40.0f 45.0, 

2 50. Oy 55. Oy 60. Oy 65. Ot 70.0, 75.0y 80. 0» 85. Oy 90. Of 95.0fl00.0f 
3105. Of 1 tO.Oy 115. Of 120. Of 125.0f 130. Of 135. Of 140. Of 145. Of 150. Of 160. Of 
41/0.0f 180. 0» 190. Of 200. Of 210. Of 220. Of 230. Of 240. 0,250. Of 260. wf 270. Of 

5280.0. 290. 0. 300.0. 350.0.400.0.500.0.600. .700. f 1000. , 1500. ,2000. , 

63000. . 4000. .5000. > 

DATA ( HN- 

l 303.67, 313.60, 323.54, 333.48, 343.41, 353.35, 378.19, 403.03, 
l 42/. 86, 452. /I, 477.56, 502.43, 527.34, 552.32, 5/7.40, 602.62, 

3 628.02, 653.64, 6/9.51, 705.66, 732.13, 758.92, 786.06, 813.55, 

4 841.40, 869.61, 898.17, 927.08, 956,33, 985.91,1015.80,1045.99, 

51 to/, 22, 1169. 49, 1232. 71, 1296, 78, 1361. 60, 1427. 10, 1493. 20, 1559. 84, 
61626. 93, 1694. 44, 1762. 30. 1830. 48 ,1898.91 ,1967.57 ,2036.43, 2382. 74, 
/2/30 . 94 , 3429 ,46*41 29 .51, 4831 . 66 , 6966 . 23 , l 0697 .2,1 4679 . 2, 

823230 . 9 . 32345 .. 41 895 . ) 
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DA I A < HP* 
l J9.68, 


59.61, 69,55, 79.49, 89.47, 99.36, 124. 20 » 149.04, 


:•> j/a, 98, i.vB,/a» 223 , 61 , /48.5H, 2 / 3 . /t« :;w,n» 324 . 90 , 351 . 22 * 

3 3/8 * 22 , 406,01 * 434. /If 464. 38, 495.09, 526.84, 559.62, 593.41, 

4 628.14. 663.75, 700.14, 737.23, 774.92, H13.10, 851.69, 890.60, 

5 969.04* 1047.84 * 1 126.58* 1204.93 , 1 282.49, 1 359 .74,1 434.03, 1511. 56, 
A 1 586 , 36 * I 660 . 49 , 1 733 . 99 , 1 806 . *5 , I 8 79 .42,1951.47, 2023 , 15 , 83 7 7 . 8 3 , 
72729 .17, 3429 .24,41 29 . 48 , 4031 . 65 , 6966 .23, 1 069 7 . 2 , 1 4679 . 7 , 

823230 .9, 32345. ,41895. > 

UA ( A< HU- , . . 

1 388.33, 398.27, 408.20, 418.14, 428.07, 438.01, 462.85, 487.6V, 

2 512.53, 537.37, 562.21, 587.05, 61 1 .89, 636.73, 661.57, 686.42, 

3 711.29, 736.18, 761.11, 786.09, 811.14, 836.28, 861.54, 886.93. 

4 912.49, 938.23, 944.18, 990.37, 1016.80, 1043.51 , 1 070.50, 1 097. 78, 
51153.27,1210.04, 1268.09, 1327.39, 1387.91 ,1449.56, 1517.26,1575.93, 
61640. 46, 1705 .76, 1771 .74, 1838 . 32 , 1905.41 ,1972.94. 2040.86,2384.38, 
72731 .52, 3429 .53,4129 ♦ 52,4831 .66,6966.23, 10697 . 2 , 14679. 2, 

823230 .9,3 2345 . , 4 1895 ♦ ) 


2 512.53, 537.37, 562.21, 

3 711.29, 736.18, 761.11, 

4 912.49, 938.23, 944.18, 


DATA(HE* 

1 49.68, 59.62, 69.57, 79.56* 89.66, 

2 195.77, 234.90, 280.97, 326.24, 371.33, 

3 538.05, 575.93, 612.56, 648.13, 682.82, 

4 815.54, 847.76, 879.77, 911.63, 943.40, 


99,96, 127.50, 159.12, 
415.35, 457.84, 498.74. 
716.78, 750.14, 783.03, 
9/5, 11 ,1006.80, 1038.52, 


51102.11,1166.03, 1230 .39,1 295 .23,1 360 .58,1426. 43 ,1492.76,1 559 . 55 , 

61 626 , 75 , 1694 . 32 , 1 762 . 23 , 1830 ♦ 43 , 1898 . 88 , 1 967 . 55 , 2036 . 42 , 2382 . 74 , 
72730 . 94 , 3429 . 46 , 4129 . 51 , 4831 ♦ 66 , 6966 ♦ 23 , 10697 . 2 , 14679 . 2 , 

823230 . 9 , 32345 . , 41895 . ) 

80 10(1,2,3,4,5) ,N 

1 CPOH-ATK’im < TI,HP,T,58» 6, NES , . 01 ) *4 . 184 
RETURN 

2 CPOH- ATK I NT ( T I , HN , T , 58 , 6»NES, .01 >*4. 184 
RETURN 

3 CPGH*ATKINT< TI»H0»T,58» o»NES, .01 >*4. 184 
RETURN 

4 CPOH*ATKINT < TI ,HE»T ,58,6 , NES , . 01 ) *4 . 184 
RETURN 

5 PERCE “PERCEN /100. 

CPOH* < ATKIN 1 ( T 1 • HO , T , 58 , 6 , NES , . 0 1 ) * < 1 . -PERCE ) + 

1ATKINT < TI . HP, T, 58, 6, NES, .01 >*PERCL >*4.184 
RETURN 
ENLi 

FUNCTION CPOS< TI ,N> 

DIMENSION T ( 60 ) , SP < 60 ) , SN ( 60 ) , SO < 60 ) , BE < 60 ) 

CQMMON/PARA/PERCEN 

CALCULATES THE ENTROPY OF THE IDEAL UAS FOR H2 DY IN 1 E.KPOI. A I INC. 
DAKEN FRUM RP 1932, UNITS OF THE TABLES ARE CAL /MOl. DEO K. i 
UNITS OF OUTPUT ARE JOULES/MOl. DEU K. THE INDEX N DETERMINES THE 
SPECIES, FOR N"1 ,PARAHYDR0GEN,N=2 NORMAL, N*3 ORTHO, N*4 EQUILTB 
N ::: 5,S0ME ORTHO-PARA MIXTURE SPECIFIED BY COMMON /PARA/, PERCE NT 
RANUE OF TEMP IS FROM 10 TO 5000K. 

DAI A< T- , „ 

1 10,0, 12.0, 14.0, 16.0, 18.0, 20.0, 25.0, 30.0, 35.0, *0,0, 4...0, 

2 50.0, 55.0, 60.0, 65.0, 70.0, 75.0, 80.0, 85.0. 90.0, 95.0.100.0, 
3105. 0 * 1 10 .0, 115 .0, 120.0, 125.0, 130.0 , 1 35 0,140.0,145.0,150.0,160.0, 

4170.0. 180.0.190.0. 200 .0.210.0.220.0. 230 ♦ 0 , 240 . 0 , 250 . 0 * 260 . 0 * 2 70 . 0 . 

5280.0. 290.0.300.0*350.0,400.0, 450.0,500.0,550.0,600.0,700. , tOoO. . 

6 1 500 . . 2000 . . 3000 . . 4000 . . 5000 . ) 

111 1.214,12.1 20 * 1 2 . 886 ,13.54. ,14.135,14. 658 .15. 766 , 1 6 . 672 , 1 7 . 438 > 
118.102,18.688,19.214,19. 693 , 20 . 1 35 , 20 . 548 , 20 . 938 , 2 1 . 3 1 0 , 2 1 . 66V , 
n: >2 01 :> t 22 . 356 , 22 . 688 , 23 . 0 14 , 23 . 334 , 2 3 . 648 , 23 . 95 / , 24 . 260 , 24 . 55 7 * 

424 . 8 A 8 , 25 . 1. 32 ,25.41 0 , 25 .681 ,25.945, 26 . 45 1,26. 92° :• 2 * . 3 79 ,27. 802 • 

528 . 20 1 , 28 . 5 7 7 , 28 . 932 , 29 . 268 , 29 . 586 . 29 . 889 ,30.1 77,30. 452 * .30 . / 1 6 » 



630.969 

'31*212 

f 32. 305 f 

739.701 

f 42 ♦ 720 

f 45 ♦ 007 f 

DAI A <80* 


Ui 5*581 

.16*486 

f 1 7* 252 f 

222 ♦ 468 

t 23.053 

>23*576' 

326.215 

t 26 ♦ 500 

'26*771' 

428.375 

f 28. 575 

f28. 769f 

530. b08 

f 31 ♦ 109 

f 31 .401 f 

633.216 

f 33.446 

f 34 ♦ 505f 

741.884 

f 44 .903 

f 4 7 ♦ 190 f 

DAT A ( SH s ~ 


115.607 

f 16 ♦ 512 

f 1 7 ♦ 278 f 

222.494 

f23.079 

f 23 ♦ 603 f 

326.283 

f 26 .582 

f 26 ♦ 868 f 

428.611 

f 28 .832 

f 29 ♦ 047 f 

531.274 

f 31 ♦ 594 

f 31 ♦ 901 f 

733.772 

f 34. 005 

f35. 073 f 

742.455 

f 45.475 

f 47. 762 f 

DATACSE- 


111.215 

f!2. 120 

f 12.887 f 

219.294 

f 20. 331 

f 21 ♦ 285 f 

325.674 

f 26. 080 

f 26 ♦ 455 f 

428.510 

f 28 .750 

f 28 ♦ 980f 

531.269 

f 31 ♦ 591 

f 31 .899 f 

633.772 

f 34. 005 

f35.073f 

742.455 

f45.475 

f47.762f 

GO T0<1f2f3f4 

f5)fN 


33 . 244 , 34 ♦ 069 . 34 . BOA , 35 ,473 , 36 .082, 37 . 1 A5 , 
51 , 221 , 53 . 839 * 55 . 9A9 ) 

1/.916, 18,501, J 9.024 ,20. 133, 21. 038, 21,804, 
24. 050, 24. 482 ,74. 880, 25. 248, 25. 591, 25.912, 
27 . 029 , 27 , 275 , 2-> . 5 1 2 , 27 , 739 , 27 . 958 » 28 . 170 , 
28. 958, 29. 143, 29. 502, 29. 84A, 30. 177, 30. 498, 
31 , A84 » 3i ♦ 958 , 32 ♦ 225,32.484 , 32 . 735,32 .97^, 
35 . 432 , 3A ♦ 253 , 3A . 989 , 37 . ASA , 38 . 2A5 , 39 . 348 , 
51.221,53. 839 , 55 . 9A9 ) 

17.941,18.527, 19.050,20. 159,21. 0A4 ,21 .830, 
24. 078, 24. 513, 24. 914, 25. 288, 25. A38, 25. 9A9, 
27. 143, 27. 407, 27. AA3, 27. 91 1,28. 151, 28. 384, 
29 ♦ 256, 29. 4A1, 29 .85A, 30. 234, 30. 595, 30. 942, 
32.197,32.483,32.758,33.025,33.282,33.531, 
36.003,36.825,37.561,38.228,38.836,39.920, 
51.221,53.839,55.969) 

13.554,14.149,14.692,15.918,17.069,18.196, 
22.145,22.911,23.592,24.198,24.740,25.229, 
26.804,27.129,27.435,27.724,27.999,28.260, 
29.203,29.418,29.828,30.216,30.584,30.934, 
32.196,32. 482 , 32 ♦ 758 ,33.024,33.282,33.531, 
36 < 003 , 36 . 825 ,37.561,38. 228 , 38 . 836 , 39 . 920 , 
51.221,53.839,55.969) 


6,NES,.01>#4.184 
6,NES, .01)*4.184 


6,NES».01>#4.184 


1 CP0S=ATKINT<TI»SP,T,60, 6»NES» ,01 )*4. 184 
RETURN 

2 CPOS=ATK INr<TI,SN,T,AO, 6,NES, .01 >*4. 184 
RE TURN 

3 CPOS~ATKINT <TI,S0,T*60> 6,NES, .01 )*4, 184 
RETURN 

4 CPOS=ATKINT <TI,SE,T,AO, 6,NES, .01 >#4. 184 
RETURN 

5 PERCE =PERCEN /100. 

CPOS= ! < ATKINT <TI ,SO,T ,60,A,NES, .01 )#< 1 ♦ -PERCE ) + 

1ATKTNT < T I , SP , T , 60 , 6 , NES , . 01 ) #PtRCE ) *4 . 184 
THE EXPRESSION FOR THE ENTROPY OF EQUILIBRIUM H2 IS INCORRECT, 

THE ENTROPY OF MIXING MUXT BE ADDED TO MAKE IT COMPLETE 

RE URN 

END 

FUNCTION ATKINT ( X , YMAT , XMAT , NELM IS , NMAX » NESSY , ACRCY ) 

THIS PROGRAM HAS BEEN CHANGED SO THAT THE OSCILLATING NATURE OF 
THE MATRIX TO BE INTERPOLATED EXISTS ONLY AT THE UPPER END OF THE 

table; 

THIS ROUTINE WILL TAKE INPUT MATRICES OF UP TO 999 ELEMENTS EACH, 
ARRANGED SO THAT THE X MATRIX(XMAT) IS IN EITHER ASCENDING OR 
DESCENDING ORDER, SELECT NMAX OF THESE P01N IS, CHOSEN SO THAT 
SUCESSIVE X VALUES OSCILAIE ABOUT THE VALUE OF THE ARGUMENT X 
UNLESS THE ENDS OF THE XMATRIX INTERFERE (IN THIS CASE THE 
OSCIL.ATQRY NATURE IS LOST BUT THE PROGRAM WTLL STILL. PERFORM AN 
INTERPOLATION) , INTERPOLATE ON THESE NMAX PAIRS OF DATA BY 
AN 0SC1I AUNG VARIABLE POINT AITKEN INTERPOLATION ALGORITHM 
EI THER UNTIL THE PERCENTAGE CHANGE IN THE INTERPOL ANT JS l F SS 

I HAN THE ACRCY ARGUMENT < THE ARGUMENT NESSY INDICATES THE 
NUMBER OF I HE POINT JUST BEFORE l HE LAST ONE CHECKED) OR UNI 11 
THE NMAX POINTS ARE All USED. IT IS SUGGESTED THAT NMAX 

BE LESS THAN 10, AND OF COURSE. LESS THAN NEIMTS. NILMIS 
INDICATES THE NUMBER OF ELEMENTS IN XMAT OR YMAT. 

II NESSY IS ZERO II INDICATES THAI THE INTERPOLATION REQUIRE MEN I 
HAS NOT BEEN SATISFIED. IE NISSY IS 1 It MEANS THAI I HE VAl (IE 01 
X LIES OUT SIDE THE RANGE OF XMAT. 


O 


anno no 


DIMENSION YMAT < 999 ) , XMAT (999) r A< 21 , 20 ) 

100 FORMAT < 42H1NTERPQLATIQN REQUIREMENT NOT SATISF 1ED( X«,E16* 8 , 1H >/33H 
1LAST 2 APPROXIMATIONS OF Y ARE < Y» yE 16,8, 1H, ,E16,8,1H) > 

200 FORMAT <55HTHIS REPRESENTS AN EXTRAPOLATION OF THE XMAT MATRIX<X*, 
1E16*8,1H) /33HN0 CALCULATION HAS BEEN PERFORMED) 

300 FORMAT (24HNELMTS IS LESS THAN NMAX) 

400 FORMAT (22HNMAX IS LARGER THAN 20) 

IF <NMAX~20) 71 ,71,69 
69 PRINT 400 
ATKINT»0*0 
RETURN 

71 IF < NMAX- NELM IS ) 75 , 75 , 73 
73 PRINT 300 
ATKINT *0 ♦ 0 
RETURN 
75 CONTINUE 

FIRST TWO SUCCESSIVE VALUES OF THE XMATRIX THAT STRADDLE THE 
VALUE X WILL BE SOUGHT 
JJl-NELMTS-l 
DO 20 1*1, JJ1 
DIF1*X~XMAT ( I ) 

DIF2«XMAT<I+1)-X 
IF<DIF1>16,15,16 

15 ATKINT«YMAT ( I ) 

NESSY =NMAX 
RETURN 

16 IF<DIF2)18,17,18 

17 ATKINT*YMAT < If 1 ) 

NESSY =NMAX 
RETURN 

18 RATION'D IF 1/ DIF2 
IF (RATIO) 20, 20, 19 

19 IMID*I 
GO TO 32 

20 CONTINUE 

AT THIS POINT ONE COULD PRINT THE FOLLOWING STATEMENT 
WRITE OUTPUT TAPE 6,200,X 
NESSY*! 

ATKINT*0*0 
RETURN 

32 CONTINUE 

NOTE THAT RATIO IS POSITIVE IF THE TWO POINTS STRADDLE X 
REGARDLESS WHICH IS LARGER 
JJJ»IHID 
JUP»XKID 
JDN-IMID 

IF< JJJ+NMAX-NELMTS+1 )98,*8, 102 
98 DO 201 J*1 ,NMAX 
JJJ*IMID+J-1 
A < 1 , J ) *XMAT < J J J ) 

201 A(2, J ) *YMAT < JJJ) 

GO TO 203 

102 DO 41 J*1 ,NMAX 

JJ* J/2 
JOE- -J-2KJJ 

C JOE IS 0 IF J IS EVEN AND 1 IF J IS ODD 
IF< J~1 )33,40,33 

33 IF < JDN-1 >34,36,34 

34 IF < JUP-NELMTS )35, 37? 35 

35 IF < JOE >37,36,37 

36 JUP® JUP+ 1 
JJJ* JUP 
GO TO 40 


fi 1 


3/ JDN»JUN~1 
JJ J^JUN 
OU TO 40 

40 A<1* J)~XMAf <JJJ> 

A<2r J)^YilAT(JJJ) 

41 CONTINUE 
203 NNN-*NMAX+1 

DO 6 NNN 
L *““ J ~ l 

I. *0 5 K®L f NMAX 

C J IS THE COL f MN NUMBER 

C K IS THE ROW NUMBER 

0A< JfK)*<A< J~ lfK)-A<J-lf J~2> >#<X~A(1 f J-2) > /< A< 1 fK > ~A < 1 f J~2>> 

1 +A^J~1fJ-2> 

IT <K-L>3f2f3 

2 IF<ABS< <A<JfL)-A< J 1 f L-l ) >/A< JfL) >-ACRCY/100*0) 7 f 7f 3 

3 CONTINUE 

5 CONTINUE 

6 CONTINUE 
NLSSY-0 

C AT THIS POINT ONE COULD PRINT OUT THE FOLLOWING STATEMENT ♦ 

C WRITE OUTPUT TAPE 6t 100 fX f A<NNNfNMAX ) f A<NNN-1 fNMAX-1 ) 

ATKINT»A(NNNvNMAX) 

RETURN 

7 NESSYwJ-l 
ATKINT~A< Jt L ) 

RETURN 

END 

SUBRlJU I J. NE !.. X h I T S < P ? 7 t It. ) 

DIMENSION 0(32) ? VP ( V ) 

UJMMON/ DA TA/G » R f GAMMA f UP f D T P f PCC * FTP t ICC f TIP t TUL f I'LL f PUP 

it < p*ur ♦ pup)uo 10 io 

J. F < ( « U T * fUL ♦ OR ♦ T ♦ L T ♦ TL.L ) GO T O 12 
IL-1 

RE I URN 

10 PR IN I Uf PUP 

It FORMAT <# THE INPUT PRESSURE IS OUT OF THE RANGE OF THIS EQUATION * 

l/* THE RANGE FUR THIS EQUATIONS IS FROM 0 TO #fF6*0f# BAR# ) 

1L*=V 
RE 1 URN 

J2 PR IN f UfILL.fTUL 

13 FURMATl# THE INPUT TEMPERATURE IS OUT OF RANGE# 

1/# T HE RANGE FOR THIS EQUATION IS FROM #fF6*2f# TO #fF6*0f# K#) 

XL 0 

RETURN 

END 

SUBROUTINE DATA NE 
PR IN I 10*j 

ION IORMAK* THE TEMPERA TURK RANGE OF NEON IS 25 TO 300 KELVIN# 

F* WITH PRESSURES TO 200 BAR# 

2/# THU MAY ENTER THIS ROUTINE WITH ANY TWO OF THE VARIABLES# 

I./# PfD UR If AND 0 I- OR THE THIRD f FOR SATURATION ENTER WITH* 

4/# I I TITER P OR If AND 0 FOR THE OTHER TWO fIF YOU WANT # 

G/* 10 CHANGE FLUIDS ENTER 0 FOR A U THREE INPUT VARIABLES# ) 

44*4 

P( -JA.S4 
TIP 20* 

I UL -200 * 

I HI 
TIL- 

I M 30.103 
PRINT 104 
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p 



k, 
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V* 

‘ 
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! 

8 i 

•' \ 
i 



|.|^ FURMAH* no YOU WAN! F NOT NEE KING I IN MS UK ME (RIO UNI I H 
I / ♦ LN I Eh: A 0 HIE ENGINEER I NO UNITS OR A I I UK ME! RIG* > 

I mI I. li-HA N 
hi AH 102? ID 

U>:\ I OKMA Hit) 

1/ IHUEKEO^hU lU IB 
I ‘KIN I .103 

lv % \ HJkhnK* I: NT I" R PRESSURE IN HAR? DENSITY IN MOL F S/I I TEE ? AND f HIRER 
l A I UK!.. IN KELVl NS*> 

34 Kl. All * ? P ? D ? I 
P*-P/l ♦ 01 32t> 

Gil IU IV 
10 EKiN I 104 

106 FUhflAH* ENTER PREbSURF IN LJtf/SQ IN? DfcNSi I’Y IN I HUH) Ff AND* 

1/* ft.MEfc.RA FURL IN DEGREES F* ) 

Rfc.AU *?P?Br! 

P^IVI4*6V5<M9 

U' : J*16«0J8463Sl/EM 

I I * r » lu*o*o>go ro iv 

1*0 -152 ♦ ‘/.I ♦8+2/3 * 1 5 
IV LON'f INUfc. 

IMP*LE'*0*0)G0 TO 14 
IF<D*LE *0*0>U0 TO 13 
IF < T ♦ LE ♦ 0 ♦ 0 >G0 TO 12 
BO TO 21 

12 t F < P ♦ LE ♦ 0 ♦ 0 ♦ OR * D ♦ l. E ♦ 0 ♦ 0 ) GO TO 15 
n=FNDTNE<RfD) 

I F < P ♦ OT ♦ PUL ♦ OR ♦ T ♦ G r ♦ r IJL ♦ OR * T * LT * TLL ) GO TO 25 
GO TO 16 

13 IFO*LE*0*0*0R*P*LE*0*0)G0 TO 15 

IF < P ♦ G !’ ♦ PUL ♦ OR ♦ T ♦ GT ♦ FUL ♦ OR ♦ T ♦ LT ♦ ILL > GO TO 25 
131 D-FNllDNECPf T) 

GO TO 16 

14 IF (I* I F ♦ 0 ♦ 0 ♦ OK ♦ I ♦ LE ♦ 0 ♦ 0 ) GO TO 15 
P=FNDPNE<D»T> 

1 F < P ♦ G F ♦ PUL * OR * T * GT * T Ml ♦ OR ♦ T ♦ LT ♦ TLL ) GO TO 25 
16 h*enihnl:<p?D'T> 

S«ENTKNE <D?T ) 

PP«P*1 *01325 
I F < I U ♦ EU ♦ 0 > GO TO 23 
PRINT 1C0?PP?D? T ?H?S 

100 FORMAT < * PRESSURE** fF7*2f * BAR ? DENSITY**? F7 *4? * MOL ES/L ITER ? IE* 

1 *MPER A TURE * * t F6 ♦ 2 ? * KELVIN* 

2/* ENTHALPY-* ?FV* 2?* JOULES/MOL ? ENTROPY** ? FA *2?* JOULES/MOL-K* ) 

GO TO 1/ 

23 PQ~P*14 ♦ 695V4V 
TlMT*l *8-459*67 
DO*D*£ M/ 16 ♦ 018463/1 
H0*H/(2*324445*EM) 

HU-3/ (4*1 8400 1*EM ) 

PRINT 101 ? PO? DO? 1 0? HO ? SO 

101 FORMA! <* PRESSURE-*? FI 0*3?* LB /SO IN? DENS I TY=* ? F10 * 4 ? * LB/CU FI* 
1/* I EMPfc NATURE** ? F8 ♦ 2 ? * F? ENTHALPY** ?F 9 *2?* BTIJ/LB* 

7/* EN fRQPY**? FA* 2 ? * BFU/LB-F*) 

GO 10 18 

15 JF <P*LI ♦0*0*AND* l*LF*0*0)G0 TO 21 
IF <P*nr*0*0>Gfl TO 20 

1F< HOT * TC*OR* I.Lf* FTP >00 TO 22 
I ' -VF'NNE ( r ) 

tOV FORMAT (* DO YOU WAN! SAIURAIE.D IIOUJ D OR VAPOR* 
l/* ENIER A i FOR VAPOR OR A 0 FUR L IUU I D* > 

PR IN I 1W 
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READ *.IP 

if!ip,gt,o>p»p-,oooi 

IF!IP,EG,0>P«P+,<)0001 
00 TO 131 

20 IF!P.GT.PC)GO TO 22 
T»VPTENE!P) 

PRINT 109 
READ id » IP 

IP< IP.GT.0)P»P~.0001 
IF ( IP»EQ.0)P»P+.00001 
00 TO 131 
22 PRINT 107 

107 FORMAT <* YOUR INPUT PRESSURE OR TEMPERATURE IS OUT OF RANGE* 

1/* OF THE SATURATION CURVE* Tl>44.4. PC®26,89. T TRIPLETS.* 

2/* TRY AGAIN*) 

GO TO 17 

21 RETURN 

25 PRINT 108 

108 FORMAT (* YOUR INPUT TEMPERATURE OR PRESSURE IS OUT OF RANGE* 

1/* TEMP MUST BE BETWEEN 25 AND 300K AND PRESSURE BETWEEN* 

2/* 0 AND 200BAR. TRY AGAIN*) 

GO TO 17 
END 

SUBROUTINE DATA HE 
PRINT 105 

105 FORMAT!* THE TEMPERATURE RANGE OF HELIUM IS 2.177 TO 1500 KELVIN* 
17* WITH PRESSURES TO 1000 BAR* 

27* YOU MAY ENTER THIS ROUTINE WITH ANY TWO OF THE VARIABLES* 

3/* P.D OR T. AND 0 FOR THE THIRD. FOR SATURATION ENTER WITH* 

4/* EITHER P OR T. AND A 0 FOR THE OTHER TWO . IF YOU WANT * 

57* TO CHANGE FLUIDS ENTER A 0 FOR ALL THREE INPUT VARIABLES*) 
TD=5.201 
PC=2 . 245 
TTP=2 .17 
PUL-1000. 

TUI.® 1500. 

TLL =2 . 177 
EM=4,0026 
PRINT 104 

104 FORMAT!* DO YOU WANT ENGINEERING UNITS OR METRIC UNITS ?* 

17* ENTER A 0 FOR ENGINEERING UNITS OR A 1 FOR METRIC*) 

READ 102.11' 

102 FORMAT! ID 

17 IF ! IU ♦ EG . 0 )G0 TO 18 
PRINT 10'. 

103 FORMAT!* ENTER PRESSURE IN BAR. DENSITY IN MOLES/LITER. AND TEMPER 
lHTURE IN KELVINS*) 

24 READ *.P.D.T 
P*P/1 .01325 
GO TO 19 

18 PRINT 106 

106 FORMAT!* ENTER PRESSURE IN LB7SG IN. DENSITY IN LB/CU FT AND* 

1/* TEMPERATURE IN DEGREES F*> 

READ * > P » D . T 
P*P7 14. 695949 
0“D* 16. 01 846371 /EM 
IF !T. EG, 0.0)00 TO 19 
T- ! T-32. >/ 1.8+273. 15 
'9 CONTINUE 

J F ! P , LE . 0 ♦ 0 ) GO TO 14 
IF !D«LE.0,0)G0 TO 13 
IF !T.LE .O.OIGO IU 12 
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GO TO 21 

12 IF < P . LE . 0 . 0 . OR . D . LE • 0 . 0 > GO TO 15 
T»FNDTHE<P»D> 

IF < P . GT . PUL. * OR . T . GT . Till... . OR . r . LT . TI L. ) GO TO 25 
GO TO 14 

13 IFtT.LE. 0*0. OR. P.L.E. 0.0)00 TO 15 

IF < P ♦ GT ♦ PUL ♦ OR ♦ T • BT ♦ TUL « OR » T ♦ L T ♦ TLL. ) GO TO 25 
131 D»F NDDME ( P » T ) 

GO TO 16 

14 IF < D . LE . 0 . 0 . OR . T » LE » 0 » 0 ) GO TO 15 
P®FNDPHE<D»T> 

IF (P.GT .PUL.QR.T *0T * TUL.QR.T .LT » TLL)GO TO 25 
16 H«ENTHHE<Di>T> 

S»ENTRHE(D»T> 

W“SGUNHE<D»T> 

CPP®CPME<D»T > 

CVV«CVPIE<DrT> 

TH®THERHE(D»T>*100. 

V®VISCHE<D»T> 

PP®P*1 .01325 

IF< XU.EO.O)GO TO 23 

PRINT 100fPF*»D>TfH»S» CPP . CUV » W » TH » V 

100 FORMAT!* PRESSURE®* . F8 . 2 . # BAR. DENSITY=*»F7.4»* MOLES/LITER* TE* 
1*MPERATURE®*»F7.2»* KELVIN* 

2/* ENTHALPY®*.F9.2»* JQULES/MQL. ENTR0PY=*.F6.2>* JOULES/MOL-K* 
3/* CP«*»F7.2»* JOULES/MOL-K f CV=*»F7.2»* JOULES/MOL-K* 

4/* SPEED OF SOUNEi®* » F8 * 2 » * METERS/SEC* 

5/* THERMAL CONDUCTIVITY®* rF7.2»* M WATS/M-K* 

6/* V r SC0SITY®*rF7.2»* MICRO G/CM-SIO 
GO TO 17 

23 PO®P* 14. 695949 
T0=T*1. 8-459. 67 
DQ=-D*EM/16 .01846371 
H0-H/(2.3244-. j#EM> 

S0=S/(4.184001*EM) 

CPP0=CPP/<4. 184001*EM) 

CVVO=CVV/ < 4 ♦ 184001 *EM ) 

W0®W*3. 280840 
THO®TH*. 0005781 76 
VO®V*. 0067196897 

PRINT 101»P0»D0»TQ»H0»S0»CFP0»CVV0»W0>TH0»V0 

101 FORMAT < * PRESSURE®* >F10. 3.* l.B/SQ IN. DENSITy?*.F10.4.* LB/CU FT* 
1/* TEMPERATURE®*. F8. 2.* F. ENTHALPY®*. F9.2. * BTU/LB* 

2/* ENTROPY- *»F6.2»* BTU/LB-F* 

3/* CP=*.F7.3»* AND CV=*.F7.3»* BTU/LB-F* 

4/* SPEED OF SOUND®* . F8 ♦ 1 * * FT/SEC* 

5/* THERMAL CONDUCT I VI IY®*»F7.2»* BTU/F T-HR-F* 

6/* VISC0SITY®*.F7 .2** LB/FT-S X Ef/*> 

GO TO 18 

15 IMP. LE. 0.0. AND. T.LE.O.O)GO TO 21 
1F(P.GT. 0.0)00 TO 20 

IF (T .GT .TC.OR.T ,I.T.TTF‘)GO TO 22 
P®VPNHE ( T ) 

PRINT 109 
READ *.IP 

IF ( IP ♦ GT » 0 ) P®P- » Ov. j 1 
IF ( IP.GE. 1) P®P F. 00001 

109 FORMAT < ♦ DO YOU WANT SATURATED LIQUID OR VAPOR* 

1/* ENTER A 1 FOR VAPOR OR A 0 FOR LIQUID*) 

GU TO 131 

20 IF < P.GT, PC) GO TO 2? 

T -VP' TEHE < F* ) 
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A 


o 


PRINT 1W 
RE All **IP 

IF ( IP* (if *0>P*>P-.0001 
IF < IP. of:, Dp *p f . 0000 1 
00 TO 111 
22 PRINT 10/ 

10/ FORMA f t * YOUR INPUT PRPSSURT 
1/# OF 1 HE SATURATION CURVE* 

2/* TRY AOAINKO 
00 10 17 
21. RETURN 
2b PRINT 10B 

10W FORMAT (* YOUR INF'UI TEMF'FRA TURF UR PRI-BBURF 
l/* TEMP MUST BE BETWFJ N 2.17/ AND lbOOK ANTI 


UR TEMPERA HIRE IB OUT 
TOb.JOl ♦ PO 2 , 2/5 * I 


i IT RANliT * 

R J PI T . 1 //* 


18 OUT OF RANGE* 
PRESSURE BE T WEEN* 


PROGRAM ASKS FOR A FLUID SELECTION* ENTER THE AP 
WILI TERMINATE THE PROGRAM*) 


2/* 0 AND tOOOBAR* TRY AGAIN*) 

GO TO 1/ 

END 

SUBROUTINE INFO 
PRINT 2 

2 FORMAT <* WHEN THE 
1PR0PRIATE NUMBER* 

2/* AN INAPPROPRIATE NUMBER 
PRINT 1 

\rSp?SIi.^t- WH cMrIc E A [! R0GRAM ASKS F0R A ASSURE* DENSITY* AND*/* I 

ANY 2 ° F THE THR|;:E AND A 0 FOR THE THIRD.*/* TH 

?L ^TF A. ? ™ Jr F ’ r, : T : ANn 0NE 0F rHE three must be 0. * 

3/# IF ALL THRLL ARE 0 THE PROGRAM ASKS FOR A NEW FLUID* ) 

PRINT 110 

UO i/2 R UnI^r^r^ A ^ CRESTED IN A DEFINITION OF THE VARIOUS* 

1/* MODIE ILAT10NS OF HYDROGEN ENTER A 1* IF NOT ENTER A 0*) 

READ **IC 

IFdC.NE.I ) GO TO 111 
CALL H2 INFO 
111 PRINT 1.12 

U2 1/Te^eS A F |. Y ?f NOT eSt T ErToJ) IN ™ F S ° URCES ° F THESE PR ° GRAMS * 
READ **IS 

IFdS.EQ.l )CAL.L SOURCE 

RE I URN 

END 

SUBROUTINE H2 INFO 
PRINT 100 

iOO FORMAT!* THE H™KOGEN MOLECULE IS MADE UP OF TWO HYDROGEN ATOMS* 

Jptw nr ful ath2o S ^ LE RELATIVb ORIENTATIONS OF THE NUCLEAR* 
rL HE ^ T0MS IN A MOLECULE OF HYDROGEN* I*E* THE SPINS* 

IN THE SAME DIRECTION (ORTHO) OR IN OPPOSITE DIRECTIONS* 
(PARA)* THE RELATIVE AMOUNTS OF ORTHO AND F*ARA HYDROGEN IN* 

r.o, A SAMFLE 0F EQUILIBRIUM HYDROGEN IS TEMPERATURE DEPENDANT*) 

F*RIN1 101 

* 01 ^i! R ^ Al ** AT r 00M TEMPERATURE EQUILIBRIUM HYDROGEN IS 75 PERCENT* 

}/* !? rt !° AND 25 PERCENT PARA AND IS CAl.LED NORMAL HYDROGEN* 

'‘■^rnlx.*' 0 KELV1N EQUILIBRIUM HYDROGEN IS 99,8 PERCENT PARA.*) 
hfc I UhN 

end 

SUBROUTINE SOURCE- 
PRINT 100 

100 FORMAT!* THIS PROGRAM WAS WRITTEN BY R.D. MCCARTY AT THE NATIONAL* 
1/* BUREAU OF STANDARDS IN BOULDER COLORADO UNDER A CONTRACT WITH* 
THE JOHNSON SPACE CENTER IN HOUSTON TEXAS? THE NASA CONTRACT* 

JSC 18 W,LIE * SCtlrr - THE OOCUWEMT*!™ !!f ™e5 * 

FRGGRA H 18 IN PREPARATION AND IS PLANNED AS A JOINT* 

NASA-NDS PUBLICATION* ) 


2 /* 

3/* 

4/* 

5/* 


2 /* 

3/* 

4/* 

U/* 


M: fURN 
t Nil 


MINI.! r l(.IN I NDDHE < P I , I J ) 

f! ’{(UVt'i I HI' HU. I lift IUUAIIUN Oh SIAIE PUR DFN8 1 1 Y IN MOl./llIER 

' p,;<E - sa,,,l?E lN atmospheres and a iemperaturk in degrees k 

p I • I 

11 t T.LT .5.2)00 III A 
PM-' 3 .1. 00 1 , 

CM I .I T. 1*3. >PM~PMtl MI) 

IFtPM.LT .P>GU I'll .10 
II U ,GT. 100. >011 TO 1 
I *t’ '■:> . 2449 + ( 1-5 , 201 4 > #1 . 76 
IE < P ♦ LT • PC > GO TO J 
PM~200.+< T-5,2>*t2*31 

D=*17 , ,199+ < < PM-PC ) / ( PM -POE t ♦ ) ) #2 , 33* 1 7 . 39V 
GO TO 7 
2 0‘*.0001 

IE ( MET. 4.2)00 TO 7 
l VB !S VIRB< T) 

RT=0 . 0820558#T 

Pl-RT/P 

D=1 . / (Pl+VB ) 

GO TO 7 

6 IFIP.LT .VPNHEtT) >60 TO 2 
DS=46. 18+ < T -2. >*4.02 
DL-DSATL ( T ) *1 000 . /4 . 0026 
DEL-DS-DL 

PM*PMELT f T ) 

I FtP.GT.PM nlO TO 30 
D=DL+DEL*P/PM 

7 DO 10 1=1, SO 
P2=FNDPHE <D»T) 

IF<ABS<P-P2)-1,E-7*P>20*20»8 

8 DP==DPDDHE < D * T ) 

CORR* ( P2-P ) /DP 

IF<ABS<C0RR>-1.E-7*D>20*20*10 
10 D=D-C0RR 
FNDDHEaO 
RETURN 


C 

C 

c 

c 


20 FNDDHE=D 
RETURN 
30 ENDDHE =0 
RETURN 
END 

FUND T I CN FNDPHE ( DD * T T > 

THE. MAIN EQUATION OF STATE ROUTII'I FOR HELIUM I 
THE INPUT IS DENSITY IN MOLES/LI ILK AND TEMPERATURE 
IN DEGREES K * THE OUTPUT IS PRESSURE* OP/DDf DP/D T* ENTHALPY* 
ENTROPY AND CV IN ATMOSPHERES AND JOULES 
DIMENSION A(26) *B<27*4) 

COMMON/ IDtX/ ID 


DA I A ( ID- 0) 

DAIAt < B < 1 ) * I -82 * 1 08 > 

1 5 . 988310 109E-9* -4 . 96530521 8 7E- 7 

13E-5* -1 . 4830691828E-4 * 3 .05961 74335E-4 * 
22E--4* 1 .5527899/12E-5* -3. 61 10403503E-5 * 
3/E-5» 1 .938145109E-5.-4, 1496408961?- 4t -5 
4E-3 >-6 . 83838U9924E-2* --2. 1 3B2474225E--2 * 
08E - 2 * 2 . 50/S/5338E-3 * 7 . 9041 608815E-2 * - 1 
6E - 7,1, 73364 10358E-6 1 -2 . 5454 1 8 7855E- 6 * 
DAIAt ( B ( I ) * 1-55*81 > --1 .4802195348F -8*4 


*-3 .81 16033499E- 6» 3 .681 1 6713 
-3 . 39081 90224E-4 * 1 , 962408024 
-1..0R39/88073E -5*4. 972810121 
. 7465772899E-4 * -4 . 3470945634 
2 ♦ 710A954908E-2*-l ,262796778 
4024724318E-4*- 2.8278987249 
0.0005) 

.1721 791 1 19E 7 * -2 . 3326553271 
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*E~7 » 4 . 00551 10B8E-7 » 1 * 09005A7964E- 5.-5. 0060952775K-5# 1 . t .ft 2765043 
2E-4#-l . 253984320 7E-4# 1 .9661380688E A# 1 , 71779376461?' -4 #2 .30510005 <3 
3E~ 4 » -9 . 65647391E-4 » -3 . 6027 735292E-5 » 1 ♦ 607994655SE-3 » -2 . 7441 763615 
4E-2# , 14739506957 » - ♦ 43559344838* 1.3447956078* -1.7040375125# .9026247 
54040# 5.68756441 1 1C-3 * -1 . 44 J814AA25E- 1 *3.3760874851 f -3# 1 .0754201 21 8 
6E--6 # -4 . 5264622308E-S # 3 . 8597388864E-5 » - . 0005 > 

OA T A ( ( B < I ) . I *28 *54>«-4» 228 745462 AE -8#4. 45293544 1 3E- 7 » - 1 .0246 1 50954 
IE-5 # 8 ♦ 5284608956E-5 #-2.51 63069255E-4 » 3 ♦ 287 77092H5E -4 » • 1 .060195758 
2E-4#~1 . 0687738074E-4 # -3 . 2120950632E-5 # 1 .4159010971 -4*1 .4725630701 
4E-3 # -2 . 61035 4941E-3 *2.0461501 117E-5# 1 . 2746996288E-3 #-2.0272929583 
5E-2 » 7. 4648036415E-2#-. 1721 7966521 . .51053439738* -.401 78202497# 

6 ♦ 26829864632 # 7 . 906601 204E-3 # - 8 . 9393485656E-2 # - . 1 5076580053 » 

72. 6882494327E- A » ~3 .379431 6832> fc:~5» -2. 449595 119SE-5# -.0005) 
DATAttBtI)#I=l#27>= 

A -1 ♦ 50968626 19E-7 # A . 4A40898904E-7 » 4 . 1 36235736 7E-5# 

1-3. 7910190353E-4 # 1 . 3806454049E-3# - 2. 50854 12058F-3# 2 . 36975A0398E-3 
2# -9 . 572646 1066E-4 # 3 ♦ 7405931828E-5#-6. 4103220333E-4 # 

3 1 ♦ 85793661 77E-3 » 7 ♦ 4007986606E-4 » 1 . 4792568 1 48E-4 # ~3 . 2531355477F-3 » 
41 .9518739286E-2# -♦ 1057161 7135 #. 33164944449# - .51 130022535# 

53 . 9940004906E- 1 # - ♦ 15555244471 » 4 . 906264031 E-3 # -2 . 61 48004377E-2 # 
63.4221685545E-2# 5 . 4159662622E-6 »-l .0687806777E-5#- 8.9484651SA9E-6# 
7-. 0025) 

KP=1 

GO TO 10 
ENTRY DPDTHE 
KP=2 

GO TO 10 
ENTRY DPDD HE 
KP=3 

BO TO 10 
ENTRY ENTRHE 
KP=4 

GO TO 10 
ENTRY ENTHHE 
KP=5 

GO TO 10 
ENTRY CVHE 
KP=6 

10 Kl-0 
KH=1 
K~1 
KK=1 
K4=l 

IF t ID«NE.O)GO TO 20 
IFtTT.GE.15. >GQ TO 20 
IFtTT.GT. iO)GO TO 30 
IF < CD. GT. 17.3987)60 TO 40 

11 1*1 
T=TT 
0=00 

8 GO TO < 9 # 1 00 » 200 » 300 # 400 » 500 ) NP 

9 02=0*0 
03*02*0 
04=03*0 
05=04*0 
GAMMA=B< 27 # 1 ) 

EX=EXP(D2*GAMHA) 

EXD3=EX*D3 
EXD5- EX*05 

M ; -l 

N=t 

At N> 05*0 * N=N+t 
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A(H;-A(N-n/T $ N-N+l 
DO 2 1-1,6 
FI -I 

A < N ) - nS#T## < . 75- f 1/4. ) 

2 N«N+1 

DO 3 1-1,4 
M = I 

A < N > =-D4# I **< 1 , 5-FI ) 

3 N m N+1 

DO 1 1*1,8 
F I « I 

A<N>«D3#I*#<t,5~FI/2. ) 

N-N+l 

1 CONTINUE 
DO 4 1=1,3 
FI® I 

A<N)=EXD3#T#*<1 ,~F!> 

4 N=N+1 

DO 5 1=1,3 
F1*I 

A<N>=EXD5*T*#<1,-FT> 

5 N=N+1 
N=N-1 
I=M 

7 P=0 . 

DO 15 J=1,N 9 

15 P=F*+ B (J,I)#A(J) 

P»P+ . 08205S8*D*T* < 1 . +VIRB ( T ) *D ) 
1F(KH.LT.1)G0 TO 413 
GO TO< 50,50,30,40)K 

100 D2=D*D 
D3=D##3 
D4=D3*D 
D5=D4*D 
DA=D5*D 
T2=T*T 
T3=T2*T 
T4=T##4 
M = I 

GAMMA=B<27,M> 

EX=EXP<D2#GAMMA> 

N=1 

R=. 0820558 
A<N)=0,0 $ N=N+1 

A<N>=<-1. >#D6/T2 * N=N+1 

DO 102 1=1, A 
FI* I 

A<N)-=D5*T#*< .75-FI/4.-1. )*( .75-FI/4. ) 

102 N=N+1 

DO 103 1=1,4 
F I = I 

A < N ) =D4#T ** < 1 » 5-F 1/1. -1. )#<1 .5 -FI) 

103 N=N+1 

DO 101 1=1,8 
F 1 =1 

A < N ) =D3#T**< 1 »5-F 1/2, -1, )#<1, 5~F 1/2, ) 

101 N=N+1 

DO 104 1=1,3 
FI* I 

A ( N ) -=EX#D3# I ## < t , -F I 1 , )*( 1 ,-Fl ) 

104 N=N+1 

DO 105 1=1,3 


fid 


11 I 

A < N ) -fc XiMiUti I ## < 1 , I t - 1 , » * < I . ||) 
loti N-NU 
N» N-t 
P 0 

MU t Hi ,1-1 f N 
11!i P P+A< i>$U< M1» 

PPHmi*< I . I 0 1 f < H ( I > * D ) Ip-MlfUMUUil ( | - *ii 
J M 

•ill 1 II ( 50 « 1i0 » 30 » 40 )|S 

200 Li;.' ii*n 
n.p ii;'#n 
(l4»H3*Ii 
P5*D4*0 
M” I. 

OAMhA-h< 2/« m > 

EX-EXP ( II'.' ) 

DEX*UAMMA*2.>MHiPX 

N--1 

k •--><) , 0820558 

AIN )- A. #05 * N N II 

A<N)=AIN~1 )/l <t> N-NH 

DO 202 1-1.6 

PI “I 

AIN)=5.#D4*T##< , /5~F 1/4. ) 

202 N=N+1 

DO 203 I~1 .4 
FI = I 

A<N)=D3#4.#T##< I .5-l'J ) 

203 N=N+1 

DO 201 1-1.8 
*1=1 

A I N > =D2*3 . * r*# (1 ♦ 5-F r /2 . ) 

201 N=N+1 

DO 204 1=1.3 
F I ■- 1 

A(N) = <DEX*D3+3.*D2*EX>*T*1t<l.-Fl ) 

204 N=N+1 

DO 205 1*1,3 
F 1 = 1 

A<N)*< DEX#D5+5 . #D4#EX > * T*# < 1 . -I l ) 

205 N=N+1 
N~N~1 

P==0 

DO 215 5=1, N 
215 P=P+A< J)#D< J.M) 

I=M 

P=P+R*T# < 1 , +2 . #D#V1RB ( T > ) 

GO TO < 50 , 50 , 30 , 40 ) Is 
300 D2=D#D 
D3=D2#D 
D4=D3*D 
N*1 

k*» 0820558 
M-'I 

GAMhA=B ( 2 / , II ) 

EX-EXP ( D2#GAMMA ) 

A(N>=0.0 * N-N+l 

A ( N ) I H4*D/5 , ) *1 #* (-2 . ) # ( I. . ) M Nil 

DO 302 1*1,6 
F 1 ; -X 


AIN) -ID4/4. )*!#*•( ,/h-l 1/4. -1 



/5--E 1/4. ) 
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302 N=N+t 

DO 303 J>1 >4 
FI=J 

A<N)*(D3/3» >* T**<1.5-FI-1. >#<1.5-H> 

303 N*N+1 

DO 301 I = 1 r 8 
FI“I 

A(N)=<D2/2* >*T**<1.5-FI/2*~1* >*U.5-FI/2. > 

301 N*N+1 

DO 304 I«1 r 3 
F J«I 

A<N> = <EX/<2.*GAHMA>>*T**<1.-FI-1* >*<1.-FI> 

304 N«N+1 

DO 3053 1 = 1,3 
FI»I 

A(N)=< D2*EX/ < 2 ♦ *GAMMA ) -EX/ < 2 . *GAMMA**2 ) ) * 

1T**< 1 ♦-FI-1 • >*< 1 »— FI > 

305 N=N+1 
N=N-1 

SINT=D*R*<VIRB<T)+T*DBDT<T> > 

DO 306 1 = 1 »N 

306 SINT=SINT+B<I»M>*A<I> 

N=21 

EX=1 . 

D2=0 

DO 310 I«l»3 
FI=I 

A(N)=<EX/<2**GAMMA> >*T#*< 1 *-FI-l*)*(l*-Fl> 

310 N=N+1 

DO 311 1=1 >3 
FI«I 

A ( N > * < D2*EX/ < 2 » *GAMMA ) -EX/ ( 2 . *GAMHA**2 ) ) * 

1T**< 1 ♦-FI-1 . )*< 1 *-FI > 

311 N«N+1 
N=N-1 

DO 312 1*21 »N 

312 SINT=SINT-B<I»M)*A<t) 

P«< 9* 371658+5 »193043*AL0G<T/4« 2144) -25 ,31469* <SINT+R*AL.OG(R*T*D> ) ) 

P»P*4,0026 

I^H 

GO TD<50»50, 30,40>K 
400 KH-0 

00 TO 9 
413 PP«P 
KH=1 
D2=D*D 
D3=D*D2 
D4=D3*D 
N*1 

R= *0820558 
H«I 

QAHHA=B(27»M> 

EX*EXP < D2BGAHHA > 

A(N)=<D4*D>/5. * N-N+l 

A(N)=<D4*D/5. )*(2./T) * N=N+1 

DO 402 1=1 ,6 
FI-I 

A(N> = (D4/4. >*(T**< .75-FT/4, >-T**< .75-FI/4. >*< .75-FI/4. > > 

402 N-N+l 

DO 403 1=1 »4 
FI«1 

A<N>*<D3/3, >*<T**<1 «5-FI >-TM< 1 ,5-Fl >*( 1 »5-FI > > 


II 


403 
401 

404 


N-N+l 

DO 401 1-1,8 
FI-I 

fiijjir <D2/2 * >,<m<i * 5 ' Fi/2 *>” m<i * s - Fi/2 *> ,<i * 5 ” Fi/2 ‘ >> 

DO 404 fair 3 
FI-I 

A < N)-(EX/ < 2**0AMMA))* <T **<l.-Fl). T * 1|( a,. F i ) » (1 ^ F1 , ) 


DO 405 1*1,3 
FI-I 


*l® 2 f 1 *QAMMA ) -fcX/ < 2 * *GAMHA**2 ) ) * < Ttt < 1 » -F I ) ~T** < 1 
1*< 1 ♦ — FI > ) 

405 N-N+l 
N-N-l 


HINT-R*T*T*D*DBDT< T> *<-l«> 

DO 406 1=1 ,N 

406 HINT*HINT+B<I,M)*A<I) 

N-21 

D2«0 

ex-i. 

DO 410 1*1,3 
FI-I 

410 ^* <EX/<2 ** 6AHMA>> * <T ** <1 ‘- FI >-T**<1.-FI)*(1.-FI>) 

DO 411 1*1,3 
FI«I 

i * < i ! - Fi > > EX/ ( 2 * * GAMMA } ~ EX/ ( 2 4 * 6AMMA ** 2 > J * ( T ** ( 1 * “ F 1 } ” Tt * ( 1 

411 N-N+l 
N-N-l 

DO 412 I«21,N 

412 HINT«HINT-B<I,H)*A<I) 

P-P* 4?0026 +5 * * 93 ° 43 * <T “ 4 ‘ 2 l^>+25, 31469* < HINT+PP/D-R*T ) 

I*M 

GO T0<50,50,30»40)K 

500 D2«D*D 
D3*D2*D 
D4-D3*D 
N*1 

R- *0820558 
H*I 

GAHMA-B<27,H) 

EX-EXP <D2*GAMhA) 

SINT- V*D*R* < 2 ♦ *DBDT ( T ) +D2DBDT2 < T > *T ) 

A(N)*0*0 * N-N+l 

A<N)»<D4*D/5.)*T**(-2,)*(-l,)*(-2.) * N-N+l 

DO 502 1-1,6 
FI-I 

502 N-Um < 04/44 > * T ** ( ‘75-PI/4*-! * )*< *75-Fl/4. )*< .75-FI/4.-1. ) 

DO 503 1-1,4 
FI-I 

503 N-N+l <D3/3 * >f *®“ F * - ** 1 *5-FI )*( 1 *5-FI-l * ) 

DO 501 1-1,8 
FI-I 

501 U-N+l * G2 ^ 2 * 1 *5-F!/ 2 *-l « )*( 1 *5-FI/2* )*<1 *5-FI/2*-l* ) 

DO 504 1-1,3 
FI-I 


Ft) 


FI) 
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A<N>»<EX/<2.*0AHHA>>*T**<1*-FI-1*>*<1»-FI>*<1.-FI-1»> 

504 N-N+l 

DO 505 1-1*3 
Fl»I 

A<N>-<D2*EX/<2,*QAHMA>-EX/<2.*GAHHA**2>>* 

1T**<1.-FI-1»>*<1*-FI>*<1«-FI-1»> 

505 N-N+l 
N-N *1 

DO 506 I-1*N 

506 SINT-8INT+B< I *M>*A< I ) 

P-8XNT 

EX-1* 

D2-0 

N-21 

DO 510 1-1*3 
FI-I 

A<N>»<EX/<2.*GANHA>>*T**<1*-FI-1.>*<1.-FI>*<1.-FI-1*> 

510 N-N+l 

DO 511 1-1*3 
FI-I 

A ( H ) « ( D2*EX/ ( 2 . *OAMHA ) -EX/ < 2 ♦ *0AHHA*«2 ) ) * 

1T**<1 *-FI-l ♦ >*<l*-FI)*<l*-FI-i « ) 

511 N-H+l 
N-N-l 

DO 512 I-21*N 

512 PaP-8<I»N)*A«I) 

P-5 ♦ 1 93043* < 3 ♦ /5 . ) *4 . 0026-P* 10 1 ♦ 3278 
I-M 

00 T0<50*50*30*40)K 
20 K-2 
1*3 

IF<XD.NE*0)X«4 
T-TT 
D-DD 
00 TO 8 
30 K-3 

00 T0<33*34*35)KK 

33 D-DD 
T-TT 
KK-2 

XF<DD»QT* 17.3987)80 TO 40 
1-1 

00 TO 8 

34 PTI-P 
1-3 
KK-3 
T-TT 
D-DD 

00 TO 8 

35 PTIII-P 

38 F«<15.-T>/5. 

P«F*PTI+<1.-F>*PTIII 
IF<KH*LT»1)G0 TO 413 
FNDPHE-P 
RETURN 

40 IF<K*E0*3>Kl-3 

00 T0<41»42*43*44)K4 

41 K-4 
1-2 
K4-2 
D-DD 

XF<Kl*EQ*0>T«TT 
00 TO 8 
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4? PIID®P 
D«* 17.398/ 

IFiT.LT .5»2Q14)D«D8ATL< T >#1000. /4.0026 

K4»3 

go TO 8 

43 PIIDC«P 
1*1 
K4*4 

00 TO 8 

44 PIDC*P 

P»PIDC+<PIID-PIIDC> 

K4*l 

IF <K1 .EQ«3)G0 TO 30 
FNDI>HE«P 
50 FNDPHE*P 
RETURN 
END 

FUNCTION CPHE<D»T) 

C C8INH<P) 

C FOR HELIUM* INPUT IS DENSITY (HOLES/LITER) AND TEMPERATURE < K ) 

C OUTPUT IS JOULES/MOL-K 
C<SEE ALSO NEXT PA IE) 

P*CVHE<D»T>+<T#<DPDTHE<B*T>#*2>/< <D#*2>*DPDDHE(D»T> > >#101.325 

CPHE»P 

END 

FUNCTION SOUNHE < D * T > 

C VELOCITY OF SOUND FOR HELIUM* INPUT IS DENSITY (MOLES/LITER) 

C AND TEMPERATURE (KELVIN) * OUTPUT IS IN METRES/SEC 

S*<<CPHE<D»T>/CVHE<D»T>>*(DP0DHE<D»T>*25311.>>#*.5 
SOUNHE“S 
END 

FUNCTION VIRB(T) 

DIMENSION A<9) »V<45) 

DATA ( A»-5 *081571 004 IE-7 *-1.1 168680862E-4 *1*1 652480354E-2 » 

1 7 ♦ 4474587998E-2 *-5.31431 74768E- 1 » -9 . 575921 9306E- 1 * 

2 3 ♦ 9374414843 » -5 . 1370239224 » 2 . 0804456338 ) 

COEFFICIENTS FROM PROGRAM 5/28/70-1630 

THIS SUB PROGRAM CALCULATES THE SECOND VIRIAL COEFFICIENT FOR 
HELIUM. THE RANGE IS FROM 2 TO 1500 DEG K. INPUT IS TEMPERATURE 
IN DEGREES KELVIN AND IU3E* IF ISUE IS 0 OR NEGATIVE THE ROUTINE 
CALCULATES B FOR THE EQUATION PV*RT<1+BD+. . . * FOR OTHER VALUES OF 
IU8E* THE ROUTINE CALCULATES THE VARIANC OF B AT THE INPUT TEMP 
UNITS ARE ATM* DEG KELVIN* AND MOLES/LI TER *4/3/69-1 253 *R.D. MCCARTY 
REVISED 2/12/70-925 
1 B«0 

DO 5 I»l*9 
FI®I 

5 B®B+T**< 1 .5-FI/2. )*A(I) 

VIRB«B 
RETURN 
END 

FUNCTION DBDT <T> 

DIMENSION A(9)*V(45) 

THIS SUB PROGRAM CALCULATES THE SECOND VIRIAL COEFFICIENT FOR 
HELIUM. THE RANGE 18 FROM 2 TO 1500 DEG K. INPUT IS TEMPERATURE 
IN DEGREES KELVIN AND IU8E* IF I8UE IS 0 OR NEGATIVE THE ROUTINE 
CALCULATES B FOR THE EQUATION PV-RTC 1+BD+. ♦ . * FOR OTHER VALUES OF 
I USE* THE ROUTINE CALCULATES THE VARIANC OF B AT THE INPUT tEMP 
DAT A ( A®-5 . 081571 0041 E -7 * - 1 . 1 1 68680862E-4 *1.1 652480354E-2 * 

1 7 . 4474587998E-2 *-5.3143 1 74768E-1 * -9 . 575921 9306E-1 * 

2 3 . 9374414843 *-5.1 370239224 * 2 . 0804456338 ) 

C UNITS ARE ATM* DEG KELVIN* ANU M0LE8/L I TER * 5/28/70- 1 630 » R « D » MCCARTY 
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1 B«0 

DO s x-1,9 

f I»I 

S B»B+T#*< .5-FI/2. )*A< I >*< 1 »5-Fl/2» ) 

DBDT-B 
RE TURN 
END 

FUNCTION D2D»DT2<T> 

%SiS 3 SSS 1 ^S^-^S2:?ir?Hg 8 |^p|*5Ux« L co K rxcxc^Po S 

!!» L ^ E ES H K E ^;“^ "s E , 2 IF°Is“ “S Si H«A?IvE ^rK.NE 

CALCULATES B FOR THE EQUATION PV*RT<1+BD +. » ♦ t [? R T 2I H fSpu? L TEMp° F 
IUSE» THE ROUTINE CALCULATES THE VARIANC OF B AT THE ™PUT TEMP 
UNITS ARE ATM» DEO KELVIN, AND MQLES/L I TER , 5/28/70“ 1 4<a>0 , R ♦ D ♦ HCCART Y 

1 B«0 

DO S 1*1,9 

S B«B+T**< ,5— Fl/2,-l»)*{l,5-FI/2,)*< *5-FI/2* >*A(I> 

D2DBDT2*B 

RETURN 

FND 

n t ufr 1 VAPOR^PRESSURE FOR HELIUM IN ATMOSPHERES GIVEN A 
TEMPERATURE IN KELVIN, THE FUNCTION REPRODUCES THE i960 HELIUM 
TEMPERATURE SCALE TO ,0001 KELVIN 

MT0“i-3.»U^2«7il41. 274,75,8.-1440.7741545. 11,74.537102. 

1- 55283 .30^818* 16421,. 56504* -325212. 82840 .3,0643. 22730* 

2 M«!S®:“!S?^^i“334417.-3707.5430856.12880.4734,l. 

1^—30048. 345554. 4^532. 267434.— 5,337.558548.32311.2,4025. 

2- 33950 ,233134,1 4028 , 474003 , -5354 , 1 038947 ,1199,0301 904 , 

3 -141444342959,9,8811553384) 

T-TT 

T*T-DELT<T) 

P«-0 

IF (T-2, 1720) 10, 10,1 
1 DO 5 1*1,10 
5 P»P+C<I)*T**<2-I) 

VP *EXP<P)/»74E+4 

VPNHE-VP 

RETURN 

10 DO 15 1*1,14 
15 P«P+D<I)*T*t<2-D 
VP *EXP<P)/,74E44 
VPNHE-VP 
RETURN 
END 

FUNCTION DELKTT) 

C ADJU8T8 THE TEMPERATURE TO A HEM SCALE 
T-TT 

DELT-,0014»002*T 

RETURN 

END 

C SOLVES^THE^VAPOR^PRESSURE EOUATION FOR TEMPERATURE GIVEN APRESSURE 

C ?Se FLuId IS HELIUM AND THE UNITS ARE ATMOSPHERES AND KELVINS 
P-PP 
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lF<P.LT * .842105)00 TO 12 
T-5.0 

PCAL“VPNHE < T > 

00 TO 13 

12 PCAL». 04973/ 

IF ( AB8 < P-PCAL ) - . 0000001 #PP > 11 > 1 i * 1 

1 T-2.1720 

13 DO 10 I»l>26 
DP®DVPNHE < T ) 

DEL 9 ( PCAL-P ) /OP 
T»T-DEL 
PCAL«VPNHE<T) 

IF<ABS<P-PCAL)-.0000001*P)11>11>2 

2 IF<ABS<0EL>-«OOOOOOl*T)ll>il>lO 

10 CONTINUE 
PRINT 100>T 

11 VPTEHE«T 
RETURN 

100 FORMAT <9 TEMPERATURE ITTERATION FAILED AY T«*»>E14.?> 

END 

FUNCTION DVPNHE < TT ) 

GIVES THE DERIVATIVE OF THE VAPOR PRESSURE FOR HELIUM GIVEN A 
TEMPERATURE IN KELVIN? PRESSURE IS IN ATMOSPHERES 
DIMENSION C( 12) >D< 14) 

DATA < C»-3 . 9394635287 > 141 ♦ 27497598 > -1640 * 7741565 > 1 1974 ♦ 557102 > 

1 - 55283 ♦ 309818 >.166219 . 56504 > -3252 1 2 . 82840 > 3980*' 3 . 22750 > 

2- 277718. 069?2>83395.204183) 

DATA < D—49 .510540356 > 651 ♦ 9236417 > -3707 ♦ 5430856 > 1 2880 . 673491 > 

1 -30048 . 545554 > 49532 . 267436 > -59337 ♦ 558548 >52311. 296025 > 

2-33950 . 233134 > 16028 . 674003 > -5354 . 1038967 >1199. 0301906 > 

3 -161. 46362959>9. 8811553386) 

P*0 

T»TT-DELT(TT) 

IF<T-2.1720)10>10>1 
1 DO 5 I-l>10 
5 P«P+C<I)*T**<1-I)*<2-I) 

DVPNHE»P*VPNHE(T) 

RETURN 

10 DO 15 I«l>14 
15 P«P+D(I)tT**(l-I)*<2-I) 

DVPNHE “P4VPNHE < T ) 

RETURN 
END 

FUNCTION T*'AN8P<DD>TT) 

C HELIUM 

C THIS ROUTINE CALCULTATE8 THERMAL CONDUCTIVITY AND VI8C08ITY 
C FOR AN INPUT OF DEGREES KELVIN AND DENSITY IN MOLES PER LITER 

C THE RANGE OF TEMPERATURE IS FROM 2 TO 2000 K 

C FOR TEMPERATURE8 DELON 300 K FORMULAS OFD VINCE ARP AND GE STEWARD 

C ARE USED. FOR TEMPERATURES ABOVE 300 THE DILUTE GAS OF A CRITICAL 

C COMPILATION FROM ENGLAND IS U8ED FOR BOTH VISCOSITY AND 
C THERMAL CONDUCTIVITY AND THE EXCE8S FUNCTIONS FROM THE ROUTINES BY 

C ARP AND STEWART) THE EXCES8 FUNCTIONS ARE CALCULATED FOR TEMP8 

C ABOVE 300 K WITH THE TEMPERATURE DEPENDENCE HELD AT 300 K 

C FOR TEMPS BELOW 300 K TO 100 K THE VISCOSITY EXCESS IS CALC 

C FROM STEWARTS ROUTINE BUT THE DILUTE GA8 VALUES ARE TAKEN FROM 
C THE ENGLISH CORRELATION FOR TEMP8 BETWEEN 100 AND 110 T 
C DILUTE GAS CALCULATION IS AVARAGED 
ENTRY THERHE 
1 D*DD 
T«TT 

RH0»D*4 » 0026E-3 


tF<T.LT.300.)00 TO 5 
TH030-V I 8CX < 300 . ) * . 0078 1 736 
TH0300«CQNZ<300«> 

DEL300«DELi: < 300 . » RHO ) 

TH0®VI8CX < ' r ) * . 0078 1 736+TH0300-T HU30 
THE«TH030f *DEL300-TH0300 
TRAN8P-THJ+THE 
RETURN 

S TRAN8P f 'CQNZ < T > *DELC < T » RHO > +CR I T I C < T * RHO ) 

OUTPUT IN MU/CM. K 
RETURN 

ENTRV VISCHE 
D a DD 
T*TT 

IF<T.LT.100.)G0 TO 10 
IF<T.LT«300* >00 TO 8 
ETAO«VISCX<T> 

ET030«VISCX<300.) 

ET0300«VlSCDT<0.0»300. > 

ETE300®VI8CDT<Dt 300. )-ET0300 
TRANSPUT A0+ETE300 
: OUTPUT UNITS ARE MICROPOISE 

RETURN 

Q IF<T.LT.110.)G0 TO 9 
ETAO-UISCX(T) 

ETEB«VISCDT<D,T)-VISCDT <0»0>T> 

TRAN8P«ET AO+ETEB 
RETURN 

9 ETA1“MISCDT<0«0»100. ) 

ETA2*VI8CX<110. ) __ 

ETA0«ETA1+ < ETA2-ETA1 ) * ( T-100 . > /10 . 

TRANSP-ETA0+VI8CDT ( D » T ) -VISCDT ( 0 . 0 * T > 

RETURN 

10 TRANSP»VISCDT<D»T> 

RETURN 
END 

FUNCTION VI8CX< T > 

C CALCULATES THE DILUTE GAS VISCOSITY FOR HELIUM 

C FOR TEMPERATURES FROM 110K ™ 300K / t . 9*S W / T 

VISCX “ 196. * T ** 0.71938 * EXP( 12.451 / T 295.67 / T 

A 4*1249 ) 

RETURN 

END 

C CALCULATE * TH£ D EXCE88^THERMAL ter ) ^OF^HEL I UM^^ * “ 

C TEMPERATURE <K> AND DENSITY (M0LE8/LITER) OF HELIUM 
Q K«KZERO*EXPP<B<T>*RHO ♦ C<T>*RH0**2> 

C THIS PROGRAM RETURNS EXPF<B<T>*RHO + C<T>*RH0**2> 

1 BETTY^" * 4^7470660612 - 5?3641468153*BB + 3.4639703698*BB**2 
2 B -1.0702455443*BB**3 + 0. 1571349306*BB**4 -0.00892140047*BB**5 

21 09006708 ♦ 187. 74174B084CC - 1281 .094 7055*CC*CC 
3 +3645.2393216*CC**3 - 3V86.6937948*CC**4 
DELC-EXP < B*RHO+C*RHO*RHO ) 

RETURN 

END 

FUNCTION CONZ(TEMP) 

C KZERO IN MILL IMATT8/CM-K » T IN RELVIN9 22 JUNE 71. 

1 p at* 4? 36 11 622 157 ♦ 1 .9250 159286* ANNE - 0. 52544 1201 65*AMNE**2 
1+ 0 . 090045763885* ANNES*3 - 0.0054773874708*ANNE**4 
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1. 20654649* TL 


C 

C 


C0NZ«EXP<PAT> 

RETURN 
END 

FUNCTION viscdt<dgc,t> 

W.Q .STEWARD, 8 DATA 23 JUNE 71 

Inront^^ 1 ^ ARE KELVIN AND HOE /LITER 
OUTPUT UNITS ARE MICR0P0I8E 

TL-ALOG<T> * R-DGC*4. 0026/1000. 

ANNE" -0.1 3531 1743/TL + 1.00347841 + 

1 & c'S,*i 4 ? S645Sll,irL * Tl - +0 ‘ 0l2s2oe41A * TL **3 
BETT Y*R*<~47*5295259/TL + 87.6799309 ~42.0741589*TL 

1F8.33128289*TL*TL -0.S89252385*TL**3> 42 *° 741j89 * Tl ' 

CAROL * R*R*<547.309267/TL - 904.870586 + 431 . 40492 SSTI 
1-81.43048S4*TL*TL + S.37008433*TL**3> 431 ‘ 404V * 8 * TL 

* 3331*08630 - 1632. 19172*fL 

\,?2™?°jji 3 * TL * TL “ 20 . 29 j6367*TL**3 > 

VI3CDT«EXP<ANNE+ BETTY+ CAROL + DAQHAR) 

RETURN 
END 

FUNCTION CRITIC<TEMP»RHO> 

ANOMALY F0R THERM. CON., SCALED FROM H-2 
T IN KELVIN* REQUIRES DENSITY IN ORAHS/CC AND CP IN IDIII rn/Mni e 
THIS^EC" OF 18 SEPT 70. I HAVE US E r^C« R ??fs H E S E 0KS OF "^l/7O 

DML-RHO/O . 0040026 
IF<T .BE. 11,83) 00 TO 11 
IF (RH0.6T.0, 12) 00 TO 11 

7 CP1-CPHE<DML,T> 

8 CP2-CPHET DHL, 11.83) 

9 CRITtC-O. 0026* < CPI -CP2)/4, 0026 

10 IF<CRlTlC)ll, 12,12 

11 CRITIC-0.0 

12 RETURN 
END 

FUNCTION DSATV<TT) 

CALCULATES THE SATURATED LIQUID AND VAPOR DENSITIES rne uri him 

,NM, Srii«5!S5 e 2C [^< > .? UTPUT 13 8H/ “ <8EE ««««» 

. I '“:5»««^^t^n3”ss^. 2 ;f 32S53 ’' 2,34747o7i2> -' <o8 °“ 582 ‘ j - 

DAT A < DC-. 06964) 

DATA! TC-5. 2014) 

T-TT 
DCAL-DC 
R-<1 .-T/TC) 


4 

5 

6 


DO 1 1-1,6 
FI-I 

1 DCAL=DCAL+GV( I >*R«*(FI/3» ) 
DSATV-DCAL 

RETURN 
ENTRY DSATL 
T-TT 
DCAL-DC 
R-U.-T/TC) 

DO 2 1-1,6 
FI-I 

2 DCAL-DCAL+GL < I ) «R** < F 1/3 . ) 
D8ATV-DCAL 

RETURN 

END 


• t 

t: 
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i' CAl CULATE^THE^MFLTING PRESSURE FOR HEUIllHt 
i SuiPUT IS PRESSURE IH RTHOSPHERES 

BSmSS . S . 18*. 31 . 7»* . M . 8 IB? . 0 . W»13 ( 
T-TT 

IFIT.LE. 5*2)00 TO 7 
PHEL T-- 1 7 * 00+ i 7 . 3* 4578T8* l ♦ ASS4 1 4 
PMLL T®PMELT* . 90064/ 1 .01325 
RETURN 
7 P»0.0 

HO 9 1®1»5 
9 pap+A(I. >#T**<I~t> 

PMELT»P*9. 90665/10. 1325 
RETURN 
LNl) 

FUNCTION F NDPNE < D » T > 

C CALCULATES THE PRESSURE 
C 


INPUT TS TEMP K » 


IN ATMOSPHERES FOR NEON» 

AND TEMPERATURE OF NEON AS INPUT 


USES DENSITY < MOL/LITER) 

1 FNDPNE»PC<DtT) 

RETURN 

END 

1 »A18»A19 
D1»D 
D2«0*#2 
D3*D2*D 
D4-D34D 
D5*D4#D 
T1*T 
T2«T*T 
T3«T2*T 
T4»t3*T 
EX»EXP<A19*D2> 

EXI1»<EX-1. )/<2.*A19> 

EXI3®EX*D*D/<2.*A1V)-<EX-1 ♦ >/<2.*A19**2> 

ENTHA ®P/D-Al*T+A4*D3/3.+A5*D4/<Tl*2. ) 

l-A6*D2*T2/2 . +A96D2/T1+3 . *A10*D2/ < T2*2 . > 

2+Ai?*Dl+2»*A13*Dl/Tl+3.*A14*Dl/T2 

3+3.*A13*EXIl/T2+4.*A16*EXIl/T3 

4+3.*A17*EXI3/T2+4.*A18*EXI3/T3 

^CALCULATES^ENTHALPY CHANGE FROM 0 TO FINITE DENSITIES ALONG 

ISOTHERM ACCORDING TO DTRI4 - TU M _, re/l fTFR 

NEON INPUT IS PRESSURE IN ATMOSPHERES » DENSITY IN MOLES/LITER 

temperature in k t Output is joules per mole. 

RETURN 
END 

FUNCTION DLIQNE(T) 

FOR NEON 
IF<T~44.4>1»2»2 
A— 8.9200910 
D«35* 44150349 
C— 64. 720906 
D®57.745619 
E— 20.25961912 
R*< 1 .0-T/44 . 4)»*< i .0/3.0) 

DENR-1 ♦ 0+A*R+B*R*R+C*Rt*3+D*R**4+E*R**5 
R8ATD-10.**DENR 

DL I ONE“RSAT D*. 48 4*1 000 . 028/20 .183 
RETURN 


C 
C 

C FOR 
C AND 


Vi* 


2 R8ATD-1. 

DL IQNE-R8ATD* . 404*1000 . 020/20 , 183 
RETURN 

UNITS ARE IN 0M=MOLEB/LITER 
END 

FUNCTION DPDDNE<DD»TT> 

DIMENSION 0< 18) 

3 COMMON Al'A2'A3'A4rA5'A6'A7'A8>A9'A10'All'A12'A13'A14'AlS'A16'A17 
1 »A18»A19 

D.1-DD 

Tl-TT 

D2nDl**2 

D3»D2*D1 

D4»D3*D1 

DS-D4*D1 

T2»T1**2 

T3»T2*T1 

T4»T3*T1 

EX-EXP < A1 9SD2 ) 

EXl a EX*Dl*2.*A19 
EX3-3 » *D2*EX+D3*EX 1 
EX5=S • *D4*EX4 D5*EX 1 
Q(l)-D5*6« *A2 
0<2>-Tl*D3*4. *A3 

6(3) a D3*4»*A4 
0 < 4 ) a D4*5 * *A5/T 1 
0 ( 5 > *3 ♦ *D2*T2*A6 
0<6>-3.*D2*Tl*A7 
0<7>-3«*D2*A8 
0<9>»3.*D2*A9/T1 
0 < 9 ) “3 ♦ *D2*A10/T2 
0<10)«DI*2.*T1*AU 
0<11>-D1*2«*A12 
0( 12>-D1*2«*A13/T1 
0<13>-D1*2.*A14/T2 
0< l4)«EX3*A15/t2 
0<15) a EX3*A16/T3 
0(16)-EX3*A17/T2 
0 < 17 > a EX5*A10/T3 
P2 a Al*Tl 
DO 11 K«l»17 
11 P2*P2+0 < K ) 

DPDDNE-P2 

CALCULATES THE FIRST DERIVATIVE OF PRESSURE WITH RESPECT TO 

DEN8ITY ACCORD I N0 TO DTRI4 

RETURN 

END 

FUNCTION ENTtDvT) 

COMMON Ai'A2'A3'A4tA5'A6'A7'AS'A9»Al0»All'A12'A13'A14'A15'A16'A17 
1»A18»A19 
Dl-D 
D2-D**2 
D3»D2*D 
D4»D3*D 
D5-D49D 
T1«T 
T2-T*T 
T3«T2*T 
T4«T3*T 

EX-EXP (A19*D2) 

EXIl a EX/<A19*2» > 

EX13-EX*D*D/ ( 2 . *A19 > - < EX ) / < 2 ♦ *A19**2 ) 
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ENT -A5#B4/ < T2*4 . ) -A3*D3/3 . ~A4*D2#T 1 

J ~A7*D2/2 . +A9*D2/ < T2*2 . M-A1QBD2/T3 
2-’All*DUA13<tDl/T2+Al4t?.*Dl/T3 
34A1Q*2.*EXI1/T3+A1A*3.*EXI1/T4 
4+A17*2.*eXI3/T3+A10#3.#EXI3/T4 

8~AI*ALQG<A1*T*D> , 

CALCULATES ENTROPY AT ONE LIMIT OF INTEGRATION* ISOTHERMAL* 
ACCORDING tO DTRI4 
RETURN 
END 

FUNCTION ENTR<D*T> 

l NEON , 

COMMON A1*A2*A3*A4»AS*A6*A7»A8*A9*A10»A11*A12»A13»A14*A15*A16*A17 
1*A18*A19 
Dl*D 
D2»D»*2 
D3»D2«D 
D4®E3*D 
D5»D4*D 
T1®T 
T2«T*T 
T3-T28T 
T4»T3#T 

EX®£XP<A19#D2> 

EXI1®<EX-1*)/<2**A19> 

EXI3»EX*D*D/ < 2 ♦ *A19 ) - < EX-1 . ) / < 2 . *A19**2 ) 

ENTR »A5*D4/ < T2*4 . > -A3*D3/3 . -A6*D2*T 1 

1- A7*D2/2 . +A9#D2/ < T2*2 ♦ ) +A 104D2/T3 

2- AU#Dl+A13*Dl/T2+A14*2**Dl/T3 
3+A15»2.*EXIl/T3+A16*3**EXIl/T4 
4+A I 7*2 . *EX I 3/T3+A 1 8*3 . *EX I3/T4 
5-Al*AL0G<Al*T*D> 

C CALCULATES ENTROPY ALONG AN ISOTHERM FROM 0 TO A FINITE DENSITY 

C ACCORDING TO DTRI4 
RETURN 
END 

C FOR NEON 

SUBROUTINE DATA N 

C SETS THE CONSTANTS TO THE EQUATION OF STATE FOR NEON 

3 COMMON A1*A2»A3*A4*A5*AA»A7»A8*A9*A10*AU*A12*A13»A14*A15*A16*A17 
1*A18»A19 
Al«. 0820535 
A2®* 280101 78E-07 
A3® » 1 4244722E-05 
A4»~ ♦ 14175051E-03 
A5»- » 28070849E-04 
AA*~ 4 10185146E-06 
A7“ * 42496352E-05 
A8® ♦ 25826979E-02 
A9», 12417383 
A 1 0®- • 62382566E-0 1 
All®* 172376S9E-02 
A12®-. 16730352 
A13®~*61781598E+01 
A1 4® ♦ 8725740 2E+02 
A15® ♦ 19355993ET-02 
A16— * 84739643E+03 
A 1 7— . 65528926E-01 
A18® ♦ 234531 13E+01 

A19—.53186241E-02 
RETURN 
END 
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FUNCTION PC(DD»TT> 

C FOR NEON>P 16 IN AT MOSPHERES » F It IN MOI-ER PER l I TER 
C AN !• n 16 IN KELVIN 
DIMENSION 0(16 » 

>1 COMMON A1 r A2 r A4 » A4 » Af*» AA» A7» AS > A9 » A10 » A1 1 » At 2 / A1 3 » A1 4 > A1 ('»* A1 4 » A1 f 
1 »A16»A19 

C CALCUALTES PRESSURE ACCORDING TP P1RI4 
PI aim 
Twr 

1C' •-111**2 
D3«D2*D1 
D4«D3*Dl 
DS«D4#Dl 
D6 ki P 5*D1 
T2-T1M2 
T3>«T2#Tt 
T4«I3*Tl 
EX“EXP(A19#D2> 

Z»Di*Al*Tl 

G< 1 )=D6*A2 

G(2)=T1*D4*A3 

G<3)«D4*A4 

G<4>«D5*A5/T1 

G<S>»D3#T2*A6 

G<6)°D3*T1#A7 

G<7> B D3#A8 

3<8> B D3#A9/T1 

G(9)=D3*A10/T2 

G< 1G>®D2*T1*A1 1. 

G<U)«D2*A12 
G< 12) a D2#A13/Tl 
G<13)®D2#A14/T2 
G(.t4>*D3*EX*A15/T2 
G<15>»D3*EX#At6/T3 
G(16>«D5*EX*A17/T2 
G< 17)*D5#EX#A18/T3 
P2»Z 

DO 11 K®1»17 
11 P2*P2+G(K> 

PC»P2 

RETURN 

END 

FUNCTION FNDTNE(PtD) 

C SOLVES THE NEON EQUATION OF STATE FOR TEMPERATURE IN KELVIN 
C FOR A INPUT OF PRESSURE IN ATMOSPHERES AND DENSITY IN MOL/L 
DATA ( DTP®61 ♦ 78 ) » < DC*23 . 93 ) 

PP«P 

DD b O 

IF(DD»GT .DC >00 TO 1 

TT“44«4 

GO TO 2 

1 TT»24 . 54+ < 1 9 . 86/ ( DTP-DC ) > * ( DTP-DD ) 

2 DO 10 I«l»10 
P2*FNDPNE ( DD » TT ) 

IF < ADS < PP-P2 ) - 1 ♦ E-7*PP ) 20 > 20 » 22 
22 DP»DPDTNE(DD/TT) 

CORR» < P2-PP ) /DP 
IF(AB3(C0RR>-t*E-S>20»20»10 
10 TT»TT- CORR 
20 FNDTNE“TT 
RETURN 
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I UNI I I UN I INK P(D» I ) 

DIMENSION 0<32>»VP<V> 

CQMMUN/DA I A/G » R » GAMMA » VP » D T P » PCC » P I P » MX* f TPfTUL » II I »PU» »OCI 
UD~D 

II - I 

1M1 r .LI . I CC) 00 10 10 

i call m ss<pp»dd» m 

I IND P-f P 
RE HJKN 

to p - *v'pn\ n ) 

iiv= i iNu imp- .0001 »i r> 

DL-MND D(P+.0001»TI > 

II <0D.LE.DV. UK* . Dll . OE . Dl. >011 rtl l 
E'N l N r tOOrDOf DL »DD 

CALL PRCSSTPP.DV. IT) 

E'IND P-PP 
U-OV 
RE TURN 

100 TURMAK* THE STATE POINT YOU HAVE SPECIFIED CORRESPONDS TO A * 

1/* DENSITY IN THE LIQUID VAPOR COEXISTENCE REGION* 

LV* THE DENSITY OF THE SATURATED VAPOR IS *»FA.4»* HOLES/LITER* 

J/* THE DENSITY OF THE SATURATED LIQUID IS *,F8.4,* HOLES/1 ITER* 
4/* AND THE INPUT DENSITY IS *»F8.4»* HOLES/LI TER* 
b/* SATURATED VAPOR IS ASSURED*) 

LND 

FUNC1 ION DPDTNE < D » T > 

L GIVES THE FIRST DERIVATIVE OF PRESSURE WITH RESPECT TO 
C TEMPERATURE FOR NEON 

C INPUT PARAMETERS ARE DENSI TY (MOLES/I. I TER > AND TEMPERATURE <K > 

DIMENSION 0(18) 

3 COMMON A1 * 42»A3* A4.Ab.AAr A7r A8. A9r A10. At l*A12rA13rA14rAlf>rA1A.A17 
1 r A18r A1V 
Dl-DD 
Tl- TT 
02 = 01**2 
D3=D2*D1 
D4-DJ*D1 
D5=D4*D1 
DA=D5*D1 
T2=T1**2 
T 3- T2* 1 1 
I4=I3*T1 
EX=EXP(A1V*D2) 

G<1>=0 
0(2)= 04* A3 
0(3)=0 

0(4) =>-D5*A5/T2 
0(5)=2.*D3*T1*AA 
G( A)=D3*A7 
G< 7>=0 

G(8)=-D3*A9/T2 
G(9>=-2.*D3*A10/T3 
G(10)=D2*A11 
G< 11 )=0 

G(12)«-D2*AIJ/T2 
6< 13)=-2.*D2*A14/T3 
0< 14)=-2.*D3*EX*A15/T3 
G<15>=-3.*D3*EX*A1A/T3 
0< 1A)=-2.*DS*A17*EX/T3 
G< 17>=-3.*D5*A18*EX/T4 
P2-D1*A1 


li 


DO It K*l*17 
P2«P2FG<K> 
DPDTNE »P2 


C 

C 

C 


RE! TURN 
PND 

FUNCTION FNDTHE<P*D> 

SOLVES THE EQUATION OF 8IATF FOR HEl HIM 

FOR A TEMPERATURE! IN KELVIN) OtVf.N A PRESSURE 
AND A DENSITY IN MOLES/LITER 

DATA<DTP=36.bl4>*<DC®17.3VV> 


IN ATM 


F'P-P 

DD-D 

IF< DD*GT «DC )OU TO 1 

TT=b.201 

GO 10 2 

1 T1*2.17/4<3.024/<DTP~DC> >*<DTP-DD) 

2 DO 10 1 = 1 » 10 
P2=FNDF‘HElDD»TT> 

1 F < ADS < PP-P2 ) ~1 *E-7*PP)20*20*lt 
11 DP*DPDTHE<DD*TT> 

C0RR=<P2-PP)/DP 

IF < ADS (CORR)-l *E-5> 20*20* 10 


10 TT*TT-CORR 
20 FNDTHE=TT 
RETURN 
END 

FUNCTION F INDT <P»D) 

C SOLVE THE THIRTY TWO TERM EQUATION OF STATE FOR TEMPERATURE 
1 iS KELVIN GIVEN A PRESSURE IN ATMOSPHERES AND A DENSITY IN 
C MOLES/LITER 

COMMON/DAT A/G » R ♦ GAMMA . VP * D i P » PCC * P TP * T CC t T TP ♦ TUL » ILL » PUL » DCC 


PP=P 


DD=0 

IF(P.GE.PCC/1. 01325)00 TO 1 
TSAT=FINDTV<FP> 

DV-F T ND IK PP- . 00001 * TSAT > 

111 -f INI* IK FP-F. 0001* TSAT > 
it vUD. GT.DV.AND.DD. lt. DL >G 0 TO 30 
f l=ISAT 
GO TO 2 

1 1 1 =TCC 

2 DO 10 1*1*10 

CALL PRESS<P2*D0*TT> 
IF(ABS<PP-P2>-1 .E-/*PP >20*20* 11 
11 CALL DPDTtDP*DD.TT> 
C0RR*(P2-PP)/DP 
IF C ADSCCORR >-l *F.~b >20*20* 10 
10 TT*TT~CORR 
20 F INDT® I T 
RE I URN 


30 FIND!* TSAI 
D*DV 

PRINT 100* DV * DL » DD . 

100 FORMA! <* THE STATE POINT YOU HAVE SPECIFIED CORRESPONDS TOt 
1/* A DENSITY IN THE LIQUID VAPOR COEXISTENCE REGION* 

DENSITY OF THE SATURATED VAPOR IS**F8*4»* MOLES/LITER* 
DENSITY OF THE SATURATED LIQUID 1B*»F8.4»* MOLES/lllER* 
INPUT DENSITY 1S*»F8.4»* MOLES/L ITER* 

SATURATED VAPOR CONDITIONS ARE ASSUMED*) 


2 /* 

3/* 

4/* 

b/* 


RE I URN 
END 
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FUNCTION CPOSNF. < T > 

C CALCULATES IHE INTEGRAL OF (CP/TTDT FROM »TO»TO T FOR NEON 

L TO IS 27,092 AND CP IS AT 2ER0 PRESSURE *5/2R 
CP =8 ♦ 31434*ALQG < T /27 « 092 ) #2 • 5 
CP08NE=CP 
RETURN 

C UNITS ARE JOULES/HOLE -K 

END 

FUNCTION CPOHNFJ < 1 ) 

C CALCULATES THE INTEGRAL OF <CPT)DT FROM »TO»TO T FOR NEON»CP AT 
C ZERO PRESSURE » »T0» IS 27.0«>2 

CP =8.3t434*<T-2/.092>*2,S 
CPOHNE=CP 
RETURN 

C UNITS ARE JOULES/NOLE-K 

END 

FUNCTION FNDDNE(PP»TT) 

C SOLVES THE EQUATION OF STATE FOR NEON FOR DENSITY (HOLFS/L ) 

C GIVEN A PRESSURE IN ATH AND A TEHPERATURF IN KELVIN 
P=PP 
T*TT 
D-20. 

IF< TT.GT.44,4>G0 TO 2 
PSAT = VPNNE < T T > 

D=,OOI 

IF < PSA I »HE.PP>GO TO 2 
D“DLIQNE(T)+1 . 

2 CONTINUE 
DO 10 1=1 r50 
P2=PC(D»T) 

IF (ABS(P-P2)-I »E-6*P >11*11*1 
1 DP*DPDDNE < D » T ) 

IF < ABS( (P-P2)/DP)-1 »E-6*D)11 »11»10 

10 D=D-<P2-P)/DP 
FNDDNE=0 
RETURN 

11 FNDDNE=D 
RETURN 
END 

FUNCTION ENTH <P»D»T) 

COMMON A1 »A2»A3» A4» A5»A6»A7»A8» A9» A10»A1 1 »A12»A13*A14t A15»A14»A17 
1»A18»A19 
D1=D 
D2-D**2 
D3=D2*D 
D4=D3*D 
D5-D44D 
T1=T 
T2=T*T 
TJ=T2*T 
T4=T3*T 
EX ! *EXP< A19*D2> 

EXI1=EX/<A19*2. > 

EX I 3=EX*D*D/ < 2 . *A 19 ) - < EX > / < 2 ♦ *Ai 9**2 ) 

ENTH =P/D+A4*D3/3.FA5*D4/< Tl*2, ) 

1- A6«D2« 12/ 2 ♦ +A9*D2/T 1 F3 . *A10*D27 < T2*2 ♦ ) 

2- FA12*DlF2.*A13*Dl/TlF3.*A14*Dl/T2 
3F3.*A15*EXIi/T2+4.*Al6*EXIl/T3 
4F3.*A17*EXI3/T244.*A18*EXI3/T3 
6+A2*D5/5 . FA84D2/2 . 

C CALCUALTES ENTHALPY AT ONE LIMIT OF INTEGRATION ACCORDING 
C TO DTRI 4 t ISO THERMAL 
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